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1.  INTRODUCTION 


1.1  Background 

Ship  stmctures  are  still  designed  deterministically  according  to  working  stress 
formats.  Structural  safety  is  quantified  by  the  margin  between  the  applied  load  and  the 
capacity  of  the  structure,  which  is  measured  by  the  safety  factor.  Since  these  formats  use 
only  one  safety  factor,  they  lack  the  flexibility  to  adjust  the  prescribed  safety  margin  to 
account  for  some  of  the  factors  which  are  critical  in  design.  These  factors  include  such 
items  as  variability  in  the  strength,  loads,  modeling  uncertainties,  and  the  likelihood  of 
various  load  combinations. 

Reliability  methods  have  been  used  in  the  development  of  reliability-based  design 
formats  for  civil  engineering  and  offshore  structures,  and  they  have  matured  enough  to  be 
used  to  design  ships  more  rationally.  Reliability  methods  take  into  account  more 
information  than  their  deterministic  counterparts  in  the  analysis  and  design  of  structural 
systems.  Such  information  includes  uncertainties  in  the  strength  of  various  structural 
elements,  uncertainties  in  loads,  and  modeling  errors  in  analysis  procedures.  Probability- 
based  design  formats  are  more  flexible  and  consistent  than  working  stress  formats  because 
they  provide  uniform  safety  levels  over  various  types  of  structures.  Designers  can  use 
these  formats,  which  are  called  load  and  resistance  factor  design  (LRFD)  formats,  to 
account  for  uncertainties  that  are  not  considered  properly  by  deterministic  formats, 
without  explicitly  performing  probabilistic  analysis. 

A  commonly  used  LRFD  format  consists  of  the  requirement  that  a  factored  (reduced) 
strength  of  a  structural  component  is  larger  than  a  linear  combination  of  factored 
(magnified)  load  effects.  In  this  format,  load  effects  are  increased,  and  strength  is 
reduced,  by  multiplying  the  corresponding  characteristic  (nominal)  values  with  factors, 
which  are  called  strength  (resistance)  factors  and  load  factors,  respectively.  The 
characteristic  value  of  some  quantity  is  the  value  that  is  used  in  current  design  practice, 
and  it  is  usually  equal  to  a  certain  percentile  of  the  probability  distribution  of  that  quantity. 
The  load  and  strength  factors  are  different  for  each  type  of  load  and  strength.  The  higher 
the  uncertainty  associated  with  a  load,  the  higher  the  corresponding  load  factor.  These 
factors  are  determined  probabilistically  so  that  they  correspond  to  a  prescribed  safety 
level.  It  is  also  common  to  consider  two  types  of  performance  function,  that  correspond  to 
collapse  and  unserviceability  requirements.  The  difference  between  working  stress  and 
LRFD  formats  is  that  the  latter  use  different  safety  factors  for  each  type  of  load  and 
strength.  This  makes  it  possible  to  take  into  account  uncertainties  in  load  and  strength, 
and  to  scale  their  characteristic  values  accordingly  in  the  design  equation.  Working  stress 
formats  cannot  do  that  because  they  use  only  one  safety  factor. 

1.2  Objective 

This  project  is  a  part  of  a  long  term  effort  to  develop  a  reliability  based  method  for  the 
structural  design  of  ship  structures.  The  main  task  in  the  development  of  a  reliability- 
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based  design  method  is  the  determination  of  the  load  factors  and  the  strength  (or 
resistance)  factors.  In  order  to  achieve  this  task  it  is  necessary  to  have  a  quantitative 
measure  of  the  various  uncertainties  that  are  inherent  in  both  the  loads  and  the  strength 
models.  Earlier  SSC  projects  have  examined  the  uncertainties  related  to  loads.  The 
objective  of  this  project  is  to  develop  and  demonstrate  a  method  for  quantifying  the  bias 
and  uncertainty  in  structural  strength  algorithms  (or  computational  models)  in  order  to 
further  the  overall  goal. 
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2.  TASKS 


In  order  to  achieve  the  project's  objective  the  following  tasks  were  proposed  and  have 
now  been  accomplished. 

1)  Develop  a  methodology  for  the  modeling  and  analysis  of  uncertainties  in 
strength  parameters.  The  methodology  should  be  suitable  for  the 
development  of  a  reliability-based  design  method  for  ship  structures. 
Strength  parameters  include  both  basic  strength  variables  and  strength 
predictors.  The  uncertainties  include  bias  and  randomness  for  the  basic 
strength  variables  (e.g.,  yield  stress,  dimensions,  sizes,  etc.),  and  model 
uncertainties  in  strength  predictors  (e.g.,  buckling  strength,  plastic 
capacity,  etc.). 

2)  Identify  the  failure  modes  of  the  principal  structural  members  of  ships. 

3)  For  the  failure  modes  that  involve  modeling  uncertainty,  review  the 
availability  of  sufficient  test  data  to  demonstrate  the  method. 

4)  On  the  basis  of  this  review,  determine  which  failure  mode  is  most  suitable 
for  this  demonstration. 

5)  For  the  selected  failure  mode  (panel  compressive  collapse)  collect  data 
about  strength  parameters  and  apply  the  method  to  assess  the  uncertainties 
in  the  strength  parameters. 

6)  Determine  further  research  needs  for  uncertainty  modeling  and  analysis  of 
strength  parameters. 
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3.  METHODOLOGY  TO  ASSESS  UNCERTAINTY  IN  STRENGTH 

PARAMETERS 

The  development  of  reliability-based  design  criteria  for  marine  structures  requires  the 
assessment  of  uncertainties  in  their  predicted  strength.  The  assessed  uncertainties  in 
strength  can  be  used  for  direct  reliability-based  design  of  marine  structures.  Also,  it  is 
essential  that  the  assessed  uncertainties  are  suitable  for  the  development  of  reliability- 
based  design  formats  of  marine  structures. 

The  main  task  in  the  development  of  a  reliability-based  design  code  is  the 
determination  of  the  load  factors  and  strength  (or  resistance)  factors  in  the  process  of  code 
calibration  (Ellingwood,  et  al.  1980).  In  code  development,  it  is  common  to  normalize  the 
mean  values  of  strength  or  load  parameters  with  respect  to  their  corresponding 
characteristic  (nominal)  values  for  the  purpose  of  convenience  and  of  increasing  the  range 
of  applicability  of  the  collected  statistical  data  to  many  design  situations.  Therefore,  for  a 
given  design  situation  the  mean  values  of  the  load  effects  and  strength  can  be  computed  by 
using  these  ratios  as  multipliers.  This  approach  is  commonly  performed  as  a  part  of 
uncertainty  analysis  of  strength  and  load  effects  (Ang  and  Cornell  1974,  and  Ellingwood  et 
al.  1980). 

White  and  Ayyub  (1987a,  1987b)  demonstrated  the  development  of  reliability-based 
design  formats  of  ship  structures  for  ultimate  hull  strength  and  fatigue  failure  modes. 
Guedes  Soares  and  Moan  (1985)  demonstrated  how  to  develop  checking  equations 
(design  equations)  for  the  midship  section  under  longitudinal  bending.  They  took  into 
account  uncertainties  in  Stillwater  and  wave  bending  moments  in  calibrating  the  load  and 
strength  factors.  Committee  V2  of  ISSC  (1985)  also  presented  an  example  of  calibrating 
load  and  strength  factors  for  structural  design  of  ship  hulls. 

The  suggested  methodology  for  the  assessment  of  strength  uncertainties  of  marine 
structures  consists  of  the  following  steps: 

1.  Determination  of  nominal  (characteristic)  strength  values. 

2.  Evaluation  of  corresponding  experimental  values  (or  improved  analytical 
values). 

3.  Computation  of  stochastic  characteristics,  which  can  include  mean  value, 
coefficient  of  variation  and  distribution  type,  of  the  ratios  of  nominal  to 
experimental  values.  The  results  of  this  step  are  an  assessment  of  bias  and 
uncertainties  in  strength  parameters.  Both  objective  and  subjective  knowledge 
have  to  be  used  in  this  analysis. 

4.  Performance  of  a  parametric  analysis  of  the  stochastic  characteristics  of  the 
strength  ratios  due  to  variations  in  the  strength  parameters  and  load 
proportions. 

5.  Development  of  a  summary  of  results  in  dimensionless  spaces. 
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The  remainder  of  this  section  provides  the  needed  background  information  for  performing 
these  steps. 

3.1  Uncertainty  Types 

Uncertainties  in  structural  engineering  systems  can  be  mainly  attributed  to  ambiguity 
and  vagueness  in  defining  the  parameters  of  the  systems  and  their  relations.  The  ambiguity 
component  is  generally  due  to  non-cognitive  sources.  These  sources  include  (1)  physical 
randomness;  (2)  statistical  uncertainty  due  to  the  use  of  limited  information  to  estimate  the 
characteristics  of  these  parameters;  and  (3)  modeling  (or  prediction)  uncertainties  which 
are  due  to  simplifying  assumptions  in  analytical  and  prediction  models,  simplified  methods, 
and  idealized  representations  of  real  performances.  The  vagueness  related  uncertainty  is 
due  to  cognitive  sources  that  include  (1)  the  definition  of  certain  parameters,  e.g., 
structural  performance  (failure  or  survival),  quality  and  deterioration  of  materials,  skill  and 
experience  of  construction  workers  and  engineers,  and  conditions  of  existing  structures; 

(2)  other  human  factors;  and  (3)  defining  the  inter-relationships  among  the  parameters  of 
interest,  especially  for  complex  systems. 

Structural  engineers  and  researchers  deal  with  the  ambiguity  types  of  uncertainty  in 
predicting  the  structural  behavior  and  designing  structural  systems  using  the  theories  of 
probability  and  statistics.  Probability  distributions  are  used  to  model  system  parameters 
that  are  uncertain.  Probabilistic  structural  methods  that  include  structural  reliability 
methods,  probabilistic  engineering  mechanics,  stochastic  finite  element  methods, 
reliability-based  design  formats,  random  vibration,  and  other  methods  have  been 
developed  and  used  for  this  purpose.  In  this  treatment,  however,  a  realization  was 
established  of  the  presence  of  a  cognitive  type  of  uncertainty.  Subjective  probabilities 
have  been  used  to  deal  with  this  uncertainty  type,  that  are  based  on  mathematics  used  for 
the  frequency-type  of  probability.  Uniform  and  triangular  probability  distributions  have 
been  used  to  model  this  type  of  uncertainty  for  some  parameters.  Bayesian  techniques 
have  also  been  used  to  gain  information  about  these  parameters,  thereby  updating  the 
underlying  distributions  and  probabilities.  Regardless  of  the  nature  of  the  gained 
information,  whether  it  is  cognitive  or  non-cognitive,  the  same  mathematical  assumption 
and  tools  were  used. 

The  cognitive  types  of  uncertainty  arise  from  mind-based  abstractions  of  reality. 

These  abstractions  are,  therefore,  subjective,  and  lack  crispness.  This  vagueness  is  distinct 
from  ambiguity  in  source  and  natural  properties.  The  axioms  of  probability  and  statistics 
are  limiting  for  the  proper  modeling  and  analysis  of  this  type  and  are  not  completely 
relevant  nor  completely  applicable.  The  modeling  and  analysis  of  vagueness  type  of 
uncertainty  in  civil  engineering  systems  is  discussed  along  with  applications  of  fuzzy  set 
theory  to  such  systems  by  Ayyub  (1991)  and  Ayyub  and  Lai  (1992).  These  types  of 
uncertainty  were  not  considered  in  this  study. 
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The  sources  of  uncertainty  in  the  strength  of  a  structure  can  also  be  conveniently 
categorized  as  either  "subjective"  or  "objective"  (e.g.,  Ang  1971).  The  subjective 
uncertainties  are  those  that  result  from  the  engineer's  lack  of  knowledge  or  information 
regarding  the  physical  phenomena  associated  with  structural  failure.  These  are  usually 
manifested  in  the  form  of  imperfect  analytical  models  which  necessarily  contain 
assumptions  in  order  to  make  for  a  tractable  solution.  A  more  descriptive  title  for  these 
types  of  uncertainty  would  be  "modeling"  uncertainties. 

Some  examples  of  the  sort  of  uncertainties  which  might  be  considered  modeling 
uncertainties"  are: 

•  Uncertainties  associated  with  simple  beam  theory  in  ship  primary  bending;  i.e., 
do  plane  sections  really  remain  plane? 

•  Uncertainties  in  the  effects  of  initial  deformations  on  buckling  strength. 

•  Uncertainties  in  the  amount  of  plating  to  consider  as  acting  as  an  effective 
flange  due  to  shear  lag  effects. 

•  Uncertainties  associated  with  using  small-deflection  plate  theory. 

The  sources  of  this  uncertainty  include  our  imperfect  knowledge  of  the  failure 
mechanisms,  the  assumptions  made  in  modeling  the  failure  mechanism,  and  possible 
numerical  errors  in  the  analysis  of  the  strength.  Each  of  these  sources  of  uncertainty  needs 
to  be  considered  when  performing  an  analysis  of  the  strength  of  a  structure. 

The  objective  uncertainties  are  those  associated  with  quantities  that  can  be  measured 
and  examined.  Examples  of  such  quantities  are  yield  strength,  fracture  toughness, 
thickness,  residual  stress,  and  initial  distortion.  If  enough  data  could  be  collected  on  these 
quantities,  the  uncertainties  could  be  quantified  by  the  statistical  parameters  determined 
from  an  analysis  of  the  data.  In  some  cases,  while  there  is  a  great  deal  of  data  available,  it 
is  not  always  in  a  useful  form  (yield  strength  -  Mansour  1984,  Galambos  1978).  In  others, 
the  expense  of  collecting  data  makes  it  unlikely  that  there  will  ever  be  enough  good  quality 
data  to  perform  a  high  quality  statistical  analysis  (residual  stresses  -  Alpsten  1972).  While 
the  description  of  these  sources  of  uncertainties  as  "objective"  is  widely  accepted,  we 
believe  that  there  is  a  certain  amount  of  vagueness  contained  in  that  identification, 
particularly  to  engineers  not  versed  in  reliability  methods.  Identifying  this  type  of 
uncertainty  as  uncertainties  in  the  "basic  parameters"  more  clearly  defines  this  genre. 

It  should  be  noted  that  classifying  types  of  uncertainties  in  this  manner  is  done  merely 
for  convenience.  It  is  quite  possible  that  as  our  knowledge  regarding  some  of  the  failure 
mechanisms  improves,  things  which  have  been  identified  as  modeling  uncertainties  could 
become  uncertainties  in  the  basic  parameters. 

The  fatigue  failure  mode  is  of  great  interest,  but  for  the  most  part  it  is  dealt  with  in 
detail  design,  after  the  principal  structural  members  have  been  sized.  Several  procedures 
have  been  used  for  assessment  of  fatigue  damage  (Wirsching  1984,  Wirsching  and  Chen 
1987),  such  as  Deterministic  method.  Spectral  method,  Weibull  model,  and  Nolte-Hasford 
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model.  In  general,  the  spectral  method  is  the  most  suitable  for  marine  structures.  As  was 
demonstrated  by  Chen  and  Mavrakis  (1988),  the  spectral  method  is  more  accurate  than 
the  Weibull  model  for  the  case  of  offshore  platforms  because  its  results  are  less  sensitive 
with  respect  to  the  variability  in  the  shape  of  the  wave  spectra  compared  to  the  results  of 
the  Weibull  model.  However,  the  spectral  method  is  also  the  most  computationally 
intensive.  Moreover,  this  method  requires  the  use  of  the  weighted  sea  method  for  extreme 
value  analysis.  It  is  likely  that  the  above  conclusions  also  apply  to  ships.  In  the 
development  of  probability-based  design  guidelines  for  ships,  we  need  to  calculate  fatigue 
reliability.  Fatigue  reliability  can  be  evaluated  by  using  Munse's  model  (Munse  et  al. 

1982),  Wirsching's  model  (Wirsching  1984),  or  advanced  second  moment  methods 
(Madsen,  Skjong  and  Moghtaderi-Zadeh  1986).  A  reliability-based  design  format  for 
fatigue  was  demonstrated  by  White  and  Ayyub  (1987b). 


3.1.1.  Inherent  Uncertainty  in  Basic  Strength  Parameters 


Almost  all  strength  parameters  have  some  levels  of  inherent  uncertainty  due  to 
physical  randomness.  Therefore,  the  exact  realization  of  the  characteristics  of  these 
strength  parameters  is  not  possible.  For  practical  purposes,  these  parameters  can  be 
characterized  with  moments  and  probability  distributions.  In  developing  reliability-based 
design  codes,  it  is  necessary  to  obtain  the  central  tendency  value  (i.e.,  the  first  moment),  a 
measure  of  dispersion,  such  as  the  variance  (the  second  central  moment),  and  the 
probability  distribution  type.  Information  on  the  probabilistic  characteristics  of  the  basic 
strength  variables  of  structural  systems,  such  as  yield  stress,  endurance  limit,  modulus  of 
elasticity,  scantlings,  slenderness,  and  initial  imperfection,  can  be  obtained  by  reviewing 
the  literature.  Several  reviewed  studies  provide  information  about  the  inherent  uncertainty 
in  basic  strength  parameters,  such  as  Galambos  and  Ravindra  (1978),  Mansour  (1987), 
Stiansen  et  al.  (1979),  Daidola  and  Basar  (1981),  Guedes  Soares  and  Soreide  (1983), 
Ellingwood  et  al.  (1980),  and  Ayyub  et  al.  (1990). 

The  inherent  variability  of  a  parameter  X  can  be  expressed  in  the  form  of  the  mean 
value  ratio,  its  coefficient  of  variation,  and  its  probability  distribution  type.  The  mean 
value  is  expressed  by  the  following  notation: 


Mean  basic  strength  ratio  = 


(3-la) 


where  X  =  sample  mean  value,  and  X^  =  nominal  or  design  value.  The  coefficient  of 
variation  (COV)  of  the  ratio  is  expressed  using  the  following  notation: 

X 

Coefficient  of  variation  of  basic  strength  ratio  =  COV( - )  (3 -lb) 

The  probability  distribution  of  the  ratio  can  also  be  of  interest,  and  needs  to  be 
determined.  Frequency  analysis  or  statistical  goodness-of-fit  tests  can  be  used  for  this 
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purpose.  The  normal  or  lognormal  distributions  are  commonly  used  to  model  strength 
parameters.  Other  distributions  were  also  used  such  as  the  Rayleigh  distribution. 

Equations  3- la  and  3- lb  provide  expressions  for  the  mean  and  coefficient  of  variation 
of  a  strength  parameter  ratio.  The  benefit  of  expressing  the  results  in  the  form  of  ratios  is 
in  providing  the  results  in  dimensionless  multipliers  that  can  be  used  for  wide  range  of 
nominal  values.  Also,  these  expressions  lend  themselves  for  comparative  analysis  for 
assessing  the  level  of  inherent  uncertainties  in  different  strength  parameters. 

3.1.2.  Prediction  Uncertainty 

As  was  discussed  in  previous  sections,  the  prediction  uncertainty  is  due  to  bias  and 
variability  in  the  predictions  of  analytical  models  due  to  their  limitations,  assumptions  and 
model  simplifications.  This  uncertainty  can  be  quantified  by  determining  the  nominal  (or 
characteristic)  strength  values  of  some  strength  measures  of  interest.  Then  the 
corresponding  experimental  values  (or  improved  analytical  values)  need  to  be  evaluated. 
Computations  of  stochastic  characteristics,  which  include  mean  value,  coefficient  of 
variation  and  distribution  type,  of  the  ratios  of  experimental  to  nominal  values  are  then 
performed.  The  results  are  assessments  of  bias  and  random  uncertainties  in  strength 
parameters.  Both  objective  and  subjective  knowledge  have  to  be  used  in  this  analysis. 

The  performance  of  a  parametric  analysis  of  the  stochastic  characteristics  of  the  strength 
ratios  due  to  variations  in  the  strength  parameters  and  load  proportions  can  provide  an 
assessment  of  the  stability  of  the  results  and  their  ranges.  Then,  the  results  in 
dimensionless  spaces  need  to  be  summarized. 

According  to  measurement  theory  (Ang  and  Tang  1984),  the  error  in  an  estimated 
quantity  can  be  considered  to  consist  of  two  types,  the  systematic  (bias)  component  and 
the  random  (variability)  component. 

The  bias  component  (B)  can  be  expressed  as 

B  =  ^  •  (3-2) 

Xq 

where  =  the  actual  strength,  and  Xq  =  the  design  value  of  the  strength  parameter. 
Sometimes,  for  convenience,  the  bias  (B)  is  broken  down  into  three  components,  the 
actual  to  experimental  bias  (Bj),  experimental  to  predicted  bias  (B2),  and  the  predicted  to 
design  bias  (B2).  The  actual  to  experimental  bias  (Bj)  is 


The  experimental  to  predicted  bias  (B2)  is  given  by 
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(3-4) 


R  - 

where  Xg  =  the  experimental  value  of  the  strength  parameter,  and  Xp  =  the  predicted 
value  of  the  strength  parameter.  The  predicted  to  design  bias  (B3)  is  given  by 


= 


Therefore,  the  total  bias  (B)  given  by  Eq.  3-2  is  the  product  of  Bj  to  B3  as  follows; 

B  =  B1B2B3  (3-6) 

The  coefficient  of  variation  of  the  bias,  COV(B),  can  be  considered  to  be  a  measure  of 
the  random  component  of  uncertainty.  Therefore  it  is  of  interest  and  needs  to  be  assessed. 
In  cases  where  the  prediction  bias  includes  several  sources,  the  total  variability  (COVj) 
can  be  computed  using  a  first  order  approximation: 


COVx  =  JCOV?  +  COV2  +  ...  +  COVn 


(3-7a) 


where  COVi  =  the  ith  variability  source.  For  example  the  coefficient  of  variation  of  the 
total  bias  COV(B)  can  be  computed  as; 


COV(B)=  VCOV^(Bi)  +  C0V^(B2)  +  C0V^(B3) 
The  above  equations  are  true  if  the  COV's  are  small. 


(3-7b) 


The  probability  distribution  type  of  the  corrected  strength  measure  by  the  bias  factor  is 
also  of  interest  and  needs  to  be  determined. 


The  statistics  of  predicted  strength  can  be  calculated  using  approximate  methods 
(second  moment)  or  exact  methods  (Monte  Carlo  simulation)  (Ang  and  Tang  1975). 
Monte  Carlo  simulation  methods  are  used  to  assess  the  random  bias  in  cases  where  the 
predicted  strength  (Xp)  is  given  in  terms  of  a  function,  such  as 

Xp=  g(Xi,  X2, ...,  Xn)  (3-8) 

where  Xj  (i  =  1,2, ...  n)  =  basic  random  variable.  The  function  can  be  explicit  as  given  by 
Eq.  3-8,  or  implicit  in  the  form  of  an  elaborate  computational  procedure  that  is  possibly 
evaluated  by  computers.  The  mean  predicted  strength  value  can  be  obtained  using  Monte 
Carlo  simulation  as  follows: 


Xp=-EXp, 
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where  Xpj  =  the  predicted  strength  in  the  to  simulation  cycle,  and  N  =  the  total  number  of 
simulation  cycles.  The  coefficient  of  variation  of  Xp  can  be  computed  as 


COV(Xp)= 


^N-1 


Xp 


(3-10) 


The  statistical  error  in  the  estimated  mean  values  (Eq.  3-9)  is  given  by 


COV(Xp)  = 


(3-11) 


The  distribution  type  can  be  determined  using  frequency  analysis,  or  statistical  goodness- 
Of-fit  tests. 


The  mean  value  and  coefficient  of  variation  of  Xp  can  be  approximately  evaluated 
using  a  first-order  Taylor-series  expansion  as  follows: 

Pp  «  g(Pl,  P2.  •••>  l^n)  (3-12a) 

and 


COV(Xp) 


n  n 


TTpii  ^  Ox-ax- 

^  dXi  dXj 


(g(Xi,  X2,  ...,  X„))^ 


(3-12b) 


where  Pi=  mean  value  of  Xj,  qj  =  the  correlation  coefficient  between  Xj  and  Xj,  and  s  = 
standard  deviation.  The  partial  derivetives  are  evaluated  at  the  mean  value. 

3.1.3.  Statistical  Uncertainty 

The  selection  of  a  method  for  quantifying  uncertainty  depends  on  the  available 
information  and  its  nature.  These  methods  are  generally  based  on  statistical  concepts, 
probability  and  Bayesian  approaches. 

In  cases  where  small  sample  sizes  are  used  to  assess  a  parameter,  the  statistical 
uncertainty  in  the  parameter  needs  to  be  quantified.  For  example,  the  sample  mean  and 
standard  deviation  can  be  computed  for  a  sample  of  size  n,  such  as,  xj,  X2, ...,  x^,  as 
follows: 
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sample  mean,  X 


1  “ 

=  rS-'i 


(3-13) 


Sample  standard  deviation,  sx 


in  _  » 
n-li=i 


(3-14) 


The  estimated  mean  (X)  is  uncertain  due  to  the  limited  sample  size  and  its  dependence  on 
the  sample.  This  statistical  uncertainty  in  the  assessed  mean  can  be  determined  as 


(3-15) 


In  equation  (3-15)  we  can  use  the  sample  standard  deviation  instead  of  the  population 
standard  deviation  if  the  latter  is  not  available.  This  coefficient  of  variation  of  the  sample 
mean  accounts  only  for  the  statistical  uncertainty.  It  does  not  account  for  other  types  of 
uncertainty  that  are  inherent,  for  example,  in  the  sampling  or  testing  programs.  In  this 
study  statistical  uncertainty  was  not  considered. 

3.1.4.  Subjective  Assessment  of  Statistical  Parameters 

For  some  parameters,  limited  information  about  them  might  require  an  analyst  to  use 
judgment  to  quantify  uncertainty.  For  example,  if  only  a  range  on  the  parameter  is 
assessed  (L  =  lower  limit  of  X ,  U  =  upper  limit  of  X),  then  the  mean  value  and  coefficient 
of  variation  of  X  can  be  determined  as 


Mean  value,  p  =  —  (L  +  U) 


(3-16) 


Coefficient  of  variation,  COV(X) 


1  rU-L 
V3lu  +  L 


(3-17) 


If  a  triangular  distribution  is  assumed  over  the  range  (L,U)  with  a  mode  at  M,  then  the 
mean  and  standard  deviation  are  given  by 


Mean  value,  p  = 


L  +  M  +  U 


(3-18) 


Standard  deviation,  ax  ~ 


+M^  +U^  -LU-LM-MU 


(3-19) 
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For  example,  if  +  two  standard  deviations  of  a  normal  probability  distribution  are  assumed 
to  be  represented  by  the  range  (L,U),  then  the  mean  and  coefficient  of  variation  for  the 
normal  probability  distribution  are  given  by 


Mean  value,  p  = 


L  +  U 
2 


(3-20) 


Coefficient  of  variation,  COV(X) 


1  U-L 
2U  +  L 


(3-21) 


3.2.  Bayesian  Techniques 

Engineers  commonly  need  to  solve  a  problem,  and  they  must  make  decisions  based  on 
limited  information  about  one  or  more  of  the  parameters  of  the  problem.  The  types  of 
information  available  to  them  can  be 

1 .  objective  information  based  on  experimental  results,  or  observations; 

2.  subjective  information  based  on  experience,  intuition,  other  previous  problems 
that  are  similar  to  the  one  under  consideration,  or  the  physics  of  the  problem. 


The  first  type  of  information  can  be  dealt  with  using  the  theories  of  probability  and 
statistics  as  was  described  in  the  previous  sections.  In  this  type,  probability  is  interpreted 
as  the  frequency  of  occurrence  assuming  sufficient  repetitions  of  the  problem,  its 
outcomes,  and  parameters,  as  a  basis  of  the  information.  The  second  type  of  information 
is  subjective  and  can  depend  on  the  engineer  or  analyst  studying  the  problem.  In  this  type, 
uncertainty  exists,  and  needs  to  be  dealt  with  using  probabilities.  However,  the  definition 
of  probability  is  not  the  same  as  the  first  type,  it  is  viewed  herein  as  a  subjective 
probability  that  reflects  the  state  of  knowledge  of  the  engineer  or  the  analyst. 


It  is  common  in  engineering  to  encounter  problems  with  both  objective  and  subjective 
types  of  information.  In  these  cases,  it  is  desirable  to  utilize  both  types  of  information  to 
obtain  solutions  or  make  decisions.  The  subjective  probabilities  are  assumed  to  constitute 
a  prior  knowledge  about  a  parameter,  with  gained  objective  information  (or  probabilities). 
Combining  the  two  types  produces  posterior  knowledge.  The  combination  is  performed 
based  on  Bayes'  theorem. 
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If  Ai,  A2,  An  represent  the  prior  (subjective)  information,  or  a  partition  of  a 
sample  space  S,  and  E  cz  S  represents  the  objective  information  (or  arbitrary  event)  as 
shown  in  Figure  3. 1,  the  theorem  of  total  probability  states  that 

P(E)  =  P(Ai)  P(ElAi)  +  P(A2)  P(E1A2)  +  ...  +  P(An)  P(ElAn)  (3-22) 

This  theorem  is  very  important  in  computing  the  probability  of  the  event  E,  especially  in 
practical  cases  where  the  probability  cannot  be  computed  directly,  but,  the  probabilities  of 
the  partitioning  events  and  the  conditional  probabilities  can  be  computed. 


Bayes'  theorem  is  based  on  the  same  conditions  of  partitioning  and  events  as  the 
theorem  of  total  probability  and  is  very  useful  in  computing  the  posterior  (or  reverse) 
probability  of  the  type  P(Ai|E),  for  i  =  1,  2, ...,  n.  The  posterior  probability  can  be 
computed  as  follows: 


P(Ai|E)  = 


_ P(Af)P(E|Ai) _ 

P(Ai)P(E|  Ai)  +  P(A2)P(E|  A2)+.  . .  +P(A„)P(E|  A^) 


(3-23) 


The  denominator  of  this  equation  is  P(E),  which  is  based  on  the  theorem  of  total 
probability.  According  to  Eq.  3-23,  the  prior  knowledge,  P(Ai),  is  updated  using  the 
objective  information,  P(E),  to  obtain  the  posterior  knowledge,  P(AilE).  Additional 
information  on  Bayesian  techniques  is  provided  by  Ang  and  Tang  (1975). 


3.2.1.  Discrete  Case 

For  an  unknown  parameter  0,  a  prior  distribution  for  the  parameters  can  be 
subjectively  determined,  and  expressed  using  a  probability  mass  function  as 

P0(ei)  =  P(0=0i)  for  i  =  1. 2, ...,  n  (3-24a) 

or,  in  an  abbreviated  form,  as 

P0(0i)  =  P(ei)  for  i  =  1,  2, ...,  n  (3-24b) 


Therefore,  the  parameter  0  is  assumed  to  take  n  discrete  values  with  probabilities  given  by 
Eqs.  3-24.  The  distribution  of  0  reflects  the  uncertainty  in  this  parameter  including  its 
randomness.  It  is  possible  to  have  a  parameter  that  is  not  random,  but  uncertain,  therefore 
requiring  its  description  through  a  distribution  as  given  by  Eqs.  3-24. 
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Now  assume  that  new  (objective)  information  8  was  obtained.  Using  Eq.  3-23,  the 
posterior  distribution  of  the  parameter  can  be  obtained  as 


P(0i|e)  = 


_ P(9i)P(s|9i) _ 

P(9i)  P(8|9i)+P(92)  P(8|92)+  ...+P(9n)  P(8|9„) 


(3-25a) 


where  P(9i|e)  =  the  conditional  probability  of  9i  given  e,  or  the  posterior  probability  for  9}; 
P(9i)  =  prior  probability  as  given  by  Eqs.  3-24;  and  P(e|9i)  =  the  probability  of  obtaining 
the  new  information  (e)  given  a  certain  value  (9i)  for  the  parameter.  The  following 
notation  for  the  posterior  distribution  is  also  common: 


P'(9i)  = 


_ P(9i)P(8|ei) _ 

P(9i)  P(8|9i)  +  P(92)  P(8|92)+  ...+P(9n)  P(8|9„) 


(3-25b) 


where  P'(9i)  =  the  conditional  probability  of  9}  given  e,  or  the  posterior  probability  for  qj. 

Using  the  prior  distribution  of  the  parameter  0  given  by  Eqs.  3-24,  the  expected  value 
of  the  parameter  can  be  computed  as 

E(0)  =  SeiPCSi)  (3-26) 

i=l 

Based  on  the  posterior  distribution,  the  expected  value  of  0  can  be  computed  as 

E(0|e)=  |;6iP'(ei)  (3-27) 

i=l 


In  many  engineering  problems,  the  parameter  0  can  be  used  to  define  a  probability 
distribution  of  a  random  variable  X.  The  probability  distribution  of  X  can  be  either  for  a 
discrete  random  variable  in  the  form  of  a  probability  mass  function,  Px(x),  or  for  a 
continuous  random  variable  in  the  form  of  a  density  function,  fx(x).  The  Bayesian 
estimation  of  the  parameter  can  be  used  to  compute  Bayesian  probabilities  that  are 
obtained  using  the  gained  information  about  the  parameters.  For  example,  the  probability 
that  X  is  less  than  some  value  Xq  can  be  computed  using  the  prior  distribution  as 

P(X<Xo)  =  f;P(X<Xo|9i)P(9i)  (3-28) 

i=l 

or 
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(3-29) 


FxW  =  2;Fx(xoiei)P(0i) 

i=l 

where  Fx(x)  =  the  cumulative  distribution  function  of  X  evaluated  at  Xq.  Using  the 
posterior  distribution  results  in  the  following  expressions: 

P(X<Xo)  =  SP(X<Xo|0i)F(ei)  (3-30) 

i=l 

or 

Fx(>^)  =  i;Fx(x„iei)P'(ei)  (3-31) 

i=l 

3.2.2.  Continuous  Case 

For  an  unknown  parameter  0,  a  prior  distribution  for  the  parameters  can  be 
subjectively  determined,  and  expressed  using  a  probability  density  function  f0(9).  The 
parameter  0  is  assumed  to  be  continuous  with  probabilities  that  can  be  computed  based 
on  its  density  fimction.  Again,  the  distribution  of  0  reflects  the  uncertainty  in  this 
parameter  including  its  randomness. 

Now  assume  that  new  (objective)  information  s  was  obtained.  Using  Eq.  3-23,  the 
posterior  distribution  for  the  parameter  can  be  obtained  as 


f0(0|E)  = 


fQ(9)  P(e|e) 

00 

Jp(e|e)fe(0)de 


(3-32) 


where  f0(0)  =  the  prior  density  function  of  0;  f0(0|e)  =  the  posterior  density  function  of 
0;  and  P(e|0)  =  the  probability  of  obtaining  the  new  information  (e)  given  a  certain  value 
for  the  parameter  (0).  The  probability  P(e|0)  is  called  the  likelihood  function  L(0).  The 
following  notations  for  the  posterior  distribution  is  also  common: 
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f  0(0)  = 


(3-33) 


f©(e)L(9) 

00 

jL(0)f0(0)d0 
—00 

where  f0(0)  =  the  conditional  density  function  of  0  given  e,  or  the  posterior  density 
function  of  0. 

Using  the  prior  density  function  of  the  parameter  0,  the  expected  value  of  the 
parameter  can  be  computed  as 

00 

E(0)=  J0f0(0)d0  (3-34) 

-00 

Based  on  the  posterior  distribution,  the  expected  value  of  0  can  be  computed  as 

00 

E(0|8)=  J0f'@(0)d0  (3-35) 

—00 


In  many  engineering  problems,  the  parameter  0  can  be  used  to  define  a  probability 
distribution  of  a  random  variable  X.  The  probability  distribution  of  X  can  be  either  for  a 
discrete  random  variable  in  the  form  of  a  probability  mass  function,  Px(x),  or  for  a 
continuous  random  variable  in  the  form  of  a  density  function,  fx(x).  The  Bayesian 
estimation  of  the  parameter  can  be  used  to  compute  Bayesian  probabilities  that  are 
obtained  with  the  gained  information  about  the  parameters.  For  example,  the  probability 
that  X  is  less  than  some  value  Xq  can  be  computed  using  the  prior  distribution  as 

00 

P(X<Xo)  =  Jp(X<xJ0)f®(0)d0  (3-36) 

-00 

or 

00 

Fx(xo)  =  jFx(Xo|0)f©(0)d0  (3-37) 

-00 
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where  Fx(xo)  =  the  cumulative  distribution  function  of  X  evaluated  at  Xq.  Using  the 
posterior  distribution  results  in  the  following  expression: 

00 

P(X<Xo)  =  Jp(X<Xo|e)f'@(0)d0  (3-38) 

-00 


or 


00 

Fx(xo)  =  jFx(Xo|0)f'@(0)d0 

-X 


(3-39) 


3.2.3  Bayesian  Statistics  -  Mean  Value  with  Known  Variance 

The  Bayesian  methods  that  were  developed  in  the  previous  sections  can  be  used  in  the 
statistical  analysis  of  data.  In  this  section,  two  cases  are  considered  to  illustrate  their  use 
in  statistical  analysis.  The  first  case  deals  with  a  random  variable  X  that  is  normally 
distributed  with  a  known  standard  deviation.  The  mean  value  of  the  random  variable  is  of 
interest,  and  is  estimated  using  Bayesian  methods.  In  the  second  case,  the  random  variable 
X  is  also  normally  distributed,  but  its  standard  deviation  is  unknown.  In  this  case,  both  the 
mean  value  and  the  variance  of  the  random  variable  are  of  interest,  and  are  estimated  using 
Bayesian  methods. 

A  random  variable  X  is  considered  to  be  normally  distributed  with  a  known  variance 
a^.  The  mean  value  of  the  random  variable  is  of  interest,  and  is  unknown.  The  prior 
distribution  of  the  unknown  mean  (|i)  is  normal  with  a  mean  value  and  variance 

Pq,  and  Oq,  respectively.  New  (objective)  information  was  obtained  by  a  sample  of  size 

n.  The  mean  value  based  on  the  sample  is  X.  We  are  interested  in  determining  the 
posterior  distribution  of  the  mean.  Using  Eq.  3-33,  the  following  expression  can  be 
established: 


f(p)  = 


f(^)  L(i^) 

CO 

jL(|i)f(n)  dg 


-CO 


(3-40) 
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Jiere  f(|A)  =  the  prior  density  function  of  p,  which  is  normal  with  mean  and  variance  of 

|io,  and  ol,  respectively,  (i.e.,  N(Po,ao));  f(|i)  =  the  posterior  density  function  of  the 
unknown  mean  p;  and  L(p)  =  the  likelihood  function  for  the  sample  of  size  n.  The 
likelihood  function  can  be  computed  as  the  product  of  n  values  of  the  density  function  of 
the  normal  distribution  with  a  mean  p  and  standard  deviation  s,  each  evaluated  at  a 
sampled  value  xj.  The  product  can  be  expressed  as 


L(p)  = 


,,n/2  n 
(271)  a 


exp 


1 


X;  -p 


^2 


V  a  J 


(3-41) 


It  can  be  shown  that  by  substituting  Eq.  3-41  in  Eq.  3-40,  the  resulting  f(p)  is  normally 
distributed  with  the  following  mean  value  and  variance,  respectively: 


nXa^+Poa^ 

2  2 
nao  +a 


(3-42a) 


i-> 


2_2 


aoO 


nao 


(3-42b) 


The  resulting  p',  and  s'  are  the  posterior  mean  and  standard  deviation  of  the  unknown 
mean  value  p .  Using  the  normal  posterior  distribution,  any  Bayesian  probabilities  of 
interest  for  the  random  variable  X  can  be  computed. 


The  prior  and  posterior  mean  values  and  variances  can  also  be  used  in  other  aspects  of 
statistical  analysis  such  as  confidence  intervals,  and  hypothesis  testing.  For  example,  they 
can  be  used  to  establish  the  following  prior  confidence  interval  on  the  mean: 

^o-^a/2^0  ^  ^  flo+Za/2<^o  (3-43a) 

Also,  they  can  be  used  to  establish  the  following  posterior  confidence  interval: 


1^'  -^a/2^'  ^  ^  1^'  +2a/2<^'  (3-43b) 

where  (1-a)  is  the  confidence  level.  In  a  similar  approach  prior  hypothesis  testing,  and 
posterior  hypothesis  testing  can  be  performed. 
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3.2.4.  Bayesian  Statistics  -  Mean  Value  with  Unknown  Variance 

The  random  variable  X  in  this  case  is  considered  to  be  normally  distributed  with  a 
unknown  mean  value  (n),  and  unknown  variance  (o^).  Both  the  mean  value  and  variance 
of  the  random  variable  are  of  interest,  and  are  unknown.  The  prior  joint  distribution  of  the 
unknown  mean  (p),  and  unknown  variance  (o^)  is  assumed  to  be  normal-Gamma  which  is 
defined  as  the  product  of  a  normal  distribution  for  the  mean  (p),  and  a  Gamma 
distribution  for  the  variance  (o^).  The  prior  information  about  the  mean  and  variance  is 

based  on  a  sample  of  size  N  with  sample  mean  and  variance  of  Xq  and  Sq,  respectively. 
New  (objective)  information  was  obtained  by  a  sample  of  size  n.  The  mean  value  and 
variance  based  on  the  sample  are  X  and  S^,  respectively.  We  are  interested  in 
determining  the  posterior  distribution  of  the  mean  and  variance.  Using  Eq.  3-33,  it  can  be 
shown  that  the  posterior  distribution  is  also  a  normal-Gamma. 


The  posterior  mean  (X'),  and  posterior  variance  (S'^)  can  be  shown  to  be 

-  ^  NXp+nX 
n' 


(3-44) 


and 


(N  -  l)So  +  NX^  +  (n  -  1)S^  +  nX^  -  n'  X  ^ 
n'-l 


(3-45) 


where 

n'=  N  +  n  (3-46) 

The  resulting  values  from  Eqs.  3-44  to  3-46  are  the  posterior  mean  and  standard  deviation 
of  the  unknown  mean  and  variance. 


3.3.  Relative  Importance  of  Strength  Parameters 

The  relative  importance  of  different  strength  parameters  needs  to  be  assessed  for  the 
purpose  of  allocating  resources  in  quantifying  unknown  uncertainties,  quality  control, 
development  of  design  changes,  and  reliability  improvement.  Parametric  analysis, 
sensitivity  factors,  and  weighted  sensitivity  factors  can  be  used  for  that  purpose.  They 
offer  some  benefits  in  certain  aspects  towards  this  objective. 
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3.3.1.  Parametric  Analysis 


Monte  Carlo  simulation  methods  can  be  used  to  assess  the  mean  value,  coefficient  of 
variation  and  distribution  type  of  the  predicted  strength  (Xp)  which  is  according  to  a 
function  of  the  type 

Xp=  g(Xi,X2....,X„)  (3-47) 

where  Xj  (i  =  1,2, ...  n)  =  basic  random  variable.  The  function  can  be  explicit  as  given  by 
Eq.  3-47,  or  implicit  in  the  form  of  an  elaborate  computational  procedure  that  is  possibly 
evaluated  by  computers.  The  mean  predicted  strength  value  can  be  obtained  using  Monte 
Carlo  simulation  according  to  Eq.  3-9.  The  coefficient  of  variation  of  Xp  can  be 
computed  using  Eq.  3-10.  The  statistical  error  in  the  estimated  mean  value  is  given  by  Eq. 
3-11.  The  distribution  type  can  be  determined  using  frequency  analysis,  or  statistical 
goodness-of-fit  tests. 

The  objective  of  parametric  analysis  is  to  investigate  the  effect  of  perturbing  the  mean 
value,  coefficient  of  variation,  or  changing  the  distribution  type  of  only  one  of  the  basic 
random  variables  on  the  mean  value,  coefficient  of  variation,  or  distribution  type  of  the 
predicted  strength  Xp. .  The  analysis  is  typically  repeated  for  every  basic  random  variables. 
The  effects  of  Xp  can  be  normalized  into  some  convenient  dimensionless  quantity,  and 
then  ranked  according  to  their  gravity. 

In  this  study,  two  methods  of  parametric  analysis  were  developed.  The  first  method 
finds  the  uncertainty  effect  of  each  basic  random  variable  when  only  the  corresponding 
variable  is  random  and  the  other  variables  are  fixed  at  values  equal  to  their  respective 
means.  For  this  purpose,  Monte-Carlo  simulation  for  each  random  variable,  while  keeping 
the  other  variables  fixed  at  their  mean  values  needs  to  be  performed.  Then  the  resulting 
mean  and  coefficient  of  variation  of  the  predicted  strength  (Xp)  can  be  calculated.  The 
results  can  be  expressed  using  the  following  parametric  coefficients  (Pij). 


COV(Xp)i  (3.48) 

|;(COV(Xp)if 
i=l 

where  Pii  is  the  parametric  coefficient  of  the  ith  random  variable  using  method  1,  and 
COV(Xp)i  is  the  coefficient  of  variation  of  the  predicted  strength  due  to  uncertainty  in  the 
ith  random  variable. 


The  second  method  calculates  the  derivatives  of  the  coefficient  of  variation  of  the 
predicted  strength  with  respect  to  the  coefficients  of  variation  of  the  basic  random 
variables.  The  sensitivities  are  found  in  reference  to  a  "base"  strength  calculated  using  the 
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statistics  of  the  basic  random  variables.  In  this  case,  the  parametric  coefficient  (P2i )  of 
the  uncertainty  in  the  ith  basic  random  variable  is  calculated  using  the  following  equation; 


ACOV(Xp  )i 
Fn 

2;(icov(Xp)i 


(3-49) 


where  P2i  is  the  parametric  coefficient  of  the  ith  random  variable  using  method  2;  and 
DCOV(Xp)j  is  the  change  in  the  coefficient  of  variation  of  the  predicted  strength  due  to  a 
change  in  the  coefficient  of  variation  of  the  ith  basic  random  variable  by  some  percent. 
The  percent  change  needs  to  be  the  same  for  all  the  basic  random  variables. 


3.3.2.  Sensitivity  Coefficients 

The  sensitivity  coefficients  are  defined  in  this  study  as  the  normalized  change  in 
predicted  strength  due  to  a  change  in  the  moments  of  the  basic  random  variables.  In  this 
method,  all  the  random  variables  are  generated  in  the  simulation  of  the  predicted  strength. 
The  following  four  cases  were  developed: 

1 .  The  sensitivity  coefficient  for  the  mean  predicted  strength  due  to  a  change  in 
the  mean  value  of  a  basic  random  variable  (C^^j)  is  computed  as 


AXpj/Xp 
AXj  /Xj 


"  '^AXpi/Xj 


AXj/Xi 


(3-50a) 
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where  AXj  =  change  in  the  mean  value  of  the  ith  basic  random  variable;  Xp  - 

the  mean  value  of  the  predicted  strength;  and  AXpj  =  change  in  the  mean  value 
of  the  predicted  strength  due  to  the  change  in  the  mean  value  of  the  ith  basic 
random  variable.  The  means  of  the  basic  random  variables  are  perturbed  one 
random  variable  at  a  time,  while  keeping  the  means  of  the  remaining  random 
variables  at  their  non-perturbed  means.  Their  means  are  also  perturbed  using 
the  same  percent  of  their  respective  means.  The  selection  of  the  percent 
change  should  be  based  on  realistic  possible  levels  for  the  means.  The  percent 

change  and  Xp  in  Eq.  3-50a  cancel  out  from  it  to  produce  the  following 
equation: 


C 


AX  Pi 


i=l 


(3-50b) 


2  The  sensitivity  coefficient  for  the  coefficient  of  variation  of  predicted  strength 
due  to  a  change  in  the  mean  value  of  a  basic  random  variable  (Copi)  is 
computed  as 


ACOV(Xp)i/COV(Xp) 

AXj/Xj 


n  f 


ACOV(Xp)i/COV(Xp) 

AXj/Xi 


^2 


(3-5  la) 


where  COV(Xp )  =  the  coefficient  of  variation  of  the  predicted  strength;  and 
ACOV(Xp)i  =  change  in  the  coefficient  of  variation  of  the  predicted  strength 
due  to  the  change  in  the  mean  value  of  the  ith  basic  random  variable.  The 
means  of  the  basic  random  variables  are  perturbed  one  random  variable  at  a 
time,  while  keeping  the  means  of  the  remaining  random  variables  at  their  non- 
perturbed  means.  The  means  are  also  perturbed  using  the  same  percent  of  their 

respective  means.  The  percent  change  and  COV(Xp)  in  Eq.  3-51a  cancel  out 
from  it  to  produce  the  following  equation; 
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ACOV(Xp  )i 

“  jH - - 

IX(ACOV(Xp)if 


(3-5  lb) 


3.  The  sensitivity  coefficient  for  the  mean  predicted  strength  due  to  a  change  in 
the  coefficient  of  variation  of  a  basic  random  variable  (G|iai)  is  computed  as 


where  ACOV(Xi)  =  change  in  the  coefficient  of  variation  of  the  ith  basic 

random  variable;  Xp  =  the  mean  value  of  the  predicted  strength;  and  AXp,  = 
change  in  the  mean  value  of  the  predicted  strength  due  to  the  change  in  the 
coefficient  of  variation  of  the  ith  basic  random  variable.  The  coefficients  of 
variation  of  the  basic  random  variables  are  perturbed  one  random  variable  at  a 
time,  while  keeping  the  COV's  of  the  remaining  random  variables  at  their  non- 
perturbed  coefficients.  The  COV's  are  also  perturbed  using  the  same  percent 
of  their  respective  non-perturbed  values.  The  coefficients  of  variation  can  be 
perturbed  by  perturbing  the  corresponding  standard  deviations  of  the  basic 

random  variables.  The  percent  change  and  Xp  in  Eq.  3-52a  cancel  out  from  it 
to  produce  the  following  equation: 


(3-52b) 


4.  The  sensitivity  coefficient  for  the  coefficient  of  variation  of  the  predicted 
strength  due  to  a  change  in  the  coefficient  of  variation  of  a  basic  random 
variable  (C^qj)  is  computed  as 
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ACOV(Xp)i/COV(Xp) 

ACOV(Xi)/COV(Xi) 


(3-53a) 


c 


OCT, 


n  ( 


ACOV(Xp)i/COV(Xp) 

ACOV(Xi)/COV(Xi) 


f 

/ 


COV(Xp)  =  the  coefficient  of  variation  of  the  predicted  strength;  and 
ACOV(Xp)i  =  change  in  the  coefficient  of  variation  of  the  predicted  strength 
due  to  the  change  in  the  coefficient  of  variation  of  the  ith  basic  random 
variable.  The  coefficients  of  variation  of  the  basic  random  variables  are 
perturbed  one  random  variable  at  a  time,  while  keeping  the  remaining  random 
variables  at  their  non-perturbed  coefficients.  They  are  also  perturbed  using  the 
same  percent  of  their  respective  non-perturbed  values.  The  coefficients  of 
variation  can  be  perturbed  by  perturbing  the  corresponding  standard  deviations 
of  the  basic  random  variables.  The  percent  change  and  COV(Xp)  in  Eq.  3- 
53a  cancel  out  from  it  to  produce  the  following  equation: 


ACOV(Xp)i 

2;(ACOV(Xp)i)^ 

lli=l 


(3-53b) 


3.3.3.  Weighted  Sensitivity  Coefficients 

In  the  previous  section,  the  sensitivity  coefficients  are  defined  as  the  normalized 
changes  in  the  moments  of  the  predicted  strength  due  to  changes  in  the  moments  of  the 
basic  random  variables.  The  first  two  cases  deal  with  the  sensitivity  coefficients  for  the 
predicted  strength  due  to  changes  in  the  mean  values  of  the  basic  random  variables  using  a 
constant  percent  change.  The  latter  two  cases  deal  with  the  sensitivity  coefficients  for  the 
predicted  strength  due  to  changes  in  the  coefficients  of  variation  of  the  basic  random 
variables  using  a  constant  percent  change.  It  needs  to  be  noted  that  in  the  first  two  cases, 
a  constant  percent  change  in  the  mean  values  can  result  in  different  percentiles  (or 
cumulative  probabilities)  for  the  different  random  variables  according  to  their  respective 
probabilistic  characteristics.  In  order  to  impose  a  constant  cumulative  probability,  the 
definition  of  the  first  two  cases  needs  to  be  revised  by  basing  the  changes  in  the  means  on 
a  constant  percent  change  of  their  respective  standard  deviations,  instead  of  a  constant 
percent  change  of  their  respective  means.  The  results  are  the  weighted  sensitivity 
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coefficients  for  the  first  two  cases.  The  latter  cases  do  not  require  any  revision.  The 
revised  first  two  cases  are  as  follows: 

1 .  The  sensitivity  coefficient  for  the  mean  predicted  strength  due  to  a  change  in 
the  mean  value  of  a  basic  random  variable  (Wppj)  is  computed  as 


W, 


AXpj/Xp 

AXj/Xj 


"  '"AXpi/Xp 


AXj  /Xj 


(3-54a) 


where  AXj  -  change  in  the  mean  value  of  the  ith  basic  random  variable;  Xp  = 

the  mean  value  of  the  predicted  strength;  and  AXpj  =  change  in  the  mean  value 
of  the  predicted  strength  due  to  the  change  in  the  mean  value  of  the  ith  basic 
random  variable.  The  means  of  the  basic  random  variables  are  perturbed  one 
random  variable  at  a  time,  while  keeping  the  remaining  random  variables  at 
their  non-perturbed  means.  They  are  perturbed  using  the  same  percent  of  their 

respective  standard  deviations.  The  percent  change  and  Xp  in  Eq.  3-54a 
cancel  out  from  it  to  produce  the  following  equation: 


AXpi 


W, 


COV(Xi ) 


AX 


Pi 


COV(Xi ) 


(3-54b) 


2.  The  sensitivity  coefficient  for  the  coefficient  of  variation  of  predicted  strength 
due  to  a  change  in  the  mean  value  of  a  basic  random  variable  (Wopj)  is 
computed  as 
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ACOV(Xp)i/COV(Xp) 


(3-55a) 


where  COV(Xp)  =  the  coefficient  of  variation  of  the  predicted  strength;  and 
ACOV(Xp)i  =  change  in  the  coefficient  of  variation  of  the  predicted  strength 
due  to  the  change  in  the  mean  value  of  the  ith  basic  random  variable.  The 
means  of  the  basic  random  variables  are  perturbed  one  random  variable  at  a 
time,  while  keeping  the  remaining  random  variables  at  their  non-perturbed 
means.  They  are  perturbed  using  the  same  percent  of  their  respective  standard 

deviations.  The  percent  change  and  COV(Xp)  in  Eq.  3-55a  cancel  out  from  it 
to  produce  the  following  equation; 


ACOV(Xp)i 


(3-55b) 
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4.  FAILURE  MODES  AND  STRENGTH  ASSESSMENT  MODELS 

4.1  Identification  of  Failure  Modes 

As  pointed  out  by  Pittaluga  (1987),  historically  the  principal  obstacles  to  achieving 
reliability-based  design  have  been  (a)  the  lack  of  accurate  and  efficient  algorithms  for 
calculating  the  limit  states,  and  (b)  the  lack  of  computer  implementation  of  these 
algorithms.  But  since  then  both  obstacles  have  been  overcome.  In  1988  SNAME 
published  Ship  Structural  Design  by  Hughes  (1988),  which  combines  and  builds  on  the 
work  of  many  ship  structures  researchers,  and  presents  the  underlying  theory  and  the 
solution  algorithms  for  all  of  the  relevant  failure  modes.  Secondly,  all  of  these  algorithms 
have  been  implemented  in  the  MAESTRO  computer  program  (Hughes  1985),  which 
requires  only  a  PC  and  is  now  widely  used  by  designers,  shipyards,  classification  societies 
and  many  other  organizations  (e.g.,  nine  navies  and  two  coast  guard  agencies). 


Table  4.1  lists  the  failure  modes  of  the  principal  members  of  ship  structures  and  gives 
the  failure  category  and  the  computational  algorithm  source  for  each  of  them. 


Failure  Modes  of  Principal 
Members 

Failure 

Category 

Computational  Algorithm 
Source 

Status  re 
Avail.  Data 

PANEL 

(SSD  denotes  Ship  Structural  Design) 

Collapse 

SSD  Sec.  14.2 

Stiffener  Flexure 

Collapse 

OK 

Combined  Buckling 

Collapse 

SSD  Sec.  13.2  - 13.4 

Insufif. 

Membrane  Yield 

Collapse 

SSD  Sec.  12.5 

N.  A. 

Stiffener  Buckling 

Collapse 

SSD  Sec.  13.1  &  15.5 

Insufif. 

Stiffener  Unserviceability 
(Initial  Yield) 

Tension,  Flange 

Yield 

Beam  Theory  &  SSD  Sec.  8.6 

N.  A. 

Tension,  Plate 

Yield 

tl 

N.  A. 

Compression,  Flange 

Yield 

II 

N.  A. 

Compression,  Plate 

Yield 

II 

N.  A. 

Plate  Unserviceability 

SSD  Sec.  9.1  &  9.2 

Yield,  plate  bending 

Yield 

N.  A. 

Local  buckling 

Unserv. 

SSD  Sec.  12.6 

OK 

Allowable  Permanent  Set 

Yield 

SSD  Sec.  9.3  -  9.5  &  (H&C  91) 

OK 

BEAM 

Collapse 

■■ 

Tripping 

Collapse 

SSD  Sec.  13.1 

Flexural-Torsional  Buckling 

Collapse 

SSD  Sec.  15.4  &  15.5 

Plastic  Hinge 

Collapse 

SSD  Sec.  16.1  &  16.2 

OK 

Unserviceability  (Initial  Yield) 

Bending 

Yield 

Beam  Theory 

N.  A. 

Web  shear 

Yield 

II  II 

N.  A. 

GRILLAGE 

Collapse 

■ 

■■ 

Overall  Buckling 

SSD  Sec.  10.2  &  13.5,6 

Plastic  Hinge 

SSD  Sec.  16.1  - 16.4 

Table  4.1  Identification  of  Failure  Modes  for  Principal  Structural  Members 
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4.2  Selection  of  Failure  Modes  for  this  Project 

4.2.1  Necessity  of  Experimental  Data 

In  general,  the  only  way  to  assess  the  accuracy  of  a  failure  theory  and  its  associated 
computational  algorithm  is  from  physical  data,  either  model  or  full  scale.  That  is,  in 
structural  engineering,  apart  from  a  few  exceptions,  it  is  inadvisable  to  try  to  assess  the 
accuracy  of  a  theory  by  means  of  another  theory.  For  ship  structures  this  important 
principle  raises  a  major  difficulty  because  there  are  many  possible  failure  modes,  and  for 
some  of  them  there  is  not  sufficient  experimental  data.  For  other  types  of  structures  the 
requisite  experiments  have  been  performed.  For  example,  for  box  girder  bridges  the  steel 
portion  (the  box),  has  a  simple  geometry  and  the  number  of  different  possible  failure 
modes  is  small;  therefore  it  has  been  easier  to  obtain  the  necessary  measurements.  Also 
this  testing  process  received  very  large  international  support  in  the  70s  because  of  a  series 
of  structural  failures  of  new  bridges. 

For  aircraft  the  large  production  numbers  make  it  possible  to  conduct  numerous  full 
scale  tests,  including  even  the  testing  to  destruction  of  a  prototype.  In  offshore  structures 
the  enormous  financial  scale  of  the  oil  industry  -  in  both  costs  and  revenues  -  has 
motivated  large  expenditure  by  the  industry  and  by  some  governments  to  pay  for  whatever 
structural  tests  were  needed.  For  example,  over  a  period  of  about  ten  years  the 
American  Petroleum  Institute  (an  industry-wide  funding  agency)  sponsored  the 
development  of  a  completely  new  and  comprehensive  family  of  failure  algorithms  for 
stiffened  cylinders  of  all  relevant  proportions  and  stiffening  arrangements.  This  included  a 
very  comprehensive  series  of  1/4  scale  collapse  tests  -  the  largest  such  testing  program 
ever  conducted  apart  from  submarine  hulls  (for  which  the  results  are  less  relevant  because 
of  the  greater  pressure,  and  are  not  available  anyway). 

4.2.2  Failure  Modes  That  Do  Not  Qualify 

For  some  failure  modes  the  limit  value  is  not  a  calculated  quantity  but  rather  a  specific 
independent  quantity  such  as  yield  stress,  or  a  specified  maximum  allowable  response 
value  under  a  nominated  "design"  load.  Two  examples  are  maximum  (elastic)  deflection 
and  maximum  allowable  permanent  set.  Since  the  limit  value  does  not  contain  any 
calculational  or  modeling  uncertainty,  these  failure  modes  are  not  relevant  for  this 
project.  Yield  stress  does  contain  some  bias  and  uncertainty,  but  it  is  obtained  from 
material  sampling,  which  is  a  separate  aspect  of  rationally-based  design.  In  Table  4. 1 
these  failure  modes  have  "N.A."  (Not  Applicable)  in  the  column  headed  Status  re 
Available  Data".  This  leaves  ten  failure  modes  that  involve  calculated  failure  values.  Of 
these,  four  were  found  to  have  sufficient  data  for  statistical  analysis,  these  are  labeled  OK 
in  the  Table.  The  other  six  do  not  have  sufficient  data,  and  they  are  labeled  "Insuff." 
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4.2.3  Selection  Criteria  and  Final  Choice 


The  final  choice  of  failure  mode  was  based  on  the  degree  of  seriousness,  as  measured 
by  the  following  criteria: 

(1)  size  or  importance  of  the  members 

(2)  sudden  failure  (e.g.,  buckling)  vs.  progressive  failure  in  proportion  to  overload 

(3)  consequence  of  failure:  does  member  collapse  or  merely  become  unserviceable 

(4)  existence  of  alternative  load  paths  (redundancy  of  member) 

Once  these  criteria  are  stated  the  choice  becomes  obvious.  Stiffened  panel  collapse 
due  to  flexural  (beam-column)  failure  of  the  stiffeners,  the  first  failure  mode,  is  more 
serious  than  the  other  three  modes  in  all  of  the  above  criteria.  The  two  modes  of  plate 
unserviceability  are  much  less  serious,  and  the  plastic  hinge  collapse  of  a  beam  is 
proportional  to  the  load.  Also,  in  a  ship  a  beam  is  usually  part  of  a  3D  framework  and  so 
there  are  several  alternative  load  paths. 
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5.  ALGORITHMS  AND  DATA  FOR  COMPRESSIVE  COLLAPSE 

OF  STIFFENED  PANELS 

5.1  Available  Algorithms 

5.1.1  "Standard"  Algorithm 

The  best  (most  thoroughly  validated  and  widely  accepted)  computational  model  for 
this  failure  mode  is  that  developed  in  the  UK  under  the  guidance  and  sponsorship  of  the 
Merrison  Committee,  which  led  to  the  UK  LRFD  Code  for  steel  box  girder  bridges, 
BS5400  (British  Standards  Institute,  1982).  This  computational  model  is  directly 
applicable  to  ship  panels  and  is  presented  for  that  application  in  Section  14.2  of  Ship 
Structural  Design  (SSD).  It  has  also  been  adopted  by  Lloyds  Register  of  Shipping  and  is 
implemented  in  the  LR  computer  program  for  panel  strength;  LR  Pass  Program  20202.  It 
is  also  implmented  in  the  MAESTRO  program.  Because  of  its  thorough  validation  and 
wide  acceptance,  we  will  herein  refer  to  this  as  the  "Standard"  algorithm. 

The  Standard  algorithm  is  presented  and  validated  in  Chapter  14  of  SSD,  and  this 
chapter  is  included  herein  as  Appendix  B.  In  this  model  each  stiffener  is  regarded  as  an 
isolated  beam-column,  with  the  plating  acting  as  one  of  the  two  flanges.  If  the  stiffener  is 
a  tee  section,  then  the  beam-column  is  monosymmetric.  Because  of  the  relatively  large 
width  of  the  plate  flange,  the  neutral  axis  of  the  beam-column  is  close  to  the  plating;  and 
hence  the  largest  flexural  stress  occurs  in  the  stiffener  flange.  Because  of  the  unsymmetry 
about  the  neutral  axis,  the  two  directions  of  primary  bending  (bending  in  the  plane  of  the 
web)  have  quite  different  consequences.  Figure  14.2  of  Appendix  B  illustrates  the  two 
cases.  When  the  bending  deflection  is  toward  the  plating,  the  flexural  stress  in  the  stiffener 
flange  is  compressive,  and  it  combines  with  the  applied  axial  compressive  stress,  so  that 
the  stiffener  flange  is  the  most  highly  stressed  location  in  the  cross  section.  Eventually  it 
reaches  the  yield  stress  and  the  beam-column  collapses.  In  this  case  we  speak  of  the 
collapse  as  being  "stiffener-induced".  In  the  terminology  of  SSD,  this  is  called  a  Mode  I 
collapse.  In  contrast,  when  the  bending  deflection  is  toward  the  stiffener,  the  compressive 
flexural  stress  now  occurs  in  the  plating  and  it  combines  with  the  applied  axial 
compressive  stress.  Eventually  this  combined  compressive  stress  causes  the  plate  to 
buckle  and  the  beam-column,  having  lost  its  major  flange,  simultaneously  collapses.  In 
this  case  we  speak  of  the  collapse  as  being  "plate-induced".  In  the  terminology  of  SSD, 
this  is  called  a  "Mode  11"  collapse. 

The  accuracy  of  the  Standard  algorithm  over  all  combinations  of  lateral  and  in-plane 
loads  is  demonstrated  in  Section  14.4  of  Appendix  B.  Figure  14.1 1  gives  a  good  overall 
summary  of  the  results. 
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5.1.2  Other  Algorithms 

Prior  to  the  publication  of  the  Standard  algorithm,  various  researchers  had  developed 
five  other  methods:  Faulkner  (1975B),  Murray  (1975),  Carlsen  (1980),  Dwight  and  Little 
(1976)  and  Horne  and  Narayanan  (1977)..  Faulkner's  method  (1975B)  is  based  on  a 
Johnson-Ostenfeld  approximation  together  with  an  effective  width  approach  for  plate 
behavior.  This  method  was  also  presented  in  Faulkner  et  al  (1973).  The  methods  of 
Murray  (1975),  Carlsen  (1980),  Dwight  and  Little  (1976),  and  Home  and  Narayanan 
(1977)  are  all  based  on  a  Perry-Robertson  formulation,  also  with  an  effective  width  model 
for  the  plating.  The  main  differences  in  these  methods  are  the  ways  they  account  for  the 
effects  of  residual  stress,  initial  imperfections,  and  eccentricity  due  to  loss  of  plate 
effectiveness. 

5.1.3  Smith  Algorithm 

A  quite  different  algorithm  is  presented  in  (Smith  1975),  which  makes  use  of  an 
extremely  detailed  model  and  large  amounts  of  iteration.  For  example,  the  web  of  each 
stiffener  is  divided  vertically  into  ten  or  more  "layers"  or  "zones",  and  a  new  solution  is 
performed  each  time  the  stress  distribution  changes  in  any  zone.  The  solution  also 
involves  incremental  stepping  along  a  series  of  idealized  stress-strain  curves  for  the 
plating.  The  total  computation  is  much  larger  than  in  the  Standard  algorithm,  and  yet  as 
shown  on  pages  484,5  of  SSD,  the  Standard  algorithm  gives  comparable  results,  provided 
that  it  is  used  in  close  association  with  a  three-dimensional  finite  element  analysis  of  the 
stmcture,  so  that  the  tme  boundary  conditions  for  the  panel  are  known.  This  is  the  way 
the  Standard  algorithm  is  used  in  the  MAESTRO  program,  and  the  finite  element  analysis 
is  not  an  extra  computational  burden  because  it  must  be  performed  in  any  case  in  order  to 
determine  what  are  the  actual  or  working  stresses  in  the  panel,  without  which  the  panel 
could  not  be  designed  or  even  assessed  as  to  its  adequacy. 

In  view  of  the  large  computational  burden  of  the  Smith  algorithm  and  the  fact  that  the 
computer  program  for  it  is  not  publicly  or  commercially  available,  the  Smith  algorithm  is 
not  suitable  for  the  purposes  of  this  project. 

5.2  Summary  of  Current  Code-based  Design  Practices 

Since  the  ultimate  goal  of  the  SSC  research  program  is  to  develop  a  probability-based 
design  procedure  for  ship  structures,  we  thought  it  would  be  helpful  to  briefly  summarize 
the  current  situation  in  regard  to  code-based  (or  rule-based)  design,  because  some 
organizations  have  already  adopted  probabilistic-based  design  guidelines,  and  others  are  in 
a  state  of  transition. 

5.2.1  AISC  Load  and  Resistance  Factor  Design  Code 

The  AISC  Load  and  Resistance  Factor  Design  (LRFD)  Specification  was  first 
introduced  in  1986.  The  specification  was  developed  under  the  leadership  of  Galambos 
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(e.g.,  Pinkham  and  Hansell  1978,  Galambos  1972,  Galambos  and  Ravindra  1978, 
Galambos  1981,  and  Ravindra  and  Galambos  1978).  The  general  format  for  the  LRFD 
specification  is 


(tiRiXriLi  (5-1) 

/=1 

where  (j)  =  strength  reduction  factor,  R  =  resistance  or  strength  according  to  some  failure 
mode,  g  =  load  amplification  factor,  L  =  load  type,  and  n  =  number  of  combined  loads. 

The  development  of  the  LRFD  code  was  based  on  a  probability-based  model  (Galambos 
and  Ravindra  1978,  Ravindra  and  Galambos  1978,  Ellingwood,  et  al.  1980,  and 
Galambos,  et  al.  1982),  calibration  with  the  1978  ASIS  Allowable  Stress  Design 
Specification,  and  expert  judgment  based  on  previous  design  experiences.  In  developing 
the  specification,  it  was  necessary  to  change  the  design  practice  from  working  stress  to 
limit  stress,  and  from  allowable  to  ultimate  which  was  also  reliability-based. 

In  preparation  for  the  development  of  the  AISC  LRFD  specification,  uncertainties  in 
strength  parameters  of  steel  structures  were  investigated  and  summarized  by  Ellingwood 
(1980).  Also,  the  study  included  other  materials  such  as  timber,  and  reinforced  and 
prestressed  concrete.  The  results  of  this  study  were  used  to  quantify  strength  uncertainties 
in  marine  structures. 

5.2.2  AASHTO  Code 

Currently,  the  AASHTO  Specification  is  being  revised  to  an  LRFD  format.  The 
National  Cooperative  Highway  Research  Program  (NCHRP)  has  published  the  third  Draft 
of  LRFD  Specifications  and  Commentary  in  1992  titled  "Development  of  Comprehensive 
Bridge  Specifications  and  Commentary."  The  document  is  marked  not  for  distribution 
pending  review  and  approval.  It  is  widely  expected  that  the  AASHTO  LRFD  code  will 
closely  follow  much  of  the  AISC  code.  Many  of  the  individuals  instrumental  in  the 
development  of  the  AISC  LRFD  code  are  involved  with  the  AASHTO  effort. 


5.2.3  API  RP  2A  Recommended  Practice 


The  eighteenth  edition  (September,  1989)  of  the  API  Recommended  Practice  for 
Planning,  Designing  and  Constructing  Fixed  Offshore  Platforms  is  not  yet  an  LRFD  code. 
The  code  is  primarily  concerned  with  the  design  of  tubular  members  and  joints.  The  code 
makes  extensive  reference  both  to  the  AISC  Code  and  to  recommended  practices  of  the 
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American  Welding  Society  (AWS),  though  it  specifically  warns  against  using  the  new 
AISC  LRFD  Code.  This  is  because  the  load  and  resistance  factors  used  in  the  AISC  Code 
are  based  on  a  calibration  with  building  design  practices.  There  is  research  currently  in 
progress  to  attempt  to  develop  load  and  resistance  factors  which  would  be  useful  in  an 
offshore  design.  In  the  area  of  fatigue  design  the  code  relies  on  extensive  studies 
performed  by  Wirsching  and  others  to  take  into  account  the  probabilistic  nature  of  the 
random  loading.  However,  the  RP  2A  fatigue  code  still  recommends  using  a  design  life  of 
two  times  the  expected  life  and  requires  that  the  cumulative  damage  ratio  remains  below  1 
for  the  design  life  of  the  structure. 

/ 

5.2.4  U.S.  Navy  Design  Data  Sheets 

In  naval  ship  structural  design  the  current  design  methods  do  not  address  failure 
modes  or  strength  values  explicitly.  Instead  the  U.  S.  Navy  uses  a  series  of  "design  data 
sheets"  which  spell  out  a  particular  sequence  of  calculations,  the  result  of  which  is  either 
the  (minimum)  required  member  sizes  or  minimum  values  of  size-related  member 
properties,  such  as  section  modulus,  from  which  the  sizes  may  be  determined.  There  are 
usually  several  such  calculation  procedures  to  be  performed,  and  each  procedure  has  been 
developed  starting  with  the  assumption  (based  on  experience)  that  for  that  particular  type 
of  ship,  one  or  two  particular  failure  modes  will  govern  the  design. 

In  the  U.S.  Navy  the  two  principal  documents  for  the  design  of  structural  members  are 
DDS  100-4  Strength  of  Structural  Members,  z.ndTyDS-\00-6  Longitudinal  Strength 
Calculations.  Both  of  these  documents  have  been  in  place  for  a  long  time,  with  the  newer 
DDS  100-4  having  an  effective  date  of  February  1979.  The  principal  references  for  DDS 
100-4  date  back  to  the  work  of  Dr.  Friedrich  Bleich  in  the  1950's  (Bleich,  1952;  Bleich 
and  Ramsey,  1958).  While  the  procedure  in  DDS  100-4  has  its  roots  in  sound  analysis, 
simplistic  modeling  of  the  interaction  of  load  components  and  arbitrary  factors  of  safety 
render  the  actual  level  of  safety  present  in  the  final  design  results  an  unknown  quantity. 
There  is  an  effort  underway  through  the  Naval  Sea  Systems  Command  to  update  and 
improve  the  Design  Data  Sheets.  It  is  very  likely  that  the  update  will  be  in  the  form  of  a 
reliability-based  design  code. 

5.2.5  American  Bureau  of  Shipping  Rules 

The  American  Bureau  of  Shipping  published  annually  the  latest  revision  to  its 
Guidelines  for  the  construction  of  steel  merchant  vessels.  Though  many  of  the  procedures 
and  formulae  in  the  publication  are  empirically-based,  there  has  been  a  steadily  increasing 
effort  to  introduce  a  probabilistic  basis  into  the  selection  of  some  of  the  coefficients  used 
in  the  design  equations.  To  help  to  make  the  rules  more  accessible  and  readable,  ABS  has 
begun  to  publish  the  rules  in  a  computer  CD  ROM  format.  With  the  search  and  retrieval 
software  available  with  the  rules,  it  is  significantly  easier  to  find  the  needed  topic  and  all 
relevant  factors  to  the  design  of  a  particular  component. 
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Their  most  recent  effort  to  update  and  improve  the  classification  procedure  has  been 
the  introduction  of  the  ABS  SafeHull  System  (Grove  and  Coulter,  1993).  Though 
available  only  for  tankers  at  the  time  of  this  report,  there  is  an  effort  underway  to  produce 
a  similar  system  for  bulk  carriers.  The  system  is  a  suite  of  computer  programs  aimed  at 
assisting  the  structural  designer  in  evaluating  design  alternatives  and  producing  a  safe 
design  suitable  for  classification.  The  strength  assessment  in  the  programs  covers  the 
failure  modes  of  yielding,  buckling,  and  fatigue.  The  SafeHull  System  is  a  computerized 
implementation  of  the  strength  assessment  procedures  published  in  the  “Guide  for 
Dynamic  Based  Design  and  Evaluation  of  Tanker  Structures”  and  the  “Guide  for  the 
Fatigue  Strength  Assessment  of  Tankers”  (ABS,  1993a  and  1993b). 

5.2.6  Lloyd’s  Register  of  Shipping  Rules 

The  complete  set  of  LR  Rules  for  ship  classification  is  now  available  on  CD  ROM. 

For  this  project  the  1989  edition  of  the  Lloyd’s  Rules  was  investigated.  The  ship 
structures  related  issues  are  contained  in  Parts  3  and  4.  The  general  organization  and  type 
of  rules  are  similar  to  those  of  ABS.  There  is  no  explicit  assessment  of  strength  in  specific 
failure  modes,  rather  the  procedure  leads  one  to  design  a  safe  structure  by  simple  formulae 
aimed  at  defining  specific  scantlings.  As  a  result,  information  about  levels  of  safety 
included  in  the  design  are  not  apparent.  However,  Lloyd’s  does  allow  for  direct 
calculations  either  through  its  LR  PASS  program  or  through  the  designer’s  own 
calculation  procedure.  In  the  latter  case  the  designer  would  be  required  to  provide  details 
of  the  calculations  and  any  programs  used  for  the  analysis. 

5.2.7  Det  Norske  Veritas  Rules 

Det  Norske  Veritas’s  role  is  similar  to  that  of  ABS  and  Lloyd’s.  Though  not  widely 
used  in  the  United  States  for  ship  design,  DnV  Rules  are  extensively  used  in  the  offshore 
field.  Most  of  the  calculations  in  the  DnV  rules  are  similar  in  nature  to  the  other 
classification  societies  -  that  is  they  are  aimed  at  assisting  the  designer  in  developing 
adequate  scantlings  without  burdening  him  with  the  analysis  of  which  failure  mode  or 
modes  is  the  most  critical.  The  basis  for  much  of  the  calculations  and  equations  in  the 
DnV  rules  can  be  found  in  the  Ships  ’  Load  and  Strength  Manual  (Lersbryggen,  1978). 
This  manual  contains  a  comprehensive  set  of  approximate  formulae  for  the  strength  of 
stiffened  panels.  These  formulae  are,  in  general,  developed  from  the  underlying  theory 
through  the  use  of  suitable  engineering  approximations.  While  very  useful,  the  manual 
does  not  provide  any  new  theory  or  approach  to  the  strength  assessment  of  stiffened 
panels. 

5.3  Available  Test  Data 

The  main  focus  of  this  study  is  the  uncertainty  in  the  strength  of  longitudinally 
stiffened  panels  under  axial  compression.  The  data  sets  used  in  this  study  are  the  those 
which  were  readily  obtainable  and  which  had  sufficient  data  for  the  purpose.  Other 
sources  of  data  were  either  found  to  be  proprietary  or  of  insufficient  scope  or  lacking 
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some  of  the  required  data.  There  may  well  be  other  sources  that  we  did  not  find  due  to 
limited  publication.  In  the  end,  sufficient  data  was  found  to  be  available  only  for  single 
bay  panels  under  purely  axial  compression.  Even  here  there  were  only  three  test  series 
which  fulfilled  the  requirements  -  the  nominally  identical  and  parametric  series  reported  in 
Faulkner  (1977)  and  the  panels  reported  in  Michelutti  (1977). 

The  tests  reported  in  Smith  (1975)  are  for  three-bay  panels,  and  in  general  they 
were  each  different  and  they  each  involved  a  different  mode  of  failure,  including  stiffener 
tripping  and  grillage  buckling.  The  only  exception  was  one  pair  of  nominally  identical 
panels.  Unfortunately  there  are  very  few  other  three-bay  panel  tests,  and  none  which  are 
sufficiently  similar  in  loading  and  geomery  to  be  grouped  with  Smith's  tests. 

Faulkner  (1977)  reported  on  four  test  programs:  the  student  series,  the  nominally 
identical  series,  the  residual  welding  stress  series,  and  the  parametric  series.  The 
nominally  identical  and  parametric  series  data  are  analyzed  in  this  study.  The  student 
series  was  deemed  inappropriate  for  use  as  a  result  of  vagueness  in  the  reporting.  The 
residual  stress  series  was  conducted  to  investigate  the  effect  of  residual  stress  on  the 
overall  strength  of  a  welded  panel,  and  consists  of  two  sets  of  panels.  This  series  was  not 
used  herein  because  the  first  set  was  stress  relieved  by  baking,  and  the  second  set  is  not 
fully  documented.  The  second  set  was  constructed  similarly  to  the  nominally  identical 
series,  but  the  welding  methods  were  varied  to  create  different  amounts  of  residual  stress. 
The  dimensions  and  material  properties  are  not  reported  for  the  panels,  nor  are  their 
statistics  aside  from  the  mean  values. 

The  nominally  identical  tests  were  conducted  to  quantify  the  strength  uncertainties 
resulting  from  normal  shipyard  construction  methods.  The  panels  were  all  constructed 
from  the  same  design  with  variation  resulting  purely  from  fabrication  and  material 
randomness.  The  panels  consisted  of  plating  reinforced  with  five  stiffeners,  oriented 
longitudinally  approximating  quarter  scale  deck  panels  of  one  bay  length.  Each  panel  was 
loaded  axially  to  compressive  failure.  In  order  to  make  the  panels  are  more  sensitive  to 
variations  in  dimensions  and  residual  stress,  they  were  designed  such  that  the  stiffener 
slenderness  parameter,  was  of  the  order  of  0.9,  and  the  plate  slenderness  parameter, 

P,  was  of  the  order  of  2. 1 .  Faulkner(1977)  presented  an  analysis  of  the  subjective  and 
objective  uncertainties  associated  with  these  tests.  In  Section  6.3.2  we  compare  our 
analysis  with  his  and  give  reasons  why  we  believe  ours  to  be  more  accurate. 

The  parametric  series  tests  were  conducted  to  examine  the  effect  on  ultimate 
strength  of  varying  slenderness  ratios.  The  slenderness  ranges  were  1, 2,3,4  for  P  and 
0.375,  0.750,  1.125  for  Xg.  They  were  of  the  same  scale  as  the  nominally  identical  series, 
constructed  from  plating  with  5  longitudinal  stiffeners,  and  of  one  bay  length.  The 
loading  was  purely  axial  and  increased  until  ultimate  failure.  The  construction  methods 
were  consistent  for  each  panel.  All  but  three  of  the  panels  had  tee  stiffeners. 
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The  tests  conducted  by  Michelutti  (1977)  consisted  of  single  bay  models  with 
combined  axial  and  lateral  pressure  loads.  Of  the  15  models  reported,  three  were  axially 
loaded  and  two  of  these  were  used  in  this  study. 

The  details  of  all  of  the  panels  and  test  programs  are  discussed  in  the  following 
chapter. 
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6.  DEMONSTRATION  OF  UNCERTAINTY  ASSESSMENT  FOR 
COLLAPSE  OF  STIFFENED  PANELS 


In  this  chapter  we  demonstrate  the  estimation  of  random  and  modeling  uncertainties  in 
the  strength  of  stiffened  panels  that  fail  under  inelastic  collapse,  using  data  from  two  series 
of  tests  reported  by  Faulkner  (1977)  and  one  series  reported  by  Michelutti  and  Murray 
(1977).  Section  6.1  describes  the  tests  and  the  panels  used  in  these  tests. 

6.1  Description  of  Panels 

6.1.1  Faulkner’s  nominally  identical  series 

The  nominally  identical  series  of  panel  tests  in  Faulkner  (1977)  were  conducted  to 
examine  strength  uncertainties  in  longitudinally  stiffened  panels  under  axial  compressive 
loading.  The  models  were  fabricated  at  a  naval  shipyard  under  normal  conditions  and 
tested  to  failure.  Of  the  13  models  tested,  4  failed  in  stiffener-induced  collapse  (Mode  I) 
and  9  failed  in  plate-induced  collapse  (Mode  II).  (For  brevity  we  will  use  the  "Mode" 
terminology;  a  description  of  the  two  modes  was  given  in  Section  5.1.1).  The  four  panels 
that  failed  in  Mode  I  initially  had  unusually  large  eccentricities,  and  they  were  straightened 
out  after  construction,  and  before  testing,  by  means  of  an  unspecified  amount  of  heating. 
We  did  not  include  these  four  panels  because  of  the  artificial  and  completely  unknown 
amount  of  bias  that  they  contain. 

Figure  6-1  depicts  a  panel  in  the  nominally  identical  series.  The  panels  were  designed 
to  represent  a  quarter  scale  deck  panel  and  were  stiffened  by  five  stiffeners  cut  from  3"xl" 
Ministry  of  Defense  (UK)  Standard  Tees.  Simply  supported  boundary  conditions  were 
simulated  by  applying  the  load  through  knife  edges  at  the  elastic  neutral  axis.  End  caps 
were  used  to  shed  the  load  from  the  knife  edge  into  the  cross  section.  The  stiffness  added 
to  the  panels  by  these  end  caps  shorten  the  simply  supported  length  for  analysis  purposes 
to  the  reported  effective  length.  The  effective  length  formulation  is  treated  in  Faulkner 
(1977). 

As  the  panels  were  nominally  identical,  any  differences  in  the  panels  were  due  to 
randomness  in  the  variables.  The  statistics  of  the  measured  material  properties  and  the 
dimensions  of  the  panels  are  reported  in  Table  6-1,  along  with  the  nominal  (design)  values. 
The  actual  scantlings  of  each  panel  were  not  reported  in  Faulkner  (1977).  All  random 
variables  in  Table  6-1  were  assumed  to  be  normal.  The  effective  length  (ag),  stiffener 
spacing  (b),  and  stiffener  flange  breadth  (bf)  were  all  considered  constant,  as  follows. 


ag  =  790.0  mm, 
b  =  183.0  mm, 
bf  =  25.4  mm. 
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We  used  effective  length  in  all  calculations.  For  simplicity,  in  the  rest  of  the  section 
we  will  refer  to  effective  length  using  the  term  "length"  and  denote  it  by  "a". 

The  plate  free  edges  had  a  width  of  b/3  to  approximate  the  same  buckling  stress  as  the 
supported  plate. 


Random  Variable 

Nominal 

Value^ 

Mean 

Standard  Deviation 

COV  (%) 

3.0  mm 

3.07  mm 

0.027937  mm 

0.91 

Depth  of  stiffener 

(d) 

32.5  mm 

32.5  mm 

0.26975  mm 

0.83 

4.4  mm 

4.90  mm 

0.01519  mm 

0.31 

6.4  mm 

5.84  mm 

0.012848  mm 

0.22 

liMliil 

8.3 

5.6 

Eccentricity  (e) 

±0.52732 

±0.52732 

100. 

Tension  zone 
parameter,  r\,  due  to 
residual  stresses^ 

2.772 

2.7095 

0.294 

10.85 

Young's  modulus 

1.5 

Notes: 

^Nominal  dimensions  from  Ministry  of  Defense  (UK)  Standard  Tee  Bar  catalog 

^The  mean  and  nominal  values  were  estimated  to  be  -a/ 1500  for  Mode  I  and  +a/1500 
for  Mode  II. 

^The  standard  deviation  of  the  eccentricity  was  assumed  to  be  equal  to  the  length  of 
the  panel  over  1500. 

4pages  393-394  of  Hughes  (1988)  address  this  parameter.  The  nominal  and  mean 
values  of  rj  are  based  on  the  assumption  that  Ci-Zay  =  0. 10,  where  Cj-  is  the  compressive 
residual  stress  and  Oy  is  the  yield  stress  of  the  plating. 

Table  6-1 .  Uncertainties  in  Geometry  and  Material  Properties  of  Stiffened  Panel 

The  compressive  residual  stress,  which  was  measured  after  fabrication,  corresponds  to 
a  tension  zone  parameter  (ri)  of  4.2.  The  residual  stress  is  not  directly  specified  in  the 
algorithm  that  we  used  to  predict  strength,  where  the  residual  stress  is  fixed  at  a  value 
approximately  equal  to  ten  percent  of  yield  stress  of  the  plating.  In  this  study  we  used  a 
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mean  tension  zone  parameter,  r|  =  2.710,  which  corresponds  to  a  residual  stress  equal  to 
ten  percent  of  the  yield  stress. 

The  statistics  of  the  panel  eccentricity  were  not  given  and  so  we  had  to  use  an 
assumed  distribution,  based  on  the  following  statement  from  Section  2.2  of  Faulkner 
(1977):  "the  stiffener  deformations  were  very  small,  usually  much  less  than  a/1000  and  in 
no  case  greater  than  a/500...".  We  therefore  assumed  mean  eccentricities  of -a/1500  for 
Mode  I  and  +a/1500  for  Mode  II,  with  a  standard  deviation  of  a/1500.  To  check  if  our 
assumption  on  the  standard  deviation  of  the  panel  eccentricity  was  reasonable,  we 
estimated  the  standard  deviation  of  the  eccentricity  of  fourteen  panels  studied  by 
Mchelutti  and  Murray  (1977).  The  standard  deviation  was  found  to  be  a/1376,  which 
indicated  that  our  assumption  was  reasonable. 

Finally,  it  is  observed  from  Table  6-1  that  the  eccentricity,  the  yield  stresses  of  the 
plate  and  stiffener  and  the  residual  stress  have  the  largest  coefficients  of  variation. 


6.1.2.  Faulkner's  parametric  series 

The  second  set  of  panels  used  in  this  study  is  the  set  referred  to  as  the  parametric 
series  in  Faulkner  (1977).  In  this  series  24  panels  were  tested  to  study  the  effect  of 
varying  slenderness  ratios  on  the  column  strength.  19  of  these  panels  failed  in  Mode  II 
and  were  used  in  this  study.  Of  the  five  not  used,  2  experienced  tripping  failure  and  3  had 
flat  bar  stiffeners  instead  of  T-shaped  stiffeners. 

The  panel  dimensions  and  properties  were  reported  for  each  panel  and  are  given  in 
Table  6-2.  The  panels  had  the  same  layout  as  the  nominally  identical  series,  with  5 
stiffeners  aligned  longitudinally  and  b/3  free-edge  plate  width.  The  stiffeners  were  also  cut 
from  3"  X  1"  Admiralty  tees.  Residual  compressive  stresses  were  measured  and  found  to 
correspond  to  a  range  of  r|  from  3  to  4.5.  For  the  strength  predictions  r|  was  defaulted  to 
2.710  as  in  the  nominally  identical  series.  The  test  procedure  was  identical  with  the 
exception  of  the  load  being  applied  at  the  midpoint  of  the  distance  between  the  initial 
(elasffc)  neutral  axis  and  the  predicted  final  neutral  axis.  The  predicted  final  neutral  axis 
location  was  based  on  an  assumption  of  40t  of  the  plating  being  effective. 

The  panel  mid  span  eccentricity  was  reported  for  each  panel.  For  the  nominal 
calculations  the  mean  eccentricity  was  set  to  one  standard  deviation  (a/1500),  and  it  was 
negative  for  Mode  I  and  positive  for  Mode  II. 
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Note:  Yield  stresses  of  stiffener  and  plate  are  expressed  in  N/mm^ 


Table  6-2:  Properties  of  Panels  in  Parametric  Series 

6.1.3  Panels  A6  and  H:  Michelutti  (1977) 

Fifteen  panels  were  tested  in  a  variety  of  loading  conditions.  Three  of  these  panels 
were  loaded  in  purely  axial  compression.  Of  these  three,  a  panel  called  "panel  J"  failed  m 
Mode  I  and  the  remaining  two  panels  called  "Panels  A6  and  H"  failed  in  mode  II.  Panel  J 
was  predisposed  to  fail  in  Mode  I  by  moving  the  point  of  load  application  closer  to  the 
free  edge  of  the  stiffeners  from  the  centroid.  The  two  other  panels,  A6  and  H,  were  not 
arranged  to  favor  either  failure  mode. 

Panels  A6  and  H  share  similar  dimensions  and  properties  with  the  most  significant 
difference  being  overall  length.  The  panels  are  stiffened  by  five  bulb  stiffeners  and  are 
loaded  at  the  elastic  neutral  axis.  The  panel  properties  are  listed  in  Table  6-3. 
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inel 

m 

m 

— 

n 

WBM 

101 

m 

HBI 

■B&l 

RuSal 

HSSSn 

■Ba 

Kuun 

■ESI 

\6 

4930 

530 

9.60 

152.0 

7.31 

28.6 

15.9 

368 

383 

H 

3450 

530 

9.47 

152.0 

7.08 

28.6 

15.9 

378 

389 

Note:  Yield  stresses  are  expressed  in  N/mm^ 

Table  6-3.  Properties  of  Panels  from  Michelutti  (1977) 

6.2  Results  for  random  uncertainty 

We  estimated  the  uncertainty  in  the  compressive  strength  of  the  stiffened  panels 
described  in  Section  6. 1  due  to  random  uncertainties  in  geometry  and  material  properties 
for  two  failure  modes:  stiffener  initiated  inelastic  buckling  (Mode  I)  and  plate  initiated 
inelastic  buckling  (Mode  II).  Section  6.2.1  presents  the  results.  We  also  estimated  the 
effect  of  each  random  uncertainty  in  geometry  and  material  properties  on  the  strength,  and 
ranked  these  uncertainties  in  terms  of  importance.  Those  results  are  presented  in  Section 
6.2.2. 

6.2.1.  Effect  of  random  uncertainties  on  strength 

Several  combinations  of  axial  and  transverse  in-plane  loads  and  lateral  pressure  were 
studied  to  cover  the  range  of  realistic  combinations  of  these  loads.  Specifically,  the  lateral 
pressure,  p,  ranged  from  -0.07  MPa  to  0.07  MPa  (positive  pressure  is  from  the  plate  to  the 
stiffener).  We  considered  four  cases  where  transverse  stress  was  0,1/16,  1/8  and  1/4  of 
the  applied  axial  stress. 

The  statistics  of  the  strength  were  estimated  using  Monte-Carlo  simulation.  In  each 
replication,  we  predicted  the  axial  strength  of  a  panel  corresponding  to  a  set  of  samples  of 
the  random  variables  using  the  Standard  algorithm  described  in  Section  5.1.1  and 
presented  in  full  in  Appendix  B. 

The  Standard  algorithm  assumes  a  fixed  design  value  for  residual  stress,  equal  to  ten 
percent  of  the  yield  stress.  In  contrast,  Faulkner's  algorithm  (1973)  explicitly  allows  for 
residual  stress.  Therefore,  in  order  to  account  for  the  uncertainty  in  residual  stress,  we 
used  two  separate  Monte-Carlo  simulations.  In  the  first  simulation  we  used  Faulkner's 
algorithm  (1973)  to  predict  strength,  and  we  estimated  the  COV  of  the  strength  due  to 
residual  stress  only,  COVresiduah  variables  at  their  mean  values.  In  the 

second  simulation,  we  used  the  Standard  algorithm  to  estimate  the  COV  of  the  strength 
due  to  all  of  the  other  factors,  COVother-  Then  we  obtained  the  total  uncertainty  in  the 
strength,  COVj,  using  the  following  equation,  which  is  derived  from  eq.(3.7): 


COVT  =  (COvJsidri  +COV^ther)‘'^  (6-1) 
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In  this  equation,  COVj  is  the  COV  of  the  strength  due  to  all  uncertainties.  For  Mode 
I,  the  residual  stress  does  not  significantly  affect  the  strength.  Therefore,  we  assumed  that 
COVx  =  COVother  for  Mode  I. 

A  panel  fails  under  Mode  I  or  Mode  II  depending  on  the  direction  of  the  eccentricity 
and  lateral  pressure.  As  a  result,  some  of  the  samples  generated  in  one  simulation  failed 
under  Mode  I  and  the  others  under  Mode  II.  We  separated  the  samples  depending  on  the 
mode  under  which  they  failed  and  estimated  the  mean  and  COV  of  each  mode  separately. 

The  mean  value  of  the  strength  for  Mode  I  was  found  to  be  247.5  MPa  and  the  COV 
10.20%.  The  corresponding  values  for  Mode  II  were  156.5  MPa  and  6.33%.  The  mean 
strength,  the  nominal  strength  and  the  mean  strength  plus  and  minus  one  and  two  standard 
deviations  are  plotted  in  Figures  6-2  to  6-9  as  a  function  of  the  lateral  pressure.  The 
nominal  strength  was  defined  as  the  smallest  of  the  strengths  of  Modes  I  and  II  that 
correspond  to  the  nominal  values  of  the  variables  in  Table  6-1 .  Figures  6-2  to  6-5  are  for 
applied  transverse  stress  equal  to  0  (case  A),  1/16  (case  B),  1/8  (case  C)  and  1/4  (case  D) 
of  the  axial  stress,  respectively.  The  corresponding  results  for  Mode  II  are  shown  in 
Figures  6-6  to  6-9. 

The  curves  for  the  mean  strength  and  the  mean  plus  and  minus  one  and  two  standard 
deviations  extend  to  the  range  of  positive  pressures  from  0  to  0.035  MPa  for  Mode  I 
(Figure  6-2)  because,  although  positive  pressures  favor  Mode  II,  the  panel  may  still  fail 
under  Mode  I  in  this  range  because  the  eccentricity  is  random.  Let  us  consider  a  panel 
that  has  failed.  For  each  value  of  the  lateral  pressure  in  the  range  from  0  to  0.035  MPa 
each  collapse  mode  has  a  conditional  probability  of  occurrence  (conditioned  on  the  event 
that  the  panel  has  failed),  which  depends  on  the  pressure  and  the  probability  distribution  of 
the  eccentricity.  The  summation  of  the  conditional  probabilities  of  occurrence  of  the  two 
modes  is  one  for  a  given  value  of  the  pressure.  For  positive  pressures  the  conditional 
probability  of  occurrence  of  mode  I  decreases  as  the  pressure  increases  because  positive 
pressure  favors  plate  induced  collapse.  The  interaction  curves  for  Mode  II  extend  to  the 
range  of  negative  pressures  for  the  same  reason. 


We  observe  the  following  from  the  Figures. 

1 .  The  standard  deviation  of  the  strength  for  Mode  I  changes  considerably  with  pressure. 
For  example,  when  there  is  no  applied  transverse  stress,  the  standard  deviation 
changes  from  about  1 1  MPa  to  27  MPa  when  the  pressure  changes  from  -0.07  MPa  to 
+0.35  MPa  (Figure  6-2). 

2.  The  range  of  variation  of  the  standard  deviation  of  the  strength  for  Mode  I  is  smaller 
for  larger  values  of  the  applied  transverse  stress  than  for  smaller  values  (Figures  6-2  to 
6-5). 
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3 .  The  variability  in  the  strength  due  to  random  uncertainty  is  almost  independent  of  the 
applied  lateral  pressure  for  Modes  II.  (Figures  6-6  to  6-9). 

4.  The  nominal  strength  is  considerably  smaller  than  the  mean  strength  for  pressures  from 
-0.07  to  -0.0175  MPa,  which  means  that  if  we  predict  the  strength  using  the  nominal 
values  of  the  random  variables  we  will  be  on  the  conservative  side  (see  Figures  6-2  to 
6-5).  We  believe  that  the  main  reason  is  that  the  nominal  yield  stress  of  the  stiffener  is 
significantly  lower  than  the  measured  value  (Table  6-1). 

5.  The  nominal  strength  is  almost  identical  to  the  mean  for  pressures  from  -0.0175  to 
0.07  MPa  (Figures  6-6  to  6-9).  This  is  because  some  of  the  nominal  values  used  to 
predict  the  strength  were  conservative  (e.g.,  plate  and  web  thickness),  and  others  were 
non  conservative  (e.g.,  plate  yield  stress)(Table  6-1),  and  yet  close  to  the  measured 
values. 


6.2.2.  Ranking  of  random  uncertainties 

We  ranked  uncertainties  in  terms  of  importance  using  the  three  methods  described  in 
Section  3.3:  Parametric  Analysis  and  the  two  methods  using  Sensitivity  Coefficients 
(unweighted  and  weighted).  The  coefficients  measure  the  effect  of  random  uncertainties 
in  strength  parameters  on  strength.  They  can  be  used  to  allocate  resources  in  design  and 
in  collecting  data  for  the  statistical  properties  of  the  random  variables.  In  the  following  we 
present  the  results; 

Parametric  Analysis 

This  method  finds  the  effect  of  each  uncertainty  on  strength  when  the  corresponding 
variable  is  random  and  the  remaining  variables  are  fixed  at  their  means.  Uncertainties  are 
ranked  using  the  corresponding  parametric  coefficients,  Pij,  that  are  calculated  from 
equation  3-48. 


Pli  = 


COV(Xp)i 

Fn 

2(COV(Xp)i)^ 

i=l 


i=l,...,n  (3-48) 


where  COV(Xp)i  is  the  COV  of  the  predicted  strength  due  to  the  uncertainty  in  the  ith 
random  variable  and  n  is  the  number  of  random  variables.  Table  6-4  presents  the  COV's 
of  the  predicted  strength  due  to  each  uncertainty  and  the  parametric  coefficient  of  each 
random  variable  for  Modes  I  and  II.  Figure  6-10  shows  the  contribution  of  each 
uncertainty  to  the  sum  in  equation  3-48,  for  Modes  I.  Figure  6-1 1  shows  the 
corresponding  results  for  Mode  II. 
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Random  Variable 

Mode  I: 
COV(Xp)i 
(%) 

Mode  I: 
Parametric 
coefficient 

Mode  II: 
COV(Xp)i 
(%) 

Mode  II; 
Parametric 
coefficient 

Plate  Thickness 

0.10 

0 

0.58 

0.09 

Depth  of  stiffener 

0.33 

0.03 

0.08 

0 

Web  thickness 

0.01 

0 

0.01 

0 

Flange  thickness 

0.04 

0 

0.01 

0 

Plate  yield  stress 

0.01 

0 

5.39 

0.87 

Stiffener  yield  stress 

4.82 

0.48 

0.01 

0 

Young's  modulus,  plate 

0.19 

0.01 

0.47 

0.08 

Young's  modulus, 
stiffener 

0.39 

0.03 

0.04 

0 

Eccentricity 

8.90 

0.88 

2.86 

0.46 

Residual  stress 

— 

' 

0.76 

0.12 

Table  6-4.  Ranking  of  uncertainties  using  Parametric  Analysis  for  Modes  I  and  II 
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The  following  are  observed: 

•  Mode  I: 

1 .  The  uncertainty  in  eccentricity  is  the  most  important  uncertainty  for  Mode  I. 

2.  The  uncertainty  in  the  yield  stress  of  the  stiffener  is  the  second  most  important 
uncertainty. 

3.  All  other  uncertainties  are  unimportant 

•  Mode  II 

1 .  The  uncertainty  in  the  yield  stress  of  the  plate  is  the  most  important. 

2.  The  uncertainty  in  the  residual  stress,  eccentricity,  Young's  modulus,  and  plate 
thickness  are  also  important. 

3.  All  other  uncertainties  are  unimportant. 

Nikolaidis  et  al.  (1992)  also  found  that  the  uncertainty  in  the  yield  stress  of  the  plating 
was  the  most  important  among  the  uncertainties  in  geometry  and  material  properties,  in 
the  case  of  a  stiffened  panel  failing  in  Mode  II,  for  which  b/t=324,  and  a/b=1.35. 

Sensitivity  coefficients 

We  calculated  the  sensitivity  coefficients  using  equations  3-51  to  3-53. 

Tables  6-5  and  6-6  show  the  sensitivity  coefficients  and  for  Mode  I  and 
Mode  II,  respectively.  Tables  6-7  and  6-8  show  the  sensitivity  coefficients  and  Caaj 
for  Mode  I  and  Mode  II,  respectively.  Tables  6-9  and  6-10  show  the  sensitivity 
coefficients  and  for  Mode  I  and  Mode  II,  respectively. 
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Table  6-5.  Mode  I:  ranking  based  on  sensitivity  of  predicted  strength  to  perturbations 
of  the  means  of  the  random  variables. 


Sm. 

0.791 

0.543 

0.261 

0.075 

0.062 

0.027 

0.019 

0.014 

0.003 

0 


TPL 


HSW 


TSW 


Residual  Stress 


TSF 


E 


E 


SigY(Stf) 


0.864 


0.445 


0.218 


0.063 


0.040 


0.037 


0.032 


0.015 


0.015 


0 


Table  6-6.  Mode  II:  ranking  based  on  sensitivity  of  predicted  strength  to  perturbations 

of  the  means  of  the  random  variables. 
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r.v. 

0.999 

0.021 

0.003 

HSW 

0.002 

TPL 

0 

IBSXSh 

0 

0 

TSF 

0 

TSW 

0 

r.v. 

^oa 

0.871 

0.491 

0.004 

1  HSW 

0.002 

0.001 

TPL 

0 

TSF 

0 

0 

1  TSW 

0 

Table  6-7.  Mode  I :  ranking  based  on  sensitivity  of  predicted  strength  to  perturbations 
of  the  standard  deviations  of  the  random  variables. 


r.v. 

mMm 

0.993 

0.116 

Residual  Stress 

0.006 

TPL 

0.004 

0.003 

HSW 

0.001 

0.001 

TSW 

0 

TSF 

0 

0 

r.v. 

^OG 

Miaasa 

0.968 

0.197 

Residual  Stress 

0.153 

TPL 

0.011 

0.008 

HSW 

0 

0 

0 

TSF 

0 

TSW 

0 

Table  6-8.  Mode  II:  ranking  based  on  sensitivity  of  predicted  strength  to  perturbations 
of  the  standard  deviations  of  the  random  variables. 

Coefficients  C^oj  show  the  effect  of  changing  the  standard  deviations  of  the  random 

variables  to  the  COV  of  the  strength.  It  is  observed  that  the  most  important  uncertainties 
for  Mode  I,  ranked  in  terms  of  importance,  are  in  the  eccentricity  and  yield  stress  of  the 
stiffener.  These  results  are  consistent  with  the  results  from  parametric  analysis.  For  Mode 
II,  the  yield  stress  of  the  plate  is  the  most  important  random  variable  while  the  eccentricity 
and  residual  stress  are  less  important,  although  they  are  not  negligible.  All  other 
uncertainties  are  negligible. 
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r.v. 

0.851 

1  HSW 

0.392 

0.253 

1  TSF 

0.161 

0.129 

1  TPL 

0.106 

0.066 

1  TSW 

0.006 

■Esxaii 

0.001 

r.v. 

Wap 

HSW 

0.729 

0.465 

1  TSF 

0.335 

0.281 

TPL 

0.238 

0.060 

0.031 

1  TSW 

0.008 

■^XSSI 

0 

Table  6-9.  Mode  I:  ranking  based  on  sensitivity  of  predicted  strength  to  weighted 
perturbations  of  the  means  of  the  random  variables. 


r.v. 

TPL 

0.800 

0.412 

0.399 

HSW 

0.138 

Residual  Stress 

0.093 

0.040 

■SSHS 

0.029 

TSF 

0.022 

TSW 

0.004 

0 

r.v.  1 

Wap 

0.840 

HSW 

0.510 

TSW 

0.132 

Residual  Stress 

0.085 

TSF 

0.072 

0.063 

0.030 

0.024 

TPL 

0.008 

BSSIESm 

0 

Table  6-10.  Mode  II:  ranking  based  on  sensitivity  of  predicted  strength  to  weighted 
perturbations  of  the  means  of  the  random  variables. 


A  summary  of  the  most  important  variables  in  each  case  is  given  in  Table  6-11  for 
Mode  I  and  Table  6-12  for  Mode  II.  The  three  most  important  variables  are  listed  in  the 
order  of  impact  on  the  predicted  strength  mean  and  coefficient  of  variation.  If  different 
values  of  the  flange  and  web  dimensions  were  used,  the  rankings  could  be  different. 
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Sensitivity 

Coefficient 

Most  important 
random  variable 

2nd  most  important 
random  variable 

3rd  most  important 
random  variable 

Stiffener  Yield 
Stress 

Stiffener  Web 
Height 

Stiffener  Young's 
Modulus 

Stiffener  Web 
Height 

Stiffener  Young's 
Modulus 

Plate  Thickness 

Eccentricity 

Stiffener  Yield 
Stress 

Stiffener  Young's 
Modulus 

Eccentricity 

Stiffener  Yield 
Stress 

Stiffener  Young's 
Modulus 

Stiffener  Yield 
Stress 

Stiffener  Web 
Height 

Stiffener  Young's 
Modulus 

Wan 

Stiffener  Web 
Height 

Stiffener  Young's 
Modulus 

Stiffener  Flange 
Thickness 

Table  6-11.  Mode  I  summary  of  results  for  sensitivity  coefficient  rankings. 


Sensitivity 

Coefficient 

Most  important 
random  variable 

2nd  most  important 
random  variable 

3rd  most  important 
random  variable 

Plate  Thickness 

Plate  Yield  Stress 

Plate  Young's 
Modulus 

Plate  Yield  Stress 

Plate  Thickness 

Stiffener  Web 
Height 

Cpa 

Eccentricity 

Plate  Yield  Stress 

Residual  Stress 

^oa 

Plate  Yield  Stress 

Eccentricity 

Residual  Stress 

Wpp 

Plate  Thickness 

Plate  Yield  Stress 

Plate  Young's 
Modulus 

Wap 

Plate  Yield  Stress 

Stiffener  Web 
Height 

Stiffener  Web 
Thickness 

Table  6-12.  Mode  II  summary  of  results  for  sensitivity  coefficient  rankings. 
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6.3  Results  for  modeling  uncertainty 

Here  the  objective  was  to  estimate  the  uncertainty  in  predicting  strength  due  to 
idealizations  and  approximations  in  theory.  We  quantified  modeling  uncertainty  by 
evaluating  the  mean  and  standard  deviation  of  modeling  bias  (equation  3-3).  We 
considered  stiffened  panels  that  are  subjected  to  purely  axial  load  because  this  is  the  only 
case  for  which  sufficient  experimental  data  is  available.  We  again  used  the  three  data  sets 
described  in  Section  6.1:  those  of  Faulkner  (1977)  and  ofMichelutti  and  Murray  (1977). 
The  results  are  presented  in  this  section.  Then,  in  Section  6.3.2,  we  compare  our  results 
with  those  of  Faulkner's  own  analysis  of  the  data,  and  with  the  results  of  another  study  by 
Bonello.  Finally,  in  Section  6.3.3,  we  demonstrate  a  Bayesian  method  for  estimating 
modeling  bias  using  both  subjective  information  (estimates  of  experts)  and  objective 
information  (measurements). 

Uncertainty  in  panel  end  rotational  restraint 

To  attempt  to  account  for  uncertainty  in  panel  end  rotational  restraint  would  introduce 
enormous  complexity  into  the  modeling  of  panel  strength.  We  believe  that  this  would  be 
counter  productive  because: 

•  There  is  not  sufficient  experimental  data  available.  To  obtain  such  data  would 
require  extensive  testing  of  a  large  number  of  multi-bay  stiffened  panels  under 
several  loading  conditions. 

•  In  typical  ship  structures  the  cross  frames  are  of  open  section  and  have  relatively 
little  torsional  stiffness.  Therefore  they  provide  only  a  small  amount  of  rotational 
restraint,  and  so  ship  panels  usually  have  little  more  than  simple  support.  With 
sufficient  lateral  pressure  panels  can  behave  as  if  they  were  clamped,  but  obviously 
such  pressure  may  not  always  be  present.  Therefore  the  usual  design  practice  for 
ship  panels  is  to  assume  simple  support,  because  this  idealization  is  the  closest  to 
reality  and  it  lies  on  the  conservative  side. 

•  In  multibay  panels  the  individual  panels  are  approximately  identical  and  are 
therefore  all  approaching  collapse  at  about  the  same  time.  Therefore  support  of 
one  panel  by  another  cannot  be  assumed. 

6.3.1  Estimation  of  modeling  uncertainty 

There  are  only  two  sets  of  data  that  are  sufficiently  complete  to  obtain  predictions  of 
modeling  uncertainty  -  the  parametric  series  and  the  nominally  identical  series  of  Faulkner 
(1977).  We  first  estimated  the  total  bias  and  then  derived  the  modeling  bias  using 
equations  3-5  and  3-6,  together  with  the  estimates  of  random  bias  from  Section  6.2.  We 
did  not  estimate  the  modeling  bias  directly  from  equation  3-3  because  this  equation 
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requires  prediction  of  the  strength  analytically  for  each  panel,  and  Faulkner  did  not 
provide  sufficient  information  on  the  location  of  the  applied  load  to  do  so. 


Parametric  series 


Table  6-2  presents  the  geometry  and  material  properties  of  the  panels.  Two  of  the 
panels  tripped  and  the  others  all  failed  in  Mode  II.  Three  of  the  latter  panels  had 
significantly  different  geometry  than  the  rest;  their  stiffeners  had  no  flange.  Therefore,  we 
only  considered  the  remaining  nineteen  panels. 

Table  6-13  presents  the  nominal  and  measured  strength,  and  the  total  bias  of  each  of 
the  nineteen  panels.  The  theoretical  strengths  are  not  presented.  The  theory  expects  the 
loading  to  be  centered  on  the  elastic  neutral  axis.  As  the  loading  was  applied  midway 
between  the  elastic  neutral  axis  and  the  neutral  axis  at  failure,  consistent  predictions  with 
measurements  were  not  forthcoming.  Attempts  were  made  to  include  loading  alignment  in 
the  panel  eccentricity  with  mixed  success.  Thus,  to  preclude  adding  to  the  modeling 
uncertainty  through  further  approximations,  we  decided  to  use  only  total  bias  results  for 
the  parametric  series  data.  Figure  6-12  shows  the  correlation  between  experiments  and 
nominal  values  of  the  ultimate  stress  divided  by  the  yield  stress.  Table  6-14  presents  the 
nominal  values  of  the  parameters  of  each  panel  that  we  used  to  predict  nominal  strength. 
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Nominal 

strength 

Moden 

(MPa) 

Experimental 
Strength 
Mode  n. 
(MPa) 

Total  bias 

203.71 

255.71 

1.26 

168.31 

186.92 

1.11 

134.81 

121.33 

131.54 

1.08 

199.50 

196.11 

0.98 

164.72 

185.25 

1.12 

131.93 

1— 

118.90 

129.78 

1.09 

188.51  ^ 

179.00 

0.95 

155.42 

162.36  ^ 

1.04 

124.89 

120.54 

0.97 

113.60 

114.24 

1.01 

165.83 

190.24 

1.15 

130.84 

134.28 

1.03 

116.88 

126.50 

1.08 

162.64 

152.85 

0.94 

128.25 

145.25 

1.13 

114.72 

117.39 

1.02 

151.92 

128.24 

0.84 

119.99 

122.26 

1.02 

108.82 

94.46 

0.87 

Table  6-13.  Estimation  of  total  bias  from  parametric  series 
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Notes: 

1.  The  following  assumptions  were  made  regarding  nominal  values: 

a.  Yield  stress  of  stiffener  and  plate  =  240  N/mm^ 

b.  Tension  zone  parameter  due  to  residual  stress  =  2.772. 

c.  The  nominal  eccentricity=a/1500 

2.  Stiffeners  of  panels  13,  17  and  21  have  no  flange 


Table  6-14:  Nominal  values  of  parameters  corresponding  to  material  properties  and 
geometry  of  Panels  in  Parametric  Series 
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The  mean  value  of  the  samples  of  total  bias  in  Table  6-13  is  1.0367,  the  standard 
deviation  0.10  and  the  COV  0.097.  The  standard  deviation  of  mean  total  bias  is  0.023 
which  is  only  about  2%  of  the  mean.  This  means  that  the  sample  size  was  large  enough 
for  estimating  the  mean  bias  accurately.  From  equations  3-5  and  3-7,  we  found  that  the 
mean  value  of  modeling  bias  is  1.05  and  its  COV  0.074.  Therefore,  the  analytical  method 
that  we  used  to  predict  strength  was  conservative,  and,  on  the  average,  underestimates 
strength  by  approximately  5%. 

Nominally  identical  series 

This  set  consisted  of  13  nominally  identical  panels.  Nine  of  these  panels  failed  in 
Mode  II  and  four  of  them  (panels  8-1 1)  failed  in  Mode  I.  Table  6-15  presents  the  stress  at 
failure  and  total  bias  of  the  first  group  of  panels.  The  nominal  failure  stress  was  found  to 
be  equal  to  158.47  N/mm^.  Using  the  results  in  Table  6-15  we  found  that  the  mean  total 
bias  was  1. 1562  and  its  COV  0.06555.  From  equations  (3-5)  and  (3-7),  we  found  that  the 
mean  modeling  bias  is  1 . 17  and  its  COV  0.0195.  We  believe  that  the  COV  of  modeling 
bias  is  unrealistically  low  because  the  panels  were  nominally  identical,  which  means  that 
the  scatter  of  the  ratio  of  experimental  results  vs.  predictions  was  too  low. 


Panel 

Measured 

Collapse  Load  (KN) 

Measured  Stress 
at  Failure  (MPa) 

Total  Bias 

1 

717 

178.5012 

1.126368 

2 

717 

178.5012 

1.126368 

3 

786 

195.6791 

1.234764 

4 

751 

186.9657 

1.179781 

5 

725 

180.4928 

1.138936 

6 

646 

160.8253 

1.104831 

7 

705  1 

175.5137 

1.107517 

12 

804 

200.1603 

1.263041 

13 

773 

192.4427 

1.214341 

Table  6-15.  Experimental  results  and  total  bias  from  nominally  identical  series 

The  COV  of  the  measurements  of  strength  of  the  four  panels  that  failed  in  Mode  II 
was  very  low.  This  is  probably  because,  as  explained  in  Section  6.1.1,  these  panels  were 
initially  very  eccentric  and  they  were  artificially  straightened  by  heating  them.  These 
unusual  statistics  seem  to  confirm  our  decision  not  to  include  these  panels. 
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Michelutti's  results 


Michelutti  and  Murray  (1977)  tested  fifteen  panels,  thirteen  of  which  were  nominally 
identical.  Three  panels  were  axially  loaded,  and  the  rest  were  subjected  to  combined 
bending  and  axial  loading.  Of  the  three  panels  that  were  axially  loaded,  two  failed  in 
Mode  II,  while  the  other  panel  was  designed  to  fail  in  Mode  I.  We  analyzed  only  the 
results  of  the  two  panels  that  failed  in  Mode  II.  The  eccentricity  of  one  panel,  called  A6 
by  Michelutti  and  Murray,  was  equal  to  a/530.  The  eccentricity  of  the  other  panel,  called 
H,  was  not  reported.  We  assumed  that  the  eccentricity  of  panel  H  was  a/1500.  Table  6- 
16  compares  predicted  and  measured  failure  stresses  for  panels  A6  and  H.  It  is  observed 
that  theory  underestimated  the  strength  of  both  panels.  Modeling  bias  was  about  1.05  for 
both  panels. 

Since  there  were  only  two  samples,  we  did  not  try  to  estimate  the  statistics  of  bias 
from  this  table.  However,  we  used  that  data  to  demonstrate  the  Bayesian  method  for 
estimating  bias,  which  will  be  described  in  the  Section  6.3.3. 


Panel 

Predicted  Stress  at 
Failure/Yield  Stress 

Measured  Stress  at 
Failure/Yield  Stress 

Modeling  Bias 

A6 

0.39 

0.41 

1.05 

H 

0.62 

0.66 

1.06 

Table  6-16.  Predicted  and  Measured  Strength  of  Panels  Studied  by  Michelutti 


6.3.2  Comparison  of  estimates  of  bias  with  results  from  other  studies 
Faulkner  09761 

Faulkner  (1976)  evaluated  the  total  and  modeling  bias  on  the  basis  of  the  results  of  the 
nominally  identical  series.  He  analyzed  the  results  of  all  the  thirteen  panels  together, 
regardless  of  the  mode  in  which  they  failed,  and  regardless  of  the  artificial  straghtening  of 
four  of  them,  all  of  which  failed  in  Mode  I.  Also  he  neglected  any  uncertainty  in 
eccentricity.  He  found  the  COV  of  the  total  bias  to  be  0. 109,  and  the  COV  of  the  random 
bias  to  be  0.057.  Then  using  equation  (3-7)  he  found  the  COV  of  the  modeling  bias  to  be 
0.093. 

Our  results,  which  were  based  on  nine  of  these  panels,  all  of  which  failed  in  Mode  II, 
are  considerably  different.  We  found  that  the  COV  of  the  total  bias  is  0.0655,  the  COV  of 
the  random  bias  is  0.0626,  and  the  COV  of  the  modeling  bias  is  0.0195. 

We  believe  that  Faulkner  overestimated  total  bias  because  he  considered  panels  that 
failed  in  Mode  I  and  II  together,  and  that  he  underestimated  random  bias  because  he 
neglected  eccentricity.  Using  these  in  equation  (3-7)  necessarily  gives  an  overestimate  of 
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the  modeling  bias.  Also  it  is  important  to  note  that,  strictly  speaking,  nominally  identical 
panels  are  inappropriate  for  estimating  modeling  bias  and  total  bias,  because  their 
similarity  will  cause  an  underestimate  of  the  CO  V  of  total  and  noimnal  bias.  Therefore, 
the  estimates  of  bias  from  the  parametric  series,  presented  in  Section  6.3.1.,  are  more 

realistic. 

The  results  discussed  in  the  remaining  of  this  subsection  are  for  failure  Mode  II  and  for 
the  case  where  panels  are  subjected  to  axial  loads  only. 

Bonello  et  al  (1992) 

Bonello  et  al  (1992)  developed  a  computer  program  for  predictirig  the  strength  of 
stiffened  multi-span  panels  using  an  inelastic  beam-column  formulation.  They  compared 
predictions  using  this  program  with  a  database  of  measurements  from  85  tests  on  single 
and  multi-span  stiffened  panels.  The  panels  were  subjected  to  axial  compression  only. 
Bonello  et  al.  (1992)  found  that  the  mean  modeling  bias  was  1.02  and  the  COV  of 
modeling  bias  0.08.  They  also  compared  the  computer  predictions  with  some  data  that 
they  called  "numerical  data"  (63  cases),  and  found  a  mean  bias  of  1.05  and  a  COY  of  0.05. 
They  attributed  the  difference  in  the  COV's  to  the  fact  that  the  COV  corresponding  to  t  e 
experiments  included  the  effect  of  scatter  due  to  experimental  errors. 

Bonello  et  al  also  calculated  the  mean  and  COV  of  modeling  bias  for  the  strength 
algorithms  used  by  various  authorities,  based  on  23  tests.  (They  did  not  indicate  whether 
these  tests  were  a  subset  of  the  85  tests  in  the  database).  Table  6-17  presents  the  results. 
The  estimates  from  the  present  study  are  also  included  in  this  Table. 


Reference  where 
method  for  strength 
prediction  is  described 

Mean  Modeling 
Bias 

COV  of  Modeling 
Bias 

Bonello  (1992) 

0.08 

BS5400  (1982) 

1.10 

0.13 

1.10 

0.16 

1.08 

0.14 

ECCS  (1990) 
(Orthotropic  plate 
approach) 

1.14 

0.15 

API  RP2V(1987) 

1.34 

0.34 

Present  study 

1.05 

0.074 

Note;  All  of  the  above  values  are  based  on  the  23  tests  except  for  those  of  the 
Present  study,  which  are  based  on  1 9  tests  from  Faulkner's  parametric  series. 

Table  6-17.  Comparison  of  modeling  uncertainties  of  prediction  methods  used  by 
various  authorities  and  results  from  present  study. 
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Guedes  Soares  and  Soreide  (1983'> 

Guedes  Soares  and  Soreide  (1983)  estimated  the  statistics  of  modeling  bias  for  four 
different  prediction  methods:  Carlsen  (1980),  Dwight  and  Little  (1976),  Faulkner 
(1975B),  and  Horne  and  Narayanan  (1977).  They  obtained  the  experimental  values  from 
three  test  series:  Home  and  Narayanan  (1976),  Horae,  Montague  and  Narayanan  (1977) 
and  Faulkner  (1977)  (of  the  latter  they  only  used  the  results  of  the  parametric  series). 
They  also  compared  predictions  from  Carlsen's  and  Faulkner's  methods  with  data  they 
called  "numerical  data",  which  were  reported  by  Carlsen  (1980). 

Faulkner's  method  (1975B)  is  based  on  a  Johnson-Ostenfeld  approximation  together 
with  the  effective  width  approach  for  plate  behavior.  This  method  was  also  presented  in 
Faulkner  et  al  (1973).  Moreover,  Faukner  (1977)  compared  predictions  of  this  method 
with  results  from  parametric  series  tests.  The  methods  of  Carlsen  (1980),  Dwight  and 
Little  (1976),  and  Horne  and  Narayanan  (1977)  are  based  on  a  Perry-Robertson 
formulation. 

As  it  happens,  Guedes  Soares  and  Soreide  (1983)  chose  to  deal  with  the  inverse  of 
modeling  bias;  that  is,  they  used  the  ratio  of  experimental/predicted  values.  Therefore  we 
repeated  their  analysis  in  order  to  obtain  the  statistics  of  modeling  bias.  Table  6. 18 
presents  these  statistics  for  Faulkner's  and  Carlsen's  methods,  which  were  derived  using 
Faulkner's  parametric  series  (Faulkner  1977).  The  bias  of  the  Standard  algorithm 
estimated  in  Section  6.3.1.  is  also  presented  in  Table  6. 18.  In  order  to  obtain  results  that 
can  be  compared  with  the  estimates  of  bias  of  the  Standard  algorithm  we  used  only  19 
panels  to  derive  the  results  in  Table  6. 18.  Specifically,  we  discarded  two  panels  that  failed 
in  tripping  in  tests  and  three  panels  whose  stiffeners  had  no  flange  (see  Table  6-13).  Table 
6.19  presents  the  modeling  bias  of  the  methods  of  Carlsen  (1980),  Dwight  and  Little 
(1976),  Faulkner  (1975B),  and  Horne  and  Narayanan  (1977)  for  each  of  the  stiffener  and 
plate  induced  collapse  modes.  It  also  presents  the  statistics  of  bias  of  these  methods 
obtained  by  analyzing  data  on  both  failure  modes  together.  The  estimates  in  Table  6. 19 
were  based  on  results  in  Horne  and  Narayanan  (1976)  and  Horae  Montague  and 
Narayanan  (1977). 


Method 

Mean  bias 

bias  COV 

Carlsen  (1980) 

1.18 

0.10 

Faulkner  (1977) 

1.01 

0.07 

Standard 

1.05 

0.07 

Table  6-18:  Modeling  Bias  of  Carlsen's,  Faulkner's  and  the  Standard  methods  for  the 
test  series  of  Faulkner. 
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Stiffener  Induced 
Collapse  (Mode  I) 

Plate  Induced 
Collapse  (Mode  II) 

Modes 

I  and  II 

Method 

Mean 

bias 

bias 

COV 

Mean 

bias 

bias 

COV 

Mean 

bias 

bias 

COV 

Carlsen 

(1980) 

1.04 

0.05 

1.09 

0.08 

1.08 

0.07 

Dwight 
and  Little 
(1976) 

0.98 

0.08 

1.04 

0.07 

1.03 

0.07 

Faulkner 

(1975B) 

0.95 

0.07 

1.00 

0.06 

0.099 

0.07 

Home 

and 

Narayanan 

(■977) 

1.06 

0.06 

1.06 

0.10 

1.06 

0.09 

Table  6-19:  Modeling  Bias  of  Carlsen's,  Dwight  and  Little,  Faulkner  and  Home  and 
Narayanan's  methods  for  the  test  series  of  Horne  et  al. 


It  is  observed  that  the  mean  bias  of  Faulkner's  method  is  closer  to  one  than  that  of  the 
other  methods  and  its  CO  V  is  small,  ranging  from  0.06  to  0.07.  The  results  of  this 
method  correlated  particularly  well  with  the  parametric  series  tests  (Table  6. 18).  On  the 
other  hand  this  method  is  not  conservative;  its  average  bias  is  less  than  one  for  some  tests, 
which  means  that  it  overestimates  strength  in  these  tests.  Carlsen's  method  is  the  most 
conservative.  The  average  bias  of  this  method  for  parametric  series  is  1.18.  The  results  of 
this  method  have  larger  scatter  than  those  of  Faulkner's  method  and  those  of  the  Standard 
method.  Dwight's  and  Little's  method  and  Horne's  method  are  also  conservative,  and  the 
latter  is  slightly  more  conservative  than  the  former.  Both  methods  have  small  scatter  — 
their  COVs  are  smaller  than  0.10. 

Guedes  Soares  and  Soreide  (1983)  also  compared  predictions  of  Faulkner's  and 
Carlsen's  methods  with  numerical  data  by  Carlsen  (1980).  They  found  that  Carlsen  s 
method  correlated  better  than  Faulkner's  method  with  the  numerical  data. 
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6.3.3  Demonstration  of  the  use  of  Bayesian  estimation 

In  Chapter  3  we  described  two  methods  for  using  both  objective  and  subjective 
information  to  estimate  modeling  bias  .  One  method  estimates  the  mean  bias  assuming 
that  the  standard  deviation  of  bias  is  known.  The  other  method  estimates  both  the  mean 
and  standard  deviation  of  bias.  In  this  section  we  provide  a  qualitative  demonstration  of 
both  methods.  Bayesian  estimation  requires  subjective  estimates  of  statistical  properties, 
which  are  normally  obtained  from  a  formal  survey  of  expert  opinion.  Since  this  was 
merely  a  demonstration  we  specified  the  values  ourselves,  and  we  varied  them 
systematically  over  a  wide  range  in  order  to  better  assess  the  effect  of  subjective  estimates 
on  the  statistics  of  bias. 

6.3 .3 . 1  Estimation  of  mean  bias 


We  used  equations  (3-42  a  and  b)  to  estimate  the  posterior  mean  of  bias  and  standard 
deviation  of  mean  bias,  respectively.  We  assumed  that  the  bias  was  normal  and  that  the 
standard  deviation  was  known.  To  demonstrate  the  method,  we  used  experimental  data 
from  the  Parametric  tests  and  from  Michelutti's  tests,  together  with  values  of  subjective 
estimates  of  mean  bias,  which  were  chosen  so  as  to  adequately  cover  the  range  of  such 
estimates. 

Parametric  series 


We  estimated  the  mean  of  total  bias  using  subjective  estimates  and  data  from  the 
parametric  series.  We  varied  the  subjective  estimate  of  mean  bias  (prior  estimate  of  mean 
bias)  from  0.9  to  1.3  and  its  standard  deviation  from  0.025  to  0.3.  The  standard  deviation 
of  total  bias  was  fixed  at  0. 1  (this  value  is  equal  to  the  standard  deviation  of  total  bias 
obtained  from  parametric  series). 

Figure  6-13  presents  the  posterior  mean  bias  as  a  function  of  the  prior  estimate  of 
mean  bias.  Two  curves  corresponding  to  the  mean  +  and  -  one  standard  deviation  of 
mean  bias  are  also  presented.  The  probability  that  the  actual  value  of  mean  bias  is  within 
the  range  defined  by  these  curves  is  68%.  The  standard  deviation  of  the  prior  estimate  of 
mean  bias  was  assumed  to  be  0. 1  in  this  figure.  The  mean  bias  obtained  from  the 
parametric  series  only  is  shown  by  the  dotted  line.  We  observe  the  following  from  Figure 
6-13; 


•  The  posterior  estimate  of  mean  bias  and  the  width  of  the  ±  one  standard  deviation 
range  are  insensitive  to  prior  knowledge. 

•  The  standard  deviation  of  the  estimate  of  mean  bias  from  the  sample,  which  as 
mentioned  in  Section  6.3.1  is  0.023,  is  equal  to  the  standard  deviation  of  the  posterior 
mean  bias. 
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The  reason  that  the  posterior  mean  and  standard  deviation  are  insensitive  to  prior 
knowledge  is  that  the  confidence  in  the  estimate  of  mean  bias  obtained  from  the  samples  is 
high  because  the  sample  size  is  large.  Therefore,  the  posterior  estimate  of  mean  bias  is  not 
affected  by  prior  knowledge. 

Figure  6-14  shows  the  posterior  mean  as  a  function  of  the  standard  deviation  of  the 
prior  mean.  It  is  observed  that  prior  knowledge  affects  the  mean  only  for  very  small 
standard  deviations  in  the  range  from  0.025  to  0. 125.  In  real  life  it  is  very  difficult  to 
guess  the  mean  bias  with  such  high  accuracy.  Therefore,  since  the  sample  size  of 
parametric  series  data  is  large,  and  the  accuracy  of  the  mean  bias  estimated  from  that 
sample  is  high,  we  do  not  need  subjective  estimates  to  estimate  the  mean  bias.  In  cases 
like  this,  where  the  sample  size  is  large,  there  is  no  point  in  using  Bayesian  estimation. 

Finally,  it  is  observed  that  the  standard  deviation  of  the  posterior  estimate  is  lower 
than  both  the  standard  deviation  of  the  sample  mean  and  the  prior  estimate  of  mean  bias. 
This  can  be  explained  by  examining  equation  (3-42b)  in  which  the  standard  deviation  of 
the  posterior  estimate  of  mean,  o'^,  is  less  than  both  the  standard  deviation  of  the  prior 
estimate  of  mean,  Oq,  and  the  standard  deviation  of  bias,  a. 

Michelutti's  data 


In  this  case  we  only  have  two  samples  of  modeling  bias.  The  sample  mean  is  1.058. 
We  assumed  that  the  standard  deviation  of  modeling  bias  was  0.078  (this  value  was 
obtained  from  the  parametric  series) 

Figure  6-15  shows  the  posterior  mean  bias  as  a  function  of  the  prior  estimate  of  this 
value.  The  standard  deviation  of  the  prior  estimate  of  mean  bias  is  0. 1  for  the  results  in 
this  figure.  In  this  case  the  posterior  estimate  of  the  mean  bias  is  affected  significantly  by 
prior  information.  This  happens  because  we  have  only  two  samples  and  therefore  there  is 
high  uncertainty  in  the  mean  bias. 

Figure  6-16  shows  the  posterior  mean  as  a  function  of  the  standard  deviation  of  the 
prior  estimate  of  the  mean.  For  small  standard  deviations,  the  posterior  estimate  of  the 
mean  is  significantly  different  than  the  prior  mean,  which  means  that  prior  knowledge 
affects  significantly  the  posterior  estimate.  On  the  other  hand,  for  high  values  of  the 
standard  deviation  of  prior  mean,  the  posterior  mean  is  almost  equal  to  the  sample  mean. 
This  happens  because  when  the  standard  deviation  of  the  prior  estimate  is  large,  there  is 
little  confidence  in  this  estimate  and  so  the  posterior  mean  is  almost  the  same  as  the 
sample  mean. 

In  conclusion,  in  cases  where  the  sample  size  is  small  (say  less  than  five),  and  there  is 
a  reliable  subjective  estimate  of  mean  bias,  the  Bayesian  approach  is  useful  because  it 
estimates  the  mean  bias  considering  both  subjective  and  objective  information.  The 
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estimate  obtained  from  the  Bayesian  method  can  be  more  reliable  than  the  sample  mean 
bias. 


6.3.3. 2.  Estimation  of  both  mean  and  standard  deviation  of  bias 


This  method  uses  equations  (3-44)  to  (3-46).  We  used  these  equations  to  demonstrate 
how  to  update  prior  knowledge  of  the  mean  and  standard  deviation  of  bias,  which  was 
obtained  from  a  sample  of  size  N,  using  a  new  sample  set.  The  above  equations  do  not 
require  the  samples  used  to  derive  the  prior  estimates  of  mean  and  standard  deviation  but 
only  the  mean,  standard  deviation  and  sample  size.  This  method  is  useful  in  cases  where 
the  original  sample  set,  from  which  prior  knowledge  of  the  mean  and  standard  deviation 
was  obtained,  is  not  available.  This  is  typical  in  real  life  because  in  many  cases  only  the 
statistics  of  bias  and  the  number  of  samples  are  published. 

The  prior  estimates  of  the  mean  and  standard  deviation  of  bias  were  assumed  to  be  1. 1 
and  0.1,  respectively.  We  also  assumed  the  following  five  sample  values  of  bias;  1.0,  0.98, 
0.95,  1.05  and  1.02.  The  mean  of  this  sample  is  1.0  and  the  standard  deviation  is  0.034. 
The  standard  deviation  of  the  mean  is  0.015.  We  calculated  the  posterior  mean  and 
standard  deviation  for  different  values  of  the  sample  size,  N,  on  which  the  prior  estimates 
were  based  in  the  range  from  2  to  100. 

Figures  6-17  and  6-18  show  the  posterior  mean  and  standard  deviation  of  bias, 
respectively,  as  a  function  of  N.  Two  curves  corresponding  to  the  mean  +  and  -  one 
standard  deviation  of  the  mean  are  also  shown.  It  is  observed  that  the  mean  and  standard 
deviation  of  bias  are  sensitive  to  N  for  small  values  of  n  but  for  large  values  they  become 
insensitive  and  converge  to  the  prior  estimates  as  N  becomes  large  (say  greater  than  50). 
Moreover,  Figure  6-17  shows  that,  as  we  expected,  the  confidence  in  the  mean  bias 
increases  as  N  increases. 
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7.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE 

RESEARCH 


7.1  Conclusions 


This  project  has  developed  and  demonstrated  a  methodology  for  estimating  the 
statistics  of  random  and  modeling  uncertainties  in  structural  strength  algorithms.  It  has 
also  presented  a  method  for  determining  the  relative  impact  of  the  basic  strength 
parameters  on  the  random  uncertainty.  Three  different  ranking  methods  were  presented 
with  a  total  of  seven  variations. 

The  methods  of  determining  the  random  and  modeling  uncertainties  were 
demonstrated  for  the  axial  compressive  failure  of  longitudinally  stiffened  panels,  for  both 
Mode  I  (plate-induced)  and  Mode  II  (stiffener-induced)  failure.  The  random  uncertainty 
was  estimated  for  Faulkner's  nominally  identical  series.  For  pure  axial  loading,  the  COV 
of  the  strength  was  found  to  be  10.2%  for  Mode  I  and  6.33%  for  Mode  II  Plots  were 
presented  of  the  interaction  between  the  axial  strength  of  the  panels  and  the  lateral 
pressure  for  Modes  I  and  II,  for  four  different  applied,  transverse  in-plane  loadings.  For 
these  load  combinations  the  random  bias  ranged  from  1.1  to  1.5  for  Mode  I  and  from  1.0 
to  1 . 1  for  Mode  II.  The  CO  Vs  ranged  from  3  %  to  1 5  %  for  Mode  I  and  from  2  %  to  9 
%  for  Mode  II. 

The  following  ranking  was  obtained  for  Faulkner's  nominally  identical  series  of  panels. 
For  stiffener  induced  collapse  (Mode  I)  the  most  important  basic  strength  parameters  are 
the  stiffener  yield  strength,  the  eccentricity  of  the  panel,  and  the  height  of  the  stiffener. 

For  plate  induced  collapse  (Mode  II)  the  most  important  variables  are  the  yield  strength  of 
the  plate,  the  thickness  of  the  plate,  and  the  eccentricity  of  the  panel. 

The  modeling  uncertainties  were  evaluated  for  Mode  II  failure  of  axially  loaded, 
longitudinally  stiffened  panels.  The  standard  algorithm  predictions  for  Faulkner's 
parametric  series  gave  a  modeling  bias  mean  of  1.05  and  a  COV  of  7.4  %.  The  standard 
algorithm  predictions  for  the  Michelutti  panels  also  gave  a  modeling  bias  of  1.05. 

A  comparison  with  other  studies  was  undertaken.  The  results  from  Bonello  (1992) 
gave  a  mean  modeling  bias  of  1 .02  and  a  COV  of  8  %.  In  a  comparison  with  other 
prediction  algorithms  for  a  different  experimental  sample  set,  Bonello  reported  mean 
modeling  biases  ranging  from  1.05  to  1.34.  Guedes  Soares  and  Soreide  (1992)  gave 
results  for  Faulkner's  parametric  series.  The  mean  modeling  bias  for  several  strength 
prediction  methods  ranged  from  1.01  to  1.18  and  the  COVs  ranged  from  7  to  10  %.  A 
variety  of  algorithm  predictions  for  other  experimental  data  are  reported.  For  Mode  I  and 
Mode  II,  the  mean  modeling  bias  ranges  from  0.99  to  1.08  and  the  COVs  ranged  from  6 
to  10  %. 
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A  Bayesian  method  for  estimation  of  the  mean  and  COV  of  modeling  bias  was 
developed  and  demonstrated.  This  method  uses  both  subjective  information  and  results 
from  comparison  of  predictions  with  tests.  The  method  is  useful  in  cases  where  limited 
experimental  data  is  available. 

The  above  results  help  gain  insights  into  levels  of  uncertainty  possible  in  structural 
strength  design  and  analysis,  whether  for  use  deterministically  or  in  a  reliability-based 
format.  However,  the  accuracy  of  the  estimates  of  random  and  modeling  uncertainty  is 
limited  by  the  availability  of  appropriate  experimental  data.  Moreover,  the  findings 
discussed  above  are  limited  to  the  specific  algorithms  and  data  sets  used  in  the  analysis. 

7.2  Recommendations  for  future  research 

7.2.1  Estimation  of  modeling  bias  for  other  failure  modes 

The  assessment  of  uncertainty  in  strength  models  can  only  be  done  when  sufficient 
experimental  data  is  available.  As  shown  in  Table  4. 1  and  as  discussed  in  Section  4.2,  of 
the  ten  failure  modes  that  involve  modeling  uncertainty,  only  four  have  sufficient  data 
available.  Of  these,  the  two  plate  failure  modes  are  merely  unserviceability  rather  than 
collapse,  and  they  are  also  very  local.  In  contrast,  of  the  six  failure  modes  for  which  data 
is  lacking,  all  are  of  the  collapse  type.  Moreover,  they  all  involve  large  portions  of 
structure,  which  makes  them  doubly  important.  Therefore  the  obvious  need  is  for 
sufficient  structural  testing  to  be  done  for  these  six  other  failure  modes. 

Of  these  six,  the  second,  third  and  fourth  are  the  most  common  for  ship  structures: 

•  Panel  collapse  due  to  stiffener  buckling,  which  can  be  either: 

-  local  buckling  of  the  web  or  flange  (easily  avoided  by  using  standard 
sections) 

-  "tripping"  (torsional  buckling  due  to  axial  compression) 

-  flexural-torsional  buckling  (torsional  buckling  due  to  bending  of  the 
stiffener,  caused  usually  by  lateral  loads  on  the  panel) 

•  Beam  collapse  due  to  tripping 

•  Beam  collapse  due  to  flexural-torsional  buckling 

These  three  have  some  features  in  common:  they  all  involve  buckling,  and  the  two 
types  of  torsional  buckling  can  occur  for  both  stiffeners  and  beams.  Therefore  we  believe 
that  the  biggest  need  in  regard  to  structural  uncertainty  is  a  carefully  planned  and 
coordinated  series  of  tests  involving  these  three  failure  modes.  Using  the  results  of  these 
tests  and  the  methodology  for  estimating  modeling  uncertainties  that  was  presented  in  this 
report  we  can  estimate  the  statistics  of  modeling  bias  of  selected  strength  analysis 
algorithms. 
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Because  tripping  is  one  of  the  failure  modes  that  we  regard  as  most  in  need  of 
attention,  we  devoted  some  extra  time  and  effort  in  this  project  to  review  the  current 
approaches  to  tripping.  This  review  is  presented  in  Appendix  A,  and  we  hope  that  it  will 
encourage  and  facilitate  the  experimental  research  that  is  required  in  this  area. 


7.2.2  Analysis  of  random  uncertainties 

The  procedure  presented  in  Section  6.1  needs  to  be  applied  to  a  series  of  panels  that 
have  different  dimensions.  These  panels  should  be  selected  in  a  way  that  their  dimensions 
cover  entire  range  of  values  encountered  in  practical  applications.  Then  the  results  should 
be  carefully  analyzed  using  the  methodology  presented  in  this  report  to  identify  the  most 
important  uncertainties  and  estimate  their  effect  on  strength.  It  is  hoped  that  this  will 
allow  general  conclusions  to  be  drawn  that  are  applicable  to  the  majority  of  panels  used  in 
ship  structures. 

The  above  procedure  should  be  also  applied  to  other  components  failing  under  the 
most  important  modes  in  Table  4.1,  and  also  to  the  entire  ship  hull  failing  under  global 
bending  to  study  random  uncertainties  in  these  modes. 

Finally,  the  statistical  properties  of  statistical  parameters  and  the  type  of  distribution  of 
random  and  modeling  bias  should  be  presented.  This  information  is  important  in  both 
reliability  analysis  and  in  developing  a  reliability  based  design  code  for  ships. 
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Figure  3.1.  Bayes' Theorem 


e  6-1.  Cross-sectional  View  of  a  Panel  in  Nominally  Identical  Series 


Figure  6“2.  Oqsc  A,  Mode  I  Failure  (Punel  with  Negative  Eccentricity).* 
Effect  of  Random  Uncertainties  on  Strength. 
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Figure  6-8.  Case  C,  Mode  II  Failure  (Panel  with  Positive  Eccentricity): 
Effect  of  Random  Uncertainties  on  Strength, 
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Figure  6-11.  Mode  II,  Parametric  Analysis:  Contributions  of  Random 
Uncertainties  to  Uncertainty  in  Strength 
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Figure  6-13.  Bayesian  Estimation  of  Mean  Bias  Using  Results  from 

Parametric  Series:  Posterior  Mean  vs.  Prior  Estimate  of  Mean 
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Figure  6-14.  Bayesian  Estimation  of  Mean  Bias  Using  Results  from 
Parametric  Series:  Posterior  Mean  vs.  Prior  Estimate  of 
Standard  Deviation  of  Mean 
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Figure  6-15.  Bayesian  Estimation  of  Mean  Bias  Using  Michelutti's  Results: 
Posterior  Mean  vs.  Prior  Estimate  of  Mean 
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Figure  6-16.  Bayesian  Estimation  of  Mean  Bias  Using  Michelutti's  Results; 

Posterior  Mean  vs.  Prior  Estimate  of  Standard  Deviation  of 
Mean 
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Figure  6-17.  Bayesian  Method  for  Mean  and  Standard  Deviation  of  Bias: 
Posterior  Mean  Bias  vs.  Number  of  Samples 
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APPENDIX  A:  REVIEW  OF  STIFFENER  TRIPPING 


A  1  Basic  Theory  of  Elastic  Tripping 

When  a  stiffener  is  subjected  to  axial  compression  there  are  a  number  of  ways  for  it  to 
fail,  including  tripping.  Tripping,  or  torsional  buckling,  differs  from  the  buckling  of  a 
column  in  three  ways:  first,  because  the  rotation  occurs  about  the  line  of  attachment 
between  the  plating  and  the  stiffener;  second,  because  the  plating  offers  some  resistance  to 
the  rotation;  and  third,  because  the  stiffener  itself  may  deflect  and  distort  -  it  isn  t  rigid 
body  rotation.  For  the  case  of  typical  ship  structures  the  dimensions  of  the  beams  are  such 
that  we  can  consider  the  problem  to  be  one  of  torsion  of  thin-walled  members.  As  such, 
the  governing  differential  equation  for  the  rotation  ^  is  (Hughes,  1988,  also  referred  to  as 
SSD) 


EIszd^-^-(GJ-aaIsp)^  +  K^<t>  =  0  (A-1) 

dx  dx 

where 

Igz  =  moment  of  inertia  of  the  stiffener  only  about  an  axis  through  the 
centroid  of  the  stiffener  and  parallel  to  the  web. 

2d  =  stiffener  web  height. 

Isp  =  polar  moment  of  inertia  of  the  stiffener  about  the  center  of  rotation. 

K(t,  =  distributed  rotational  restrain  which  the  plating  exerts  on  the 
stiffener. 

If  the  ends  of  the  stiffener  are  regarded  as  simply  supported,  the  solution  for  (t)(x)  is  a 
buckled  shape  in  which  the  rotation  (j)  varies  sinusoidally  in  m  half-waves  over  the  len^h  a 
of  the  stiffener.  The  value  of  the  applied  in-plane  stress  which  causes  tripping  according 
to  elastic  theory  is  denoted  as  cTa,T-  The  solution  to  the  above  differential  equation  is  the 
minimum  value  for  the  applied  in-plane  stress  Oa?  that  satisfies  the  following 


Elszd^  +  (GJ  -  +K4, (Oa,  m)  =  0  (A-2) 

a  ^ 

The  critical  mode  for  tripping  corresponds  to  whichever  integer  value  of  m  gives  the 
minimum  value  of  Ca.  Note  that  the  expression  for  K<j,  is  written  as  a  function  of  both  the 
number  of  half-waves  m  and  the  critical  stress  Oa-  It  is  dependent  on  the  applied  stress 
because  if  the  stress  is  high  enough  to  cause  local  plate  buckling,  the  value  of  K,j,  can  go 
to  zero  or  even  become  negative.  The  dependence  on  m  is  because  the  amount  the  value 
of  K,j)  diminishes  with  Oa  depends  on  whether  the  number  of  half-waves  in  the  stiffener 
matches  the  number  in  the  buckled  plating.  The  other  major  contributor  to  the  rotational 
restraint  provided  by  the  plating  is  the  flexural  rigidity  of  the  plating. 
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The  above  discussion  serves  to  show  that  the  tripping  of  a  stiffener  involves  three 
important  variables  which  interrelate  in  a  rather  complex  fashion:  Oa ,  m ,  and  K<j,.  The 
critical  value  of  Ga  depends  on  K<j,  and  m ,  the  value  of  m  depends  on  the  magnitude  of  K<j, 
relative  to  EIsz ,  and  K,;,  depends  on  both  Oa  and  m. 

A.2  AISC  Approach  to  Tripping 

The  AISC  Code  looks  at  tripping  as  part  of  the  special  case  of  I-shaped  beams 
subjected  to  strong  axis  bending.  This  sort  of  loading  is  considered  in  the  marine  industry 
to  be  a  case  of flexural-torsional  buckling  rather  than  lateral  torsional  buckling.  In  the 
flexure  case  the  flange  of  the  beam  is  put  in  compression  due  to  the  out-of-plane  load 
causing  a  large  bending  moment.  This  is  of  great  concern  to  bridge  and  building  structural 
engineers  because  for  these  structures  the  predominant  loads  are  out-of-plane  loads  which 
produce  bending.  They  do  not  generally  experience  the  very  large  in-plane  types  of  loads 
that  a  ship's  designer  must  be  concerned  with.  However,  the  basic  analysis  procedure  is 
the  same  and  the  LRFD  approach  adopted  by  the  AISC  should  be  carefully  considered. 

The  strength  requirement  for  tripping  is  considered  to  be  a  part  of  the  overall  strength 
requirement  of  the  stability  of  I-shaped  beams.  The  AISC  procedure  considers  the  full 
range  of  strength  from  laterally  stable  beams  to  situations  where  tripping  causes 
considerable  strength  reduction.  According  to  LRFD-F2  ,  the  strength  requirement  may 
be  stated  as 


‘t>bMn  ^  Mu 


(A-3) 


where 

(l)b  =  strength  reduction  factor  for  flexure  =  0.90 
Mn  =  nominal  moment  strength. 

Mu  =  factored  service  load,  typically  given  as  a  linear  combination  of  Dead 
Load,  Live  Load  and  Wind  or  Snow  Load,  e.g.  1.2D  +  1.6L  +  0.5S. 


Figure  A.  1,  taken  from  Salmon  and  Johnson  (1990),  shows  the  effect  of  laterally 
unbraced  length  Lb  on  the  flexural  torsional  buckling  strength.  The  cases  indicated  as 
number  1  and  2  are  not  of  direct  interest.  They  represent  situations  where  the  beam 
reaches  its  full  plastic  moment,  and  other  methods  must  be  used  for  analysis.  The  last 
three  cases  are  of  more  concern. 
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Figure  A.  1  AISC  Nominal  Strength  of  Compact  Sections 


Cases  3  and  4  are  for  the  situation  where  the  nominal  moment  strength,  Mn,  occurs  in 
the  inelastic  region.  The  difference  between  the  two  cases  is  that  case  3  is  for  "compact" 
sections.  Most  rolled  shapes  are  compact  sections,  or  sections  whose  flange  dimensions 
preclude  local  buckling  of  the  compression  flange. 

Of  more  concern,  particularly  for  a  comparative  analysis  is  case  5.  This  is  the  case 
where  the  nominal  moment  strength  Mn  is  the  elastic  buckling  strength.  Here  the  value  of 
the  nominal  moment  strength  is  given  as 


Mjj  -  Mgj  -  C], 

Lb 


n  E 


V  Lb  y 


Cwfy  +EIyGJ 


where 

Cb  =  factor  to  account  for  moment  gradient. 

Lb  =  laterally  unbraced  length. 

Cw  =  torsional  warping  constant  =  If  d2/2. 

If  =  moment  of  inertia  of  one  flange  about  the  y-axis. 
ly  =  moment  of  inertia  about  the  y-axis. 


(A-4) 


The  major  difference  between  the  AISC  formula  and  the  case  for  ships  is  that  typical 
ship  structure  is  stiffened  plating.  With  stiffened  plating  the  restraint  caused  by  the  plating 
must  be  included  in  the  analysis,  and  for  this  reason  the  expression  for  M^  in  equation  (A- 
4)  is  not  applicable.  However,  the  form  of  the  strength  requirement  in  equation  (A-3)  is  of 
interest.  The  AISC  LRFD  has  cast  the  basic  equation  in  the  form  which  is  consistent 
across  the  spectrum  of  column  collapse  mechanisms.  The  LRFD  format  allows  for  a 
constant  load  reduction  factor  and  a  consistent  means  for  applying  loading  (factored 
loads).  Using  the  combination  of  the  two  factors  allows  the  code  developers  to  control 
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the  level  of  safety  implied  in  the  equations.  What  is  hidden  is  all  of  the  work  that  went 
into  calibration  of  the  code  for  the  types  of  variables  and  their  levels  of  uncertainty. 

A.3  AASHTO  Approach  to  Tripping 

The  AASHTO  approach  to  tripping  analysis  is  veiy  simple.  It  consists  of  using  the 
AISC  Working  Stress  Design  approach  for  what  naval  architects  consider  to  be  flexural- 
torsional  buckling.  The  AASHTO  LRFD  code  is  expected  to  be  very  much  like  the  AISC 
LRFD  except  that  the  load  factors  will  be  more  specific  to  the  typical  highway  bridge 
loading  (i.e.  average  daily  truck  traffic).  Because  there  is  little  difference  with  the  AISC 
code  and  due  to  the  non-applicability  to  a  majority  of  typical  ship  structures,  the 
AASHTO  code  will  not  be  considered  fiirther. 

A.4  API  Approach  to  Tripping 

The  American  Petroleum  Institute  guidelines  does  much  the  same  as  the  AASHTO 
code  does  in  regards  to  tripping;  it  defers  to  the  AISC  approach.  However,  the  API 
guidelines  specifically  state  not  to  use  the  AISC  LRFD  approach  because  the  load  and 
resistance  factors  have  been  calibrated  for  typical  building  practice,  not  for  the  harsher 
offshore  environment. 

A.5  U  S.  Navy  Design  Data  Sheets  Approach  to  Tripping 

The  U.S.  Navy  DDS  100-4  addresses  lateral  torsional  buckling  in  terms  of  providing 
support  for  the  stiffener  flange  at  some  intermediate  points  along  its  span.  There  is  no 
specific  design  check  on  the  likelihood  of  tripping.  Figure  A.2  is  taken  from  DDS  100-4 
to  show  the  relationship  between  material  type,  flange  breadth,  and  stiffener  depth. 
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A.  MAXIMUM  permissible  RATIO  OF  SPAN  TO  FLANGE 
WIDTH  FOR  TEE  STIFFENER  WITHOUT  NTERMEDIATE 
LATERAL  SUPPORT 


RATIO  OF  FLANGE  WIDTH  TODEPTHOF  TEE  STIFFENER, t.^/d 

Figure  A.  2  Maximum  Permissible  Ratio  of  Span  to  Flange  Width  for  Tee 
Stiffener  without  Lateral  Support  (DDS  100-4) 

If  the  actual  stiffener  span  to  flange  breadth  ratio  is  greater  than  the  value  of  K3  found 
in  figure  A.2  but  less  than  1.75  then  a  support  must  be  provided  to  the  flange  at  the 
midpoint  of  the  span.  If  the  ratio  is  greater  than  1.75  K3,  support  must  be  provided  at 
distances  of  not  more  than  0.75  Kabf.  This  is  all  based  on  the  assumption  that  the  ratio  of 
the  stiffener  web  thickness  to  flange  thickness  is  about  0.6  and  that  the  sections  are 
compact. 

The  problem  with  an  approach  like  this  is  that  there  is  no  way  to  incorporate  different 
geometries  or  materials.  The  designer  is  stuck  with  a  cook  book  approach  which  offers  no 
means  of  allowing  innovation.  While  this  approach  has  thus  far  produced  safe  designs,  no 
one  really  knows  what  level  of  safety  is  provided. 

A.  6  ABS  Approach  to  Tripping 

The  American  Bureau  of  Shipping  1991  Rules  for  Building  and  Classing  Steel  Vessels 
has  an  even  more  basic  approach  to  preventing  tripping.  In  Section  9.3.5,  Web  Frames  - 
Stiffener  and  Tripping  Brackets,  the  rules  require  the  flanges  of  the  deep  webs  to  be 
supported  at  intervals  of  "about"  10  feet  when  the  flange  breadth  on  either  side  of  the  web 
exceed  8  inches.  The  same  requirement  appears  in  Section  9.5.3,  Side  Stringers,  and 
Section  1 1.7.4,  Deck  Girders.  The  only  place  in  the  Rules  where  those  values  differ 
significantly  is  in  Section  29B.  15.  lb  for  slanted  frames  on  Ice  Class  1 AA  ships.  Then  the 
spacing  is  to  be  not  more  than  5 1  inches. 
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Again,  no  specific  accounting  for  tripping  is  provided.  This  is  probably  due  to  ABS's 
experience  with  typical  ship  structures  which  has  shown  tripping  not  to  be  a  failure 
problem.  But  once  again,  there  is  no  flexibility  in  these  rules  for  innovation  in  concepts  or 
materials.  There  is  also  no  means  of  determining  the  level  of  safety  implied. 

A.  7  Proposed  Model  for  Tripping 

In  order  to  develop  a  probabilistic  approach  to  ship  structural  design,  an  explicit  model 
needs  to  be  developed  for  each  failure  mode.  These  models  can  then  be  used  as  limit- 
states  to  determine  the  levels  of  safety  implied.  The  AISC  LRFD  approach  to  developing 
a  design  code  is  a  sound  one.  The  limit  state  can  be  expressed  in  terms  of  moments  or 
stresses  and  can  be  a  simple  equation  similar  to  equation  (A-3).  In  the  case  of  ship 
structures  the  proposed  approach  is  to  have  the  limit  state  expressed  in  terms  of  applied 
in-plane  stress.  Then  one  only  needs  to  compare  the  critical  buckling  stress  for  tripping  to 
the  applied  stress. 

SSD  presents  a  model  for  tripping  that  consists  of  a  two  step  solution.  Because  of  the 
interdependence  of  the  three  important  variables  as  shown  in  Section  A.  1  of  this  report,  m 
the  number  of  half-waves  in  which  the  stiffener  will  trip,  needs  to  be  determined.  A  (non¬ 
integer)  approximate  value  for  m,  denoted  as  r,  is  given  by 


a  4DCr 

r  =  -4 - ^ 

’^VEIszd  b 


(A-5) 


where 

a  =  span  length  of  the  stiffener, 
b  =  breadth  of  the  stiffener  flange 
D  =  plate  flexural  rigidity,  Et3/(12(l-u2)). 

Cr  =  factor  by  which  the  plate  rotational  restraint  is  reduced  by  web 
bending  which  is  given  by 


Cr  = 


1 


d 

b 


t  =  thickness  of  the  plating, 
tw  =  thickness  of  the  stiffener  web 


For  each  of  the  integer  values  of  m  above  and  below  r,  the  following  equation  is  used 
for  determining  Oaj. 
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^a.T  = 
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0  4DC, 
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2 

m 


where 

Isp  =  polar  moment  ofinertia  of  the  stiffener. 

When  the  value  for  r  found  from  equation  (A-5)  is  very  close  to  1, 
equation  (A-6)  can  be  used 


^a,T  =  ■ 


fsp 


2C,b^l 


a 


n  b 


(A-7) 


11  (A-6) 


a  simpler  form  of 
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APPENDIX  B:  EXPLANATION  AND  VALIDATION  OF  THE 
STANDARD  ALGORITHM  FOR  PANEL  COLLAPSE 
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APPENDIX  B:  EXPLANATION  AND  VALIDATION  OF  THE 
STANDARD  ALGORITHM  FOR  PANEL  COLLAPSE 


(Chapter  14  of  Ship  Structural  Design) 
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ULTIMATE  STRENGTH  OF  STIFFENED  PANELS 


As  explained  at  the  beginning  of  Chapter  12.  the  stiff¬ 
ened  panels  in  ship  structures  are  usually  sufficiently 
sturdy  that  (orj^  >  ay.  Hence  the  mode  of  com¬ 
pressive  collapse  is  not  elastic  panel  buckling,  but 
rather  a  more  complicated  inelastic  process,  for  which 
we  use  the  words  “failure”  or  “collapse”  instead  of 
buckling.  Nevertheless,  there  are  several  aspects  of 
elastic  buckling  which  have  their  counterpart  in  the 
inelastic  compressive  collapse  of  plating.  For  exam¬ 
ple,  in  parallel  with  elastic  buckling  there  are  several 
possible  levels  and  modes  of  inelastic  plate  collapse.  If 
the  panel  is  transversely  stiffened,  then  plate  failure  by 
itself  constitutes  collapse,  and  the  theory  and  methods 
of  Section  12.7  are  sufficient.  If  the  panel  is  cross- 
stiffened  there  are  two  possible  levels  of  collapse: 

1,  Interframe  collapse,  that  is,  collapse  of  stiffened 
panels  between  transverse  frames. 

2.  Gross  panel  collapse,  involving  both  longitudinal 
and  transverse  stiffeners. 

In  general,  cross-stiffened  panels  should  be  designed 
such  that  interframe  collapse  occurs  before  gross  panel 
collapse  because  the  latter  involves  a  larger  portion  of 
structure  and  is  likely  to  be  more  catastrophic.  Also, 
the  structural  proportions  of  ship  panels  are  usually 
such  that  gross  panel  collapse  is  inelastic,  and  there¬ 
fore  it  is  an  extremely  difficult  and  complex  task  to 
calculate  the  collapse  load  with  sufficient  accuracy  for 
design  purposes.  Therefore  the  best  procedure  is  to 
design  the  panel  such  that  gross  panel  collapse  cannot 
occur  before  interframe  collapse.  This  can  be  achieved 
by  ensuring  that  the  elastic  gross  panel  buckling  stress 
calculated  by  the  discrete  beam  method  of  Section 
13.6  is  substantially  larger  than  the  elastic  interframe 


buckling  stress  given  by  (13.1.8),  and  that  both  of 
these  buckling  stresses  are  substantially  larger  than  the 
yield  stress. 

In  practice  the  scantlings  that  a  panel  must  have  in 
order  to  satisfy  the  requirements  of  other  modes  of 
failure  are  usually  sufficient  to  meet  both  of  these 
requirements,  but  they  should  always  be  checked 
nonetheless.  With  this  design  philosophy  it  is  not  nec¬ 
essary  to  consider  inelastic  gross  panel  collapse  and 
therefore  this  chapter  deals  only  with  inelastic  inter- 
frame  collapse. 


14.1  BOUNDARY  CONDITIONS,  LOAD 
TYPES,  AND  COLLAPSE  MODES 

Boundary  Conditions  for  Litimate  Strength 
Analysis 

In  a  deck  or  side  it  is  clear  that  the  rotational  restraint 
provided  to  the  panels  by  the  transverse  trames  is 
relatively  small  and  is  best  ignored.  In  the  bottom  and 
in  other  regions  having  a  large  lateral  load  each  panel 
receives  some  rotational  restraint  from  the  adjacent 
panel,  and  in  principle  this  would  give  a  larger  ulti¬ 
mate  compressive  strength  than  if  the  panel  were  sim¬ 
ply  supported.  However,  in  most  cases  it  is  possible 
for  the  compressive  load  to  occur  with  little  or  no 
lateral  load.  Also  there  are  many  situations  in  which  a 
bottom  panel  is  far  from  being  clamped.  Figure  14.1^ 
shows  a  common  example,  in  a  tanker  the  end  panel  of 
an  empty  cargo  tank  adjacent  to  a  full  cargo  tank  has 
only  moderate  rotational  restraint  at  one  end  and  is 
subjected  to  a  destabilizing  moment  at  the  other  end. 
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PANELS 
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Figure  14.1  Occurrence  of  nonciamped  conditions  in  bottom 
panels. 

Although  other  panels  may  have  more  rotational  re¬ 
straint,  it  is  obvious  that  this  panel  will  govern  the 
design  of  the  bottom  structure;  it  would  be  un¬ 
economical  and  impractical  for  the  panels  to  differ 
from  one  frame  to  another.  In  other  ship  types  the 
same  situation  arises  if  a  cargo  hold  is  flowed  due  to 
underwater  damage.  It  also  arises  in  a  double  bottom 
in  which  one  compartment  is  filled  with  fuel  or  ballast 
and  the  adjacent  compartment  is  empty.  Therefore  un¬ 
less  a  panel’s  rotational  restraint  is  guaranteed  to  be 
permanently  present  under  aU  conditions  (which  is 
quite  rare)  it  should  not  be  counted  when  calculating 
the  ultimate  strength  of  the  panel.  In  general  the  safest 
and  best  procedure  is  to  regard  the  panel  as  simply 
supported.  However,  the  other  effects  of  the  lateri 
load — the  maximum  deflection  and  the  maximum 
bending  moment — should  not  be  ignored  since  they 
decrease  the  panel’s  ultimate  compressive  strength.  In 


this  chapter  the  deflection  and  the  bending  moment 
due  to  the  lateral  load  (acting  alone)  are  denoted  by  So 
and  A/q.  Although  the  panel  ends  are  being  taken  as 
simply  supported  the  theory  presented  in  this  chapter 
does  not  require  that  8o  and  Mq  be  the  simply  supported 
values.  They  can  be  either  these  values  or,  if  desired 
and  if  it  is  practicable  to  calculate  them,  they  can  be 
the  largest  values  which  will  occur,  allowing  for  the 
actual  amount  of  rotational  restraint  that  is  present  in 
that  worst  case. 

Although  simply  supported  boundary  conditions 
generally  correspond  to  small  or  moderate  values  of 
lateral  load,  this  chapter  covers  the  simply  supix)rted 
case  over  the  full  range  of  lateral  loads,  from  zero  up 
to  the  plastic  hinge  value,  and  acting  on  either  side  of 
the  plating. 

Basic  Load  Types  and  Associated  Mechanisms  of 
Collapse 

For  longitudinally  stiffened  panels  with  simply  sup¬ 
ported  edges  there  are  three  basic  types  of  loads: 

1 .  Lateral  load  causing  negative  bending  of  the  panel 
(i.e.,  of  the  plate-stiffener  combination). 

2.  Lateral  load  causing  positive  bending  of  the  panel. 

3.  In -plane  compression. 

Throughout  this  chapter,  bending  moment  in  the  panel 
is  positive  when  it  causes  compression  in  the  plating, 
and  in-plane  loads  are  positive  if  compressive. 

Each  of  these  three  basic  loads  has  associated  with 
it  one  or  more  possible  mechanisms  of  collapse. 
Whenever  in-plane  compression  occurs  in  combina¬ 
tion  with  one  of  the  first  two  load  types  there  is  inter¬ 
action  between  the  specific  mechanisms  of  collapse. 
As  we  shall  see.  the  various  interactions  give  rise  to 
three  distinct  modes  of  collapse.  The  specific  collapse 
mechanisms  associated  with  each  of  the  basic  load 
types  follow  (see  Fig.  14.2). 

Negative  Bending 

With  this  type  of  load  the  plating  is  in  tension  and 
hence  cannot  buckle.  Also,  since  the  neutral  axis  is 
close  to  the  plating,  it  is  the  stiffener  flange  which  first 
reaches  yield  (point  A  in  Fig.  14.2).  As  the  bending 
moment  increases,  yielding  spreads  through  the  stiff¬ 
ener.  The  plating  then  begins  to  yield  in  tension  and, 
if  there  is  no  instabilit}'  in  the  stiffener,  collapse  even¬ 
tually  occurs  when  the  bending  moment  reaches  Mp, 
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the  value  at  which  a  plastic  hinge  is  formed  (see  Sec¬ 
tion  16.1).  This  corresponds  to  point  B  in  Fig.  14.2. 
However,  since  the  stiffener  is  in  compression  it  is 
possible  that  stiffener  buckling  might  occur.  Thus,  for 
a  negative  bending  load  the  plate-stiffener  combina¬ 
tion  may  not  be  capable  of  withstanding  the  full  value 
of  Mp  for  the  section  (point  C  in  Fig,  14.2). 

Positive  Bending 

With  this  type  of  load  the  plating  is  in  compression, 
but  since  the  neutral  axis  is  very  close  to  the  plating  the 
plate  will  not  be  heavily  stressed  until  the  load  be¬ 
comes  very  large,  by  which  time  the  entire  stiffener 
has  yielded  in  tension  and  collapse  is  imminent.  The 
plate  then  begins  to  yield  in  compression,  and  collapse 
finally  occurs  as  the  result  of  a  plastic  hinge,  with  the 
plate  reaching  compressive  yield.  In  fact,  for  most 
ship  panels  the  plastic  neutral  axis  lies  within  the  plat¬ 


ing  thickness,  and  so  some  portion  of  the  plate  s  thick¬ 
ness  is  in  tension.  This  stabilizes  the  plate  and  prevents 
it  from  buckling  and  hence  for  most  ship  panels  the 
plate -stiffener  combination  is  capable  of  absorbing  al¬ 
most  the  full  value  of  the  plastic  bending  moment  Mp 
(point  D  in  Fig.  14.2). 


Uniform  In-plane  Compression 

By  “uniform  compression”  we  mean  an  applied  load 
which  causes  a  uniform  shortening  of  the  panel,  and  it 
is  assumed  that  the  load  is  not  (directly)  influenced  by 
the  panel’s  response.  The  most  common  example  is 
the  hull  girder  stress.  Strictly  speaking,  the  loading  on 
the  panel  is  an  applied  uniform  strain  Sa^  but  for  the 
present  we  shall  deal  with  the  average  applied  stress 
CTa  (or  average  distributed  load  AO- 

Under  this  type  of  load  the  panel  is  essentially  a 
group  of  identical  (and  almost  independent)  un- 
symmctric  columns,  each  consisting  of  a  stiffener  and 
a  plate  flange  of  effective  width  Z?,.  If  the  panel  is 
sufficiently  slender  these  columns  will  buckle  elas¬ 
tically  and  the  ultimate  strength  of  the  panel  can  be 
calculated  by  the  methods  of  Chapter  13.  But  for  typ¬ 
ical  ship  panels  the  plate’s  response  becomes  inelastic 
prior  to  collapse  because  the  plate’s  load  versus  end 
shortening  curve  is  the  cTa  versus  Sa  curves  of  Fig. 
12.24.  This  makes  the  collapse  mechanism  decidedly 
more  complicated,  but  a  method  is  presented  in  this 
chapter  which  overcomes  the  complications  and  gives 
an  accurate  prediction  of  the  collapse  value  of  this  type 
of  load.  Basically  the  method  consists  in  regarding  the 
plate  portion  of  die  column  as  being  made  of  a  differ¬ 
ent  material,  having  a  reduced  elastic  modulus,  and 
then  using  the  basic  column  approach  of  Section  11.2 
to  calculate  the  collapse  load,  making  careful  allow¬ 
ance  for  initial  eccentricity.  In  this  approach, 
“collapse”  is  considered  to  occur  when  the  stress  in  the 
extreme  fiber  of  the  compression  flange  of  the  column 
reaches  the  failure  value  for  the  material  of  that  flange. 
For  the  stiffener  flange  this  is  simply  the  yield  stress, 
but  for  the  plate  flange  it  is  some  other  value  which 
allows  for  the  nonlinearity  of  the  cTa  versus  Sa  curve. 
Thus  there  are  two  different  values  of  in-plane  collapse 
load,  depending  on  which  flange  is  the  cornpression 
flange,  that  is,  depending  on  the  direction  of  the  buck¬ 
ling.  These  two  collapse  loads  correspond  to  points  E 
and  F  on  the  interaction  diagram  of  Fig.  14.2.  In 
general  these  two  loads  are  not  the  same  and  either  of 
them  can  be  the  smaller.  Since  almost  any  panel  can  be 
in  a  zero  lateral  load  condition,  both  values  must  be 
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calculated  and  the  lower  value  taken  as  the  collapse 
load. 


Combined  Loads 

In  ship  structures  the  stiffened  panels  are  usually  sub¬ 
jected  to  a  combination  of  bending  and  in-plane  loads. 
Also,  in  many  ship  panels  the  lateral  bending  can  be 
either  positive  or  negative.  For  example,  a  bottom 
panel  in  a  tanker  is  subjected  to  an  upward  or  a  down¬ 
ward  pressure,  depending  on  the  depth  of  cargo  oil  in 
that  space  relative  to  the  draft.  With  combined  loads 
there  is  an  interaction  among  the  three  basic  mech¬ 
anisms  outlined  above.  This  interaction  is  rather  com¬ 
plex  and  will  be  discussed  fully  in  the  next  three  sec¬ 
tions,  But  at  this  point  it  will  be  helpful  to  present  at 
least  a  qualitative  overall  view  of  the  simation,  and  this 
will  now'  be  done  with  the  aid  of  Fig.  14.2. 

The  first  and  most  basic  point  is  that  with  a  lateral 
load  the  plate-stiffener  "column”  becomes  a  "beam 
column,”  and  the  basic  column  ultimate  strength  the¬ 
ory  of  Section  11.2  is  replaced  by  the  beam  column 
theory  of  Section  11.3.  In  other  words,  collapse  still 
occurs  as  the  result  of  “failure  of  a  flange,”  but  .the 
bending  moment  and  deflection  caused  by  the  lateral 
load  must  be  taken  into  account.  Because  of  the  bend¬ 
ing,  either  flange  could  be  the  failure  flange,  de¬ 
pending  on  the  sign  of  the  bending  moment,  and  the 
failure  could  be  either  tensile  or  compressive.  Thus,  in 
principle,  there  are  four  modes  of  collapse;  in  practice, 
one  of  the  four — tensile  failure  of  the  plating — never 
occurs  because  the  neutral  axis  is  so  close  to  the  plat¬ 
ing.  Among  the  other  three  modes,  the  precise  namre 
of  the  collapse  mechanism  differs,  depending  on 
which  flange  is  the  compression  flange,  that  is,  de¬ 
pending  on  the  sign  of  the  bending  moment. 

For  greater  generality  we  shall  allow  for  the  possi¬ 
bility  of  the  stiffeners  having  a  different  value  of  yield 
stress  from  the  plating.  For  this  purppose  we  shall  use 
the  symbols  cry^  and  ayp.  We  will  also  need  to  refer  to 
the  mean  value  of  yield  stress  for  the  entire  section, 
which  is 


cry  = 


As 


O'Yvbt 

bt 


(14.1.1) 


Since  the  in-plane  load  usually  arises  from  hull 
girder  bending,  it  is  usually  more  fundamental  and 
more  variable  than  the  lateral  load.  Hence  it  is  custom¬ 
ary  to  define  "ultimate  strength”  in  terms  of  (Tq  and  to 
treat  the  value  of  the  lateral  load  as  an  independent 


parameter.  That  is,  the  ultimate  strength  of  a  panel  is 
the  value  of  a  a  that  causes  collapse,  in  the  presence  of 
a  specified  value  of  lateral  load.  In  an  interactive  col¬ 
lapse  diagram  such  as  Fig.  14.2  the  ultimate  strength 
is  the  horizontal  distance  to  the  curve  at  the  height 
corresponding  to  the  lateral  load.  Obviously  this  con¬ 
vention  is  not  appropriate  for  panels  in  which  the  lat¬ 
eral  load  is  always  the  dominant  load.  For  such  panels 
the  ultimate  strength  is  defined  in  terms  of  the  lateral 
load  and  can  be  calculated  by  the  plastic  hinge  theory 
of  Section  16.1*  But  for  the  general  case  of  a  panel 
subjected  to  combined  loads  the  ultimate  strength  the¬ 
ory’  must  be  sufficiently  general  to  handle  any  combi¬ 
nation  of  load.  The  theory  presented  in  this  chapter 
meets  this  requirement. 

Compression  Plus  Negative  Bending  (Fig.  14.3) 

Mode  I:  Compression  Failure  of  the  Stiffener.  With 
this  combination  of  loads  the  stiffener  flange  is  the 
compression  flange;  therefore  collapse  occurs  as  the 
result  of  compressive  failure  of  the  stiffener  flange, 
either  by  compressive  yield  or  by  buckling.  Let  us  first 
consider  compressive  yield.  When  there  is  any  appre¬ 
ciable  amount  of  axial  compressive  stress  it  di¬ 
rectly  increases  the  compressive  bending  stress  in  the 
stiffener  and  decreases  the  tensile  bending  stress  in  the 
plating.  Therefore  compressive  yielding  in  the  stiff¬ 
ener  commences  sooner  and  progresses  more  rapidly, 
whereas  plate  yielding  is  delayed.  The  result  is  that  the 
section  can  no  longer  achieve  a  plastic  hinge  condi¬ 
tion,  as  it  could  when  the  load  was  purely  lateral  (point 
B  in  Fig.  14.2).  Instead,  the  stiffener  reaches  its  limit 
of  stress  absorption  soon  after  the  yielding  has  pene¬ 
trated  through  the  stiffener  flange,  and  when  the  stiff¬ 
ener  is  thus  "neutralized”  the  effective  moment  of  in¬ 
ertia  of  the  section  becomes  very  small,  since  only  the 
plating  is  effective.  The  result  is  that  once  the  stiffener 
flange  is  fully  yielded  collapse  occurs  soon  afterward. 
Moreover,  this  mode  of  collapse  is  generally  quite 
sudden.  In  most  cases  the  precise  cause  of  the  failure 
is  the  formation  of  a  local  plastic  mechanism  in  the 
flange  such  as  that  shown  in  Fig.  14.4.  It  is  quite  local 
because  the  flange  becomes  fuUy  plastic  at  the  point 
where  the  bending  moment  is  a  maximum.  This  mech¬ 
anism  involves  large  local  sideways  deflection  and 
rotation,  permitting  the  stiffener  to  shorten  and  thereby 


*A  cross- stiffened  or  gross  panel  subjected  to  lateral  loads  is  termed 
a  grillage.  This  type  of  structure  can  be  designed  in  isolation, 
separate  from  the  hull  module,  and  is  therefore  not  particular  to  ship 
structural  design. 
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because  it  causes  tension  in  the  pbtmg  toward  the 

plating- 


Figure  14.3  In-plane  compression  and  negative  bending. 


escape  from  carrying  the  axial  compressive  load.  The 
local  sideways  deflection  of  the  flange  is  combined 
with  rotation  about  the  line  of  attachment  and  hence 
this  type  of  failure  is  sometimes  taken  to  be  tripping, 
but  \his  is  not  correct  since  the  cause  of  the  sideways 
deflection  is  a  local  plastic  mechanism  rather  than 
overall  flange  buckling. 

As  noted  earlier,  when  the  loading  is  purely  lateral 
the  plate-stiffener  combination  can  reach  a  plastic 
hinge  condition  before  collapse  occurs.  Hence  for  load 
combinations  which  involve  only  a  small  amount  of 
in-plane  compression  the  collapse  curve  can  depart 
from  the  "‘stiffener  flange  yield”  curve  and  eventually, 
at  zero  axial  load,  it  can  reach  the  plastic  hinge  col¬ 
lapse  point.  However,  the  curve  in  this  region  is  very 
steep,  such  that  a  small  underestimate  of  the  axial  load 
would  cause  a  serious  overestimate  of  the  collapse 
value  of  the  lateral  load.  For  this  reason  the  plastic 
hinge  mechanism  of  coUapse  should  only  be  used 
when  it  is  certain  that  in-plane  compression  cannot 
occur  (e.g.,  platforms,  ramps,  machinery  flats,  and 
other  laterally  loaded  stiffened  panels).  The  calcu¬ 
lation  of  ultimate  strength  under  this  type  of  loading  is 
dealt  with  in  Chapter  16.  Therefore,  for  stiffened  pan- 


Flgure  14.4  Local  plastic  collapse  mechanism  in  a  stiffener 
flange. 


els  which  are  subject  to  in-plane  compression  it  is  best 
to  use  the  “stiffener  flange  yield”  mechanism  for  all  of 
the  Mode  I  region,  even  thougli  for  very  small  in-plane 
loads  this  approach  is  conservative.  Hence  the  Mode  I 
collapse  curve  begins  at  point  A  in  Fig.  14.2. 

Alternatively,  the  stiffener  may  fail  by  tripping.  If 
the  amount  of  bending  is  small  the  elastic  tripping 
stress  (7a,r  is  given  by  the  formulas  of  Section  13. 1.  If 
it  is  not  small  then  the  stiffener  may  undergo 
flexural-torsional  buckling,  and  providing  the  stiffener 
remains  elastic  the  value  of  cTas  niay  be  calculated  by 
elastic  nonlinear  frame  analysis,  which  is  presented  in 
Section  15.5.  The  analysis  is  simple  and  rapid  because 
the  buckling  is  of  bifurcation  type  and  requires  only  an 
eigenvalue  analysis.  However,  as  mentioned  earlier, 
for  most  ship  panels,  buckling  is  not  purely  elastic  and 
the  simple  bifurcation/eigenvalue  approach  cannot 
provide  the  true  collapse  value,  cTa.u,  either  with  or 
without  Mo.  To  obtain  an  accurate  value  of  cr^.^  would 
require  the  fully  nonlinear  finite  element  method  (in¬ 
cremental  and  iterative),  which  is  discussed  in  Section 
15.6,  but  this  involves  too  much  computation  for  such 
a  local  and  frequently  occurring  calculation  (every 
panel)  especially  in  a  design  procedure.  Instead,  there 
is  a  simple  alternative  suggested  by  Smith  [1].  Since 
the  maximum  compressive  stress  in  the  flange  must  be 
kept  below  (Tys  in  order  to  avoid  a  local  plastic  mech¬ 
anism,  stiffener  buckling  can  be  avoided  simply  by 
requiring  that  the  value  of  cTas  obtained  from  the  elas¬ 
tic  theory  (Section  13.1  or  15.5,  as  appropriate)  must 
be  well  above  ctys. 

Since  collapse  occurs  when  the  stiffener  flange 
stress  reaches  yield  (or  earlier  if  <Ta.r  ^  <^rs)  Ai®  beam- 
column  remains  fully  elastic  prior  to  collapse.  There¬ 
fore  the  combination  of  compression  and  negative 
bending  is  relatively  straightforward  to  analyze,  re¬ 
quiring  only  the  ordinary  beam  column  theory  of  Sec¬ 
tion  11.3. 
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Compression  Plus  Positive  Bending  (Fig.  14.5) 

With  this  combination  of  loads  the  plating  is  the  com¬ 
pression  flange  and  there  are  two  distinct  modes  of 
collapse,  depending  on  the  relative  magnitude  of  the 
bending  and  the  in-plane  loads: 

Mode  II:  Compression  Failure  of  the  Plating.  With 
small  and  moderate  lateral  loads,  collapse  occurs  when 
the  stress  in  the  plating  reaches  the  failure  value  for  the 
‘^special  material’’  which  is  used  to  represent  the  plat¬ 
ing.  As  will  be  shown,  the  use  of  a  secant  modulus  to 
represent  the  plate’s  stress-strain  relationship  allows 
the  beam-column  to  be  treated  as  if  it  remained  en¬ 
tirely  clastic  prior  to  collapse. 

Mode  III:  Combined  Failure  of  Stiffener  and  Plat¬ 
ing.  With  large  positive  lateral  loads  the  bending 
causes  large  tensile  stresses  in  the  stiffener.  These  are 
reduced  by  the  in-plane  applied  compressive  stress, 
but  for  larger  and  larger  lateral  loads  the  stiffener  un¬ 
dergoes  more  and  more  tensile  yielding,  and  collapse, 
when  it  finally  occurs,  is  due  to  a  combination  of 
compressive  failure  of  the  plating  material  and  tensile 
yielding  of  the  stiffener. 

In  the  next  section  we  investigate  each  of  these 
three  modes  of  collapse. 


14.2  MODES  OF  COLLAPSE 


Mode  I  Collapse:  Compression  Failure  of  the 
Stiffener 

The  combination  of  in-plane  compression  and  nega¬ 
tive  bending  is  illustrated  in  Fig.  14.3.  With  this  com¬ 


bination,  collapse  occurs  due  to  compression  failure  of 
the  stiffener  flange.  That  is.  collapse  occurs  when  the 
total  stress  (cr^  plus  the  magnified  bending  stress) 
throughout  the  thickness  of  the  flange  has  reached  the 
“failure  value”  ct/t.  The  failure  value  is  either  the  yield 
stress  ays  or  the  elastic  tripping  stress  aaj,  whichever 
is  less.  Stated  mathematically 

CT/r  =  MIN{cry,,  aaj}  (14.2.1) 


where  is  obtained  from  the  method  of  Section  13.1 
or  15.5. 

Strictly  speaking,  as  the  yield  zone  penetrates  from 
the  outer  surface  to  the  midthickness  of  the  flange 
there  is  some  nonlinearity,  but  since  the  half-thickness 
is  small  relative  to  the  web  depth  the  effect  is  negli¬ 
gible.  Therefore,  the  failure  process  can  be  regarded  as 
entirely  elastic  and  can  be  completely  described  by  the 
beam-column  theory  of  Section  11.3.  Accordingly, 
the  total  stress  in  the  midthickness  of  the  flange  is 


Of  =  Oc 


M^f  ,  <7^A(5o  +  A)y;r 
I  I  ^ 


(14.2.2) 


in  which  A  and  I  =  cross-sectional  area  and  moment 
of  inertia  of  the  beam-column. 
(Note:  Since  the  plating  is  in  ten¬ 
sion  the  plate  flange  is  assumed 
to  be  of  full  width  b.  Therefore 
A  =  Aj  +  bt.) 

Mq  and  So  =  maximum  bending  moment  and 
maximum  deflection  due  to  the 
lateral  load  acting  alone.  It  is  as¬ 
sumed  that  these  maximum  val¬ 
ues  occur  at  the  midlength  of  the 
beam-column. 


Deflection  6^  :  positive  because 
it  is  tov^ard  the  stiffener. 


T 


(posrti>^) 


Bending  mornent  :  positive  because  it 
oduses  compression  in  the  plating 


Figure  14.5  In-plane  compression  and  positive  bending. 
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A  =  initial  eccentricity  of  the  beam- 
column:  a  typical  maximum  per¬ 
missible  value  for  welded  panels 
is  a/750,  where  a  is  the  length  of 
the  panel. 

yf  =  distance  from  the  centroidal  axis 
of  the  cross  section  to  the  mid¬ 
thickness  of  the  stiffener  flange. 
<p  —  magnification  factor  due  to  the 
axial  compressive  stress  cr;,: 


where 


R  = 


f  /? 

I' Vt  ^ 


(14.2.5) 


Therefore,  if  Of-  =  (Tyj  the  equacion  for  the  interactive 
collapse  curve  in  the  Mode  I  collapse  region  of  Fig. 
14.2  is 


(14.2.6a) 


(Te 


with  R  as  given  previously.  The  influence  of  Mq  is 
contained  in  the  bending  parameter  ji.  U  ctf  -  cTaj 
then  the  equation  becomes 


where  cte  is  the  Euler  buckling  stress  for  the  beam- 
column: 


CTe  = 


tt-EI 

Aa“ 


The  collapse  or  “ultimate'’  value  of  applied  stress 
(Ta^u  is  the  value  of  cr^  at  which  CTf  -  <Tf,  and  so  we  now 
proceed  to  make  these  substitutions  and  solve  for 
As  in  Section  11.3,  we  make  use  of  the  following 
nondimensional  parameters: 


where  p  = 


_  Mpyf 
I(Tf 


(14.2.3) 


In  terms  of  these  parameters  the  magnification  factor  is 

^  I  -  \-R 
Equation  14.2.2  becomes 


(1  -  R  -  tJL)il  -  X-R)  =  r)R  (14.2.4) 


and  the  solution  for  the  strength  ratio  R  is 


= /?—  {\A.2.6b) 

\(TyJi  <Ty, 

This  equation  is  based  on  the  following  sign  con¬ 
vention,  which  is  used  throughout  this  chapter. 

Stress:  Positive  if  compressive. 

Bending  moment:  Positive  when  the  plating  is  in 
compression. 

Lateral  deflection]  Positive  toward  the  stiffeners 
Eccentncity  (*•«•.  positive  upward  in  Fig. 

J  14,3).  Note:  For  a  simply 
supported  panel,  positive 
bending  causes  positive  lat¬ 
eral  deflection. 

Distance  from  neutral  axis:  Positive  toward  the 
plating. 

In  simply  supported  panels  that  arc  subjected  to  nega¬ 
tive  bending,  the  collapse  mode  is  always  Mode  I  tor 
any  significant  amount  of  negative  bending,  that  is, 
throughout  the  lower  half  of  the  interaction  diagram  of 
Fig.  14.2.  To  make  the  diagram  more  useful  the  loads 
are  nondimensionalized:  values  of  maximum  bending 
moment  are  normalized  by  dividing  by  iVfp,  and  values 
of  applied  in-plane  compressive  stress  are  divided  by 
ay.  TTius  the  point  corresponding  to  collapse  under  a 
pure  bending  load  is  the  point  —  1  on  the  bending 
moment  axis,  providing  that  buckling  does  not  precede , 
flange  yield.  As  the  amount  of  in-plane  compression 
increases,  the  amount  of  bending  which  is  required  to 
cause  collapse  is  lessened.  The  exact  shape  of  the 
interactive  collapse  curve  depends  on  the  scantlings 
and  geometry  of  the  panel.  If  the  failure  mechanism  is 
budding,  the  ultimate  strength  will  be  less  than  it 
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would  be  for  failure  by  yielding,  and  the  entire  curve 
will  lie  closer  to  the  ongin,  as  shown  by  the  dashed 
line  in  Fig.  14.2. 


Mode  n  Collapse:  Compression  Failure  of  the 
Plating 


We  next  consider  the  combination  of  in-plane  com¬ 
pression  (Ja  and  positive  bending  moment  Mo  as  shown 
in  Fig.  14.5.  If  the  bending  is  small  or  moderate  in 
magnitude,  collapse  occurs  due  to  compression  failure 
of  the  plate  flange.  If  the  plating  remained  pjerfectly 
elastic,  with  no  buckling  or  nonlinearity,  the  total  ap¬ 
plied  stress  acting  on  the  plating,  CTpo,  would  be  given 
by  the  usual  beam-column  formula 


O-a 


A 

<Vl(#+  5o)yp 
— 1 — * 


(14.2.7) 


As  shown  in  Section  12.6,  the  compressive  collapse  of 
welded  plating  is  a  complex  inelastic  process.  The 
relationship  between  load  and  end  shortening,  that  is, 
between  the  applied  stress  cr^  and  the  average  strain 
Cfl,  is  described  by  the  curves  of  Fig.  12.24.  For  our 
present  application  the  applied  stress  in  the  plating  is 
(Tpa  and  so  this  symbol  replaces  cFo  of  Fig.  12.24.  As 
an  illustration,  and  for  later  reference,  a  typical  curve 
of  (jpa  versus  Ea  from  Fig.  12.24  is  reproduced  in  Fig. 
14.6.  The  earlier  figure  showed  that  for  most  plates 
(j3  >  1)  the  relationship  between  apa  and  Eq  becomes 
nonlinear  well  before  collapse.  Both  figures  also  illus¬ 
trate  what  we  mean  by  “plate  failure,”  namely,  the 
point  in  the  versus  Ea  curve  at  which  the  plating  has 
lost  most  of  its  in-plane  stiffness.  In  geometric  terms 
it  is  the  point  at  which  the  tangent  of  the  curve  has 
become  significantly  less  than  the  original  value,  £. 
Strictly  speaking,  the  failure  point  is  the  peak  of  the 
curve.  However,  for  plates  of  low  and  intermediate 
slenderness  there  is  a  discernible  knee  in  the  curve 
(unless  the  initial  distortion  5p  is  very  large)  and  this  is 
the  simplest  and  safest  choice  of  the  failure  point.  For 
very  slender  plates  the  failure  point  is  either  the  peak 
of  the  curve  or  the  point  at  which  the  slope  has  de¬ 
creased  to  some  specified  value  (e.g.,  ^E). 

Since  the  plating  is  actually  being  subjected  to  an 
imposed  uniform  strain,  collapse  is,  strictly  speaking, 
measured  by  the  value  of  average  strain  corresponding 
to  this  loss  of  stiffness:  Eo  ==  Cuii,  as  shown  in  Fig. 
14.6.  However,  it  is  more  convenient  to  deal  in  terms 
of  stress.  For  later  use  we  define  ap  as  the  equivalent 
elastic  failure  stress,  that  is,  the  value  of  cr^  which 


Figure  14.6  Typical  curve  of  load  versus  end  shortening  for 
weided  steel  plating. 


would  correspond  to  Cuti  if  the  plate  maintained  a  linear 
relationship  between  apa  and  £a. 

For  most  plates  the  slope  of  the  curve  begins  to 
decrease  well  before  the  point  of  “plate  failure.” 
Therefore  the  true  value  of  applied  stress  which 
corresponds  to  plate  failure,  is  significantly  less  than 
o>.  This  also  means  that  during  the  load  sequence  that 
leads  to  plate  failure,  the  average  plate  stiffness  is 
significantly  less  than  the  elastic  or  material  stiffness 
£.  The  average  stiffness  over  the  entire  range  of  load¬ 
ing  is  given  by  the  secant  modulus  =  TE.  which  is 
the  slope  of  the  line  joining  the  origin  to  the  point  of 
plate  failure,  as  shown  in  Fig.  14.6.  For  our  purpose 
this  straight  line  can  be  used  in  place  of  the  actual 
curve  because  we  simply  want  to  be  able  to  predict 
plate  failure,  and  we  are  not  interested  in  the  detailed 
plate  response  prior  to  that.  With  this  approach  we  are, 
in  effect,  representing  the  plating  as  if  it  were  made  of 
a  different  material  which  has  a  smaller  value  of 
Young’s  modulus,  and  for  which  the  curve  of  apa  ver¬ 
sus  Ea  is  perfectly  linear,  right  up  to  failure.  As  we 
shall  see,  the  latter  property  is  of  very  great  advantage. 

The  mathematical  expression  for  the  secant  modu¬ 
lus  was  given  in  (12.6.4)  which  for  convenience  is 
reproduced  here 


r  = 


(14,2.8) 
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where 


Also,  it  was  noted  in  Section  12.6  that  the  ultimate 
strength  curve  has  the  same  shape  as  the  curve  for  T 
and  can  be  represented  sufficiently  accurately  by  sim¬ 
ply  subtracting  0.1  from  T.  That  is 


^  =  r-o.i 

a-Yp 

=  0.25(1.6  +  ^ 


with  ^  again  given  by 


(14.2.9) 


f  = 


-h 


2.75 

/3^ 


We  now  consider  the  composite  beam  column  com¬ 
prised  of  a  stiffener  and  a  plate  flange  of  width  b,  as 
shown  in  Fig.  14.7.  The  total  cross-sectional  area  is 
A  =  A,  +  bt.  The  loads  are  an  applied  compressive 
stress  (Ta  and  a  lateral  load  which  causes  a  positive, 
bending  moment  of  maximum  value  Mq,  The  plate 
flange  has  a  different  elastic  modulus  than  the  stiffener 
and  hence,  before  we  can  apply  beam  theory,  we  must 


Figure  14.7  Use  of  transformed  section  to  represent  combined 
stiffener-and-place  as  an  equivalent  eiasuc  beam  column. 


transform  this  composite  section  into  an  equivalent 
section  which  has  a  uniform  elastic  modulus,  as  was 
done  for  the  hull  girder  in  Section  3.6.  We  choose  the 
stiffener  as  the  reference  material,  and  hence  the  cross- 
sectional  area  of  the  plating  is  multiplied  by  the  trans¬ 
formation  factor  r.  In  the  present  application,  to  avoid 
any  confusion  about  plate  thickness  we  apply  T  to  the 
flange  width  b  instead  of  to  the  plate  thickness.  There¬ 
fore,  in  the  transformed  section  the  plate  flange  has  a 
transformed  width  btr  =  Tb.  The  total  area  is  now 
A,r  =  A^  +  bjy  the  neutral  axis  is  at  the  centroid  of 
the  transformed  section,  and  the  effective  moment  of 
inertia  of  the  beam-column  is  which  is  the  second 
moment  of  the  transformed  section.  The  strain  distri¬ 
bution  is  linear  and  continuous  but  the  stress  distribu¬ 
tion  is  discontinuous  because  in  the  plating  the  stress 
is  cTpa  =  Tcr%y  where  0%  is  the  stress  obtained  by 
applying  ordinary  beam  theory  to  the  transformed  sec¬ 
tion  (see  Fig.  14.7).  It  should  be  noted  that  cTpa,  the 
actual  value  of  stress  in  the  plating,  acts  across  the 
entire  width  of  the  plate  flange;  the  transformed  width 
brr  is  merely  a  device  for  obtaining  the  neutral  axis 
position  and  moment  of  inertia  of  the  composite  sec¬ 
tion. 

The  axial  stress  in  the  transformed  section  is 


As  4-  bt 
+  bj 


(14.2.10) 


and  this  is  the  value  of  stress  at  the  neutral  axis  of  the 
section.  If  we  now  apply  the  ordinary  beam  column 
theory  as  in  (14.2.7)  to  the  transformed  section  we 
obtain  the  uncorrected  plate  stress  o*^ 


,  A^’oyp.o-  ,  (ra.gAo-(A  +  8<i)y  o.cr  I 

0%  =  a-a.0.  +  -7 —  -r - - - <i> 

ttr  ^tr 

(14.2.11) 


and  the  true  value  of  (Xpa  is  then 

o-pa  =  Ta%  (14.2.12) 

From  beam-column  theory  we  know  that  collapse 
of  the  beam-column  occurs  when  the  stress  in  the 
compression  flange  (here  (Xp^)  reaches  the  failure  value 
for  that  material  (cr^a,  as  shown  in  Fig.  14.6).  The 
ultimate  strength  of  the  beam-column  is  defined  as  the 
value  of  applied  axial  stress  in  the  beam-column  (o'a.n-) 
which  is  sufficient  to  cause  the  above  collapse  condi¬ 
tion.  That  is,  (cTa.Juii  is  defined  as  the  value  of 
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which  satisfies 


(Tpu  =  T(t% 


(14.2.13) 


where  cr^  is  as  given  previously.  Rather  than  multi- 
plving  every  term  of  cr^  by  T,  it  will  be  more  con¬ 
venient  to  bring  T  to  the  left-hand  side  and  define  an 
equivalent  plate  failure  stress  (jf  \ 


__ 


as  was  foreshadowed  earlier  (see  Fig.  14.6). 

Thus  far  we  have  considered  only  longitudinal 
stress  in  the  plating.  If  there  is  a  transverse  stress  (Toy 
(positive  compressive)  and  shear  stress  r,  the  ultimate 
strength  of  the  plating  is  less  than  the  uniaxial  value 
a>«.  If  these  other  two  stresses  are  not  large  (less  than, 
say,  0.3  crj  the  value  of  (Tpu  can  be  corrected  by 
applying  a  reduction  factor  to  The  reduction  factor 
to  account  for  r  is  given  by  (12.6.3) 


Tr-  a/  1  “  3 


The  interaction  between  cr^i  and  (Toy  is  discussed  in 
Section  12.6,  and  is  shown  to  depend  on  the  plate 
slenderness  /3  and  aspect  ratio  a/b.  The  interaction  is 
quite  complex,  being  actually  beneficial  in  the  case  of 
sturdy  plates  (/3  <  1)  and  becoming  linear,  and  detri- 
mentk,  in  the  case  of  slender  plates.  For  design  pur¬ 
poses  we  want  a  simple  relationship  which  covers  the 
worst  cases  without  being  unduly  pessimistic  in  other 
cases.  The  worst  case  is  the  linear  interaction  given  by 
(12.6.16)  corresponding  to  slender  plates.  If  we  adopt 
this  then  the  reduction  factor  to  be  applied  to  dpu  is 

^ay  1 

in  which  is  given  by  (12.6.12)  assuming  that 
a/h  >  1.  We  now  apply  these  two  reduction  factors 
to  (Tpu,  and  for  simplicity  of  notation  we  do  so  by 
redefining  (Xf,  the  equivalent  failure  stress  for  the  plat¬ 
ing,  such  that  it  includes  these  factors.  We  also  take 
advantage  of  the  fact  that  crp„  can  be  represented  in 
terms  of  I.  as  in  (14.2.9):  =  (T  -  0.1)<ryp.  The 

definition  of  oy  becomes 


(Tf  =  - ~ - O-ypy  1  -  3 


and  (Tf  is  an  artificial  or  equivalent  failure  stress  for  the 
plating,  such  that  the  plate  fails  when  the  uncorrected 
plate  stress  ar^  reaches  ctf-  Substituting  for  ct*  from 
(14.2.11)  we  obtain  the  equation  governing  the  col¬ 
lapse  condition; 

Moyp^n  ,  ( ~F  ^)yp.fr  , 

(Tf  “  (trj^D,)ult  T  7  7 

Itr  ^77 

(14.2.15) 

In  parallel  with  the  procedure  adopted  for  Mode  I 
collapse,  we  define  the  following  nondimensional 
parameters 

„  (t7a.ir)ult 


Where  p.  =  yj- 


(A  +  5o)yn.r- 


(14.2.16) 


_  Mpyp.rr 
^  ItrO^F 


1  - 


In  terms  of  these  parameters,  (14.2.15)  reduces  to  the 
same  nondimensional  beam-column  equation  as  for 
Mode  I,  given  in  (14.2.4) 

(1  -  M  -  /?)(!  -  X^R)  =  vR 

and  the  solution  for  the  strength  ratio  R  or  (Ru,  to 
indicate  the  collapse  mode)  is 


“  2 


4 


(14.2.17) 


where 


1  4-  TJ 

=  1  -  At  + 


(14.2.14) 


This  equation  is  based  on  the  sign  convention  defined 
earlier  and  hence  5o,  and  A  are  all  positive  quan¬ 
tities. 

Induced  Load  EccENTRicm'  Due  to  Decreased 
Plate  Stiffness 

In  investigating  panel  collapse  of  a  stiffened  panel 
structure  we  usually  wish  to  find  and  analyze  the  first 
panel  to  collapse.  Hence  the  panels  which  are  immedi- 
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ately  adjacent  in  the  direction  of  the  load  should  be 
treated  as  fully  effective.  But  the  main  reason  for  treat¬ 
ing  them  as  such  is  that  this  is  a  more  severe  condition 
for  the  subject  panel.  The  reduced  stiffness  of  the 
plating  has  been  allowed  for  by  defining  a  transformed 
section  in  which  the  neutral  axis  is  displaced  slighdy 
further  from  the  plating.  By  definition  the  neutral  axis 
is  the  line  of  action,  within  the  cross  section  of  the 
beam  column,  of  the  total  axial  force 
P(P  =  AcTa  =  A^<ra.a-)  wlfich  is  transmitted  by  the 
cross  section.  But  the  line  of  action  of  the  externally 
applied  force  is  the  centroid  of  the  total  area  A.  which 
means  that  the  applied  load  has  a  slight  positive  eccen¬ 
tricity,  and  this  creates  an  extra  positive  moment  in  the 
section  and  therefore  decreases  the  ultimate  strength. 
The  eccentricity  is  equal  to  the  distance  moved  by  the 
neutral  axis,  which  is 


(14.2.18) 


where 


1  -t-  fip  (1  +  hp)A' 

Nonuniform  Applied  Compressive  Stress 


If  the  panel  forms  part  of  a  hull  girder  the  overall  depth 
of  the  hull  girder  is  sufficiently  large  that  the  hull 
girder  stress  cr^  can  be  regarded  as  constant  over  the 
height  of  the  panel.  If,  however,  the  overall  structure 
is  a  shallower  box  girder,  such  as  a  double  bottom  or 
a  barge  or  a  simple  box  girder  beam,  then  the  applied 
stress  cTa  has  a  trapezoidal  distribution  over  the  height 
of  the  panel  and  this  also  induces  some  additional 
eccentric  bending.  The  eccentricity  is  again  positive 
(providing  that  the  stiffeners  are  on  the  inside  of  the 
box  girder  flanges)  and  is  given  by 


J_ 

AH 


(14.2.21) 


where  h  is  the  distance  from  the  midplane  of  the  plat¬ 
ing  to  the  centroid  of  the  stiffener,  and  A,  is  the  sec- 
tionai  area  of  the  stiffener. 

The  effect  of  Ap  is  to  add  another  term  to  the  right- 
hand  side  of  (14.2.15).  Unlike  other  forms  of  eccen¬ 
tricity,  Ap  does  not  undergo  magnification  since  it  is 
the  final  value  of  the  neutral  axis  eccentricity,  at  the 
instant  of  failure.  Therefore  the  additional  term  is 

(CTa  .jr 

Irr 


where  A  and  /  refer  to  the  total  panel  cross  section  (not 
the  transformed  section)  and  H  is  the  distance  from  the 
centroidal  axis  of  the  total  panel  cross  section  to  the 
neutral  axis  of  the  overall  box  girder,  that  is,  to  the 
axis  of  zero  stress  in  the  box  girder,  taking  into 
account  any  axial  stress  which  may  be  present  in  it. 
Since  \h  is  an  initial  eccentricity  it  may  simply  be 
added  to  A,  and  there  is  no  need  for  a  new  term  as  there 
was  for  Ip, 

The  steps  in  calculating  the  Mode  II  ultimate 
strength  may  be  summarized  as  follows: 


Accordingly,  we  define  an  additional  nondimensional 
parameter  analogous  to  the  geometric  eccentricity 
parameter  tj 


(14.2.19) 


1.  For  the  given  lateral  load,  calculate  Mq  and  do- 

2.  From  (14.2.8)  calculate  T  and  then  calculate  ay 
from 


The  nondimensional  equation  governing  collapse  be¬ 
comes 


[1  -  M  -  A'/?)  =  rjR 


for  which  the  solution  is 


Ra 


1  ^ 

(1  4*  T7p)A- 


(14.2.20) 


with  cray,u  given  by  (12.6.12). 

3.  For  a  transformed  section  having  a  plate  flange  of 
width  Tb,  calculate  Arr,  yp,a-y  Itry  Ptry  and  Ap.  (Note: 
Eqs.  8.3.6  are  helpful  for  this  task. 

4.  Using  (14.2.16)  and  (14.2.19)  calculate  the  pa¬ 
rameters  A,  7],  7]p,  and  /x. 

5.  Calculate  the  strength  ratio  from  (14.2.20). 

6.  Calculate  the  collapse  load  (ra,u  from  (14.2.10) 


fri.u 
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which  is,  in  terms  of 


for  this  collapse  mode  is  convex  upward  and  that  it 
always  intersects  the  Mode  n  curve  at  some  distance 
above  point  G.  That  is,  the  Mode  U  collapse  analysis 
becomes  inaccurate  only  gradually,  in  proportion  to 
the  extent  to  which  Mo  exceeds  (Mo)c»  such  that  the 
collapse  involves  some  yielding  of  the  stiffener  flange. 
Therefore  we  choose  a  straight  line  between  points  D 
and  G  as  the  interaction  curve  for  Mode  III  collapse. 
This  choice  is  safe  since  the  actual  collapse  points 
alwavs  lie  slightly  above  this  line,  and  it  greatly 
simplifies  the  analysis  of  Mode  ffl  coUapse.  The  anal¬ 
ysis  only  requires  finding  point  G,  the  load  combina¬ 
tion  which  causes  simultaneous  compression  failure  of 
the  plating  and  tensile  yield  of  the  stiffener  flange.  The 
latter  occurs  when  the  total  flange  stress  reaches  the 
tensile  yield  stress,  which  in  the  sign  convention  of 
this  chapter  is  ^orys-  Therefore  the  condition  for  flange 
yield  is 


Mode  m  Collapse:  Combined  FaUure  of  Stiffener 
and  Plating 

If  the  positive  bending  moment  Mo  is  large  it  creates  a 
large  tensile  stress  in  the  stiffener.  This  is  partly  can¬ 
celed  by  the  applied  compressive  stress,  but  if  the 
bending  moment  is  sufficiently  large  there  will  be  ten¬ 
sile  yielding  in  the  stiffener  flange  as  well  as  com¬ 
pressive  failure  of  the  plating.  This  further  reduces  the 
effectiveness  of  the  section,  and  so  the  plate  failure 
analysis  of  the  previous  section  is  no  longer  sufficient. 
Tlie  point  on  the  Mode  II  curve  where  this  simulta¬ 
neous  failure  begins  to  occur  is  labeled  G  on  the  inter¬ 
action  diagram  of  Fig.  14.2.  This  simultaneous  failure 
occurs  in  all  combinations  of  collapse  loads  in  which 
the  bending  moment  exceeds  (Mo)g-  Therefore  in  this 
region  of  the  interaction  diagram  there  is  a  new  mode 
of  collapse— Mode  IQ— and  a  new  segment  of  the 
collapse  curve  extending  from  point  G  up  to  point  D, 
where  the  load  is  a  pure  bending  load.  For  this  loading, 
the  stiffener  can  achieve  full  yield  and  as  this  occurs 
the  effective  neutral  axis  moves  toward  the  plating. 
For  most  panel  proportions  it  eventually  enters  the 
plating,  putting  it  partly  in  tension.  This  stabilizes  the 
plating  such  that  it  can  absorb  almost  the  full  value  of 
compressive  yield  stress.  Hence  for  a  purely  lateral 
load  collapse  does  not  occur  until  the  panel  has 
reached  a  plastic  hinge  condition. 

A  precise  calculation  of  the  collapse  loads  for  Mode 
HI  collapse  would  be  extremely  difficult  because  of  the 
complicated  interaction  between  the  simultaneous 
stiffener  tensile  yielding  and  plate  compressive  fail¬ 
ure.  Experiments  [2]  show  that  the  interaction  curve 


Moyfrr  (trfl,rr)uJt^rr(5o  A)V/'.n'  • 

-cry,  =  (cra.rr)uli  +  ,  “  T  • 

Itr 

,  Apyy.fr 


(14.2.23) 


This  equation  corresponds  to  the  curve  GH  in  Fig. 
14.2.  To  facilitate  the  solution  for  the  value  of  (cTa.Juii 
we  again  define  a  set  of  nondimensional  parameters. 
These  differ  from  those  of  Mode  n  in  that  the  flange 
distance  yy  ^.  replaces  Vp.fn  which  means  that  both  ry  and 
Tjp  will  have  negative  values.  Also,  since  the  failure 
stress  is  now  negative  (ay  =  ^cfys)  the  strength  ratio 
R  will  likewise  have  a  negative  value,  and  this  means 
that  the  expression  for  the  magnification  factor  must  be 
altered  to 


cf>  = 


1 

1  + 


In  order  to  distinguish  the  nondimensional  parame¬ 
ters  which  correspond  to  the  curve  GH  from  those  of 
Mode  n,  we  shall  use  a  subscript  GH.  Accordingly, 
the  parameters  are 


Rgh 


(^a.rr)ult 


Acw  ""  \  r- 


(A  -r  So)yy,rr 


Vgh 
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(14.2.24)  The  result  will  be  positive  since  Rqh  is  negative.  R^h 
is  a  function  of  (iWo)^  since 


= 


^fr(  ^Ys) 


I  ^  X^Rc 


The  condition  for  flange  yield  (i.e.,  the  equation  of 
curve  GH)  is  then 


[l  —  Rcni^  Vp.gh)](^  ^chRgh) 


and  solving  for  Rch  gives* 
p  _  ^GH  ^  Ugh  .  i 

2  ’  V  4  ^  (1  -  V,,CH)yGH 


where 


r  _  ^  f^H  _  1  Vp.GH  Vgh 
1  Vp,GH  (1  Vp,Gh)^GH 


The  other  condition  defining  point  G  is  the  Mode  II . 
condition.  Therefore,  from  ( 14.2.20)  the  second  equa¬ 
tion  for  determining  point  G  is 


1  ~  Mu 

(1  '1'  7?/7.u)Au 


where 


(14.2.26) 


r  r=  ^  1  Vp.a  +  Vu 

°  1  +  Vp.a  (1  Vp.iD^h 

in  which  the  nondimensional  parameters  are  those  of 
(14.2.16)  and  (14.2.19).  From  these  expressions  for 
Rch  and  Rq  we  can  obtain  two  equations  for  (cr^.u)G  as 
a  function  of  (Mo)q,  For  the  stiffener  tensile  yield  con¬ 
dition  the  equation  for  {(Ta,u)G  in  terms  of  Rch  is  similar 
to  (14.2.22),  namely 


__  (i^())cyf.tr 
-  ~  ’ 

IA-o-ys) 

The  second  equation  for  (cTa.Jc  and  (A/o)g  comes  from 
the  plate  failure  condition,  for  which  (14.2.22)  gives 


(14.2.28) 

^Y  )  G  CTjrV  A  / 


=  r]cHRGH  which  R^  is  a  function  of  (iWo)c,  since 

(.^O^Gyp.cr 


The  required  collapse  loads  (o-a^Jc  and  (Mo)g  are 
obtained  by  solving  (14.2.27)  and  (14.2.28)  simulta¬ 
neously,  using  an  appropriate  numerical  technique.  By 
equating  the  right-hand  sides  of  these  equations  we 
obtain  a  single  implicit  equation  for  (iV/o)c* 


-cry, /a 


^Y/G  CTy 


RnCTp  —  ^RghO'ys 


(14.2.29) 


The  solution  must  be  performed  iteratively,  and  the 
best  approach  is  to  begin  with  a  trial  value  of  (Mo)g. 
calculate  the  corresponding  value  of  (5o)c  (which  de¬ 
pends  on  the  type  of  lateral  load),  then  evaluate  /?a  and 
Rgh^  aiid  hence  the  left  and  right  sides  of  ( 14.2.29).  In 
general  these  will  not  agree  and  it  will  be  necessary  to 
choose  another  value  of  (A/qIg  and  repeat  until  die 
agreement  is  satisfactory.  Some  numerical  technique 
(Newton’s  method.  Secant  method,  etc.)  may  be  used 
to  obtain  each  new  value  of  (Mo)c’ 

Since  the  other  end-point  of  the  Mode  HI  collapse 
line  is  (cr^.u/^K  ==  0*^  X^olMp  =  I)  the  equation  for  the 
ultimate  strength  of  a  panel  in  the  Mode  III  collapse 
regime  is 


^a.u  I  _  ~  Mq  ( (0"a.u)G 

(Ty  )  m  ^P  ”  (^o)g  \  ^Y 


(14.2.30) 


Rgh  (14.2.27)  in  which  (o-a^Jc  is  obtained  from  either  of  two  expres- 


•We  here  indicate  both  roots  of  the  quadratic  because  the  solution 
for  point  G  can  involve  either  root,  depending  on  Mo  and  the  panel 
proportions. 


(t^a,ii)G 


o-A  — 
\A 


(14.2.31) 
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or 


Effect  of  Initial  Eccentricity 

If  there  is  no  lateral  load  the  initial  eccentricity  A 
becomes  the  crucial  factor  in  determining  which  mode 
of  collapse  occurs.  The  load  eccentricity  Ap  which  is 
induced  by  the  shift  of  the  neutral  axis  at  plate  failure 
is  always  posidve;  that  is.  it  is  always  such  as  to 
encourage  a  Mode  11  collapse .  Therefore  if  the  initial 
geometric  eccentricity  A  is  positive,  then  the  theory 
necessarily  predicts  a  Mode  El  collapse  because  both  of 
the  eccentricities,  A  and  Ap,  are  positive.  Of  course 
there  are  other  possible  factors,  such  as  stiffener  im¬ 
perfections  and  eccentricity  of  applied  load,  which  can 
override  this  tendency  and  cause  a  Mode  I  collapse. 
But  it  is  a  fact  that  panels  which  have  a  positive  initial 
eccentricity  and  are  subjected  to  a  purely  in-plane  load 
usually  collapse  in  Mode  n.  The  converse  is  not  true; 
that  is,  a  panel  with  a  purely  axial  load  and  with  a 
small  or  moderate  negative  initial  eccentricity  may 
collapse  in  either  Mode  I  or  Mode  II.  This  is  because 


the  induced  load  eccentricity  Ap  is  always  positive,  and 
hence  it  may  eclipse  a  negative  value  of  A. 

Also,  a  negative  initial  eccentricity  can  give  an 
apparent  strengthening  against  a  Mode  n  collapse.  A 
negative  deflection  relieves  the  plate  of  some  com¬ 
pression,  and  therefore  a  small  or  moderate  negative 
initial  eccentricity  (not  large  enough  to  cause  a  Mode 
I  collapse)  may  delay  plate  failure  and  hence  give  a 
small  but  measurable  increase  in  the  Mode  n  collapse 
strength. 

This  strengthening  effect  is  not  restricted  to  the  case 
of  pure  in-plane  loading.  For  any  lateral  load  it  will 
occur  whenever  A  is  of  opposite  sign  to  5o,  whether 
positive  or  negative.  However,  the  effect  is  significant 
only  for  small  lateral  loads,  when  5^  is  of  the  same 
order  of  magnitude  as  A.  In  this  case  tj  is  very  small 
since  it  is  proportional  to  the  sum  of  these  two  eccen¬ 
tricities.  The  effect  which  this  has  on  the  collapse 
curves  is  shown  in  Fig.  14.8.  exaggerated  for  clarity'. 
Some  examples  of  this  effect  in  actual  panels  are 
presented  in  Section  14.4. 

14.3  ULTIMATE  STRAIN  OF  STIFFENED 
PANELS 

As  mentioned  earlier,  in  overall  structural  collapse  due 
to  hull  girder  bending  the  loading  condition  on  each 
plate  panel  is  acmally  a  uniform  applied  strain,  e.. 
rather  than  a  uniform  applied  stress,  and  it  is  more 
correct  to  express  plate  collapse  in  terms  of  the  value 
of  strain  in  the  plate  when  it  collapses,  £a.u-  The  same 
is  true  for  the  stiffener-plate  combination  and  so  it  is 
necessary  to  have  an  expression  for  Ec.u  fot  this  type  of 
member.*  The  expression  depends  on  the  degree  of 
yielding  in  the  cross  section,  which  in  turn  depends  on 
the  collapse  mode.  For  Mode  I  the  cross  section  is  still 
elastic  and  so  the  ultimate  strain  is 


with  o-a.u  given  by  (14.2.6). 

For  Mode  n  the  stiffener  is  again  elastic  since  the 
upper  end  of  the  Mode  El  range  is  stiffener  tensile 
yield.  The  plating  is  inelastic  but  it  has  been  replaced 
by  a  pseudoelastic  plate,  having  an  effective  elastic 
modulus  Thus  the  stiffener  and  the  plate  are  like 
two  parallel  elastic  rods  for  which  the  axial  stiffnesses 
are  EAJL  and  E,Ap/L.  Since  they  are  in  parallel  and 
must  have  the  same  axial  shortening,  the  total  stiffness 
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is  the  sum  of  the  individual  stiffnesses.  Expressing  the 
total  stiffness  in  terms  of  the  total  ar^  A  =  Ap  A, 
and  a  mean  value  of  secant  modulus  E,  gives 

EA  _  EA,  EAp 

and  by  noting  that  A„.  =  A,  +  TAp  we  have 

E,  =  —E  (14.3.2) 

A 


=3 

=^3 


=0 


The  ultimate  strain  is  then 

A  (TaM 
Arr  E 

in  which  cra,u  comes  from  the  procedure  given  at  the 
end  of  the  discussion  of  Mode  II  collapse,  during 
which  Arr  is  also  evaluated. 

For  Mode  in  the  degree  of  yielding  in  the  stiffener 
is  not  known  precisely  and  the  safest  and  easiest  course 
is  to  regard  the  stiffener  as  being  fuUy  yielded,  with  an 
effective  tangent  modulus  of  zero.  The  mean  value  of 
secant  modulus  for  the  panel  is  then 


(a)  Panei  Geometry 


W  ,  W 


i  I 


(14.3.3) 


(14.3.4) 


(b)  Loading 

Figure  14.9  Experimental  setup  of  Monash  tests. 


and  the  ultimate  strain  is 


A  a-g.u 
Ap.<r  E 


(14.3.5) 


each  stiffener  at  the  a/3  positions,  giving  a  bending 
moment  Wo  =  Wa/3  which  is  constant  over  the  mid¬ 
dle  one-third  of  the  panel.  For  this  load  the  maximum 
deflection  is 
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= 


23H^a^ 

648£/ 


(14.4.1) 


Single  Span  Panels 

In  order  to  verify  the  beam  column  approach  presented 
in  the  previous  sections,  we  now  compare  the  collapse 
loads  predicted  by  this  method  with  the  measured  col¬ 
lapse  loads  of  a  series  of  fifteen  panels  which  were 
tested  to  destruction  under  various  combinations  of 
lateral  and  in-plane  load  by  Michelutti  and  Murray  at 
Monash  University  in  Melbourne  [2,3,4].  The  experi¬ 
mental  semp  is  shown  in  Fig.  14.9.  The  bending  mo¬ 
ment  was  provided  by  two  point  loads  W  applied  to 


Also  the  plastic  hinge  load  Wp  is  related  to  Mp  by 


Mp 


Wpa 

3 


(14.4.2) 


Thirteen  of  the  panels  were  of  the  same  (nominal) 
overall  dimensions,  namely: 

a  =  3.20  m,  B  =  2.44  m,  b  =  0.533  m, 

=  0.152  m,  tp  =  9.7  ±  0.3  mm, 

=  7.3  —  0.2  mm 
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One  panel  (panel  A4  in  the  temunology  of  Ref.  2)  was 
slightly  shorter  (a  =  3.05  m)  but  was  geometnc^ly 
similar  to  the  others.  Only  one  panel  (panel  A6)  had 
any  significant  difference  in  geometry;  panel  A6  was 
54%  longer,  having  a  =  4.93  m. 

The  principal  difference  among  the  panels  was  in 
regard  to  yield  stress.  Within  each  panel  the  stiffeners 
and  plating  had  different  yield  stresses,  and  the  rauo 
varied  from  panel  to  panel.  In  nine  panels  o-y,  exceeded 
ayp,  the  three  largest  ratios  being  1.25  (panel  Q)  and 
1.19  (panels  A5  and  N)  and  in  six  panels  o-y.  was  less 
than  rtyp,  the  two  smallest  ratios  being  0.86  (panel  L) 
and  0.91  (panel  A3).  Details  are  given  in  Table  14.1. 
including  the  value  of  initial  eccentricity-  A /a  and  the 
values  of  various  parameters.  The  table  gives  the  ex¬ 
perimental  coUapse  loads  and  the  corresponding  nor- 
miized  values  a^J^y  and  W  jWp.  This  is  foUowed 
by  the  theoretical  normalized  values  obtained  from  the 
theory  of  the  preceding  sections. 

Although  fourteen  of  the  panels  were  vutually  iden¬ 
tical  in  geometry,  the  variation  in  the  ratio  of  yield 
stress  means  that  each  panel  has  a  different  collapse 
curve.  For  each  panel  the  complete  curve  was  calcu¬ 
lated  using  the  particular  values  of  plate  and  stiffener 
thickness  and  of  initial  eccentricity  of  that  panel,,  in 
order  to  make  the  comparison  of  theory’  and  experi¬ 
ment  more  precise.  Some  examples  are  shown  in  Fig. 
14.10.  These  have  been  chosen  to  cover  the  full  range 


of  collapse  load  combinations  and  to  illustrate  the  va¬ 
riety  among  the  curves.  For  clarity  most  of  the  curves 
are  not  plotted  in  full  but  only  in  the  vicinity  of  the 
experimental  collapse  load  for  that  panel. 

Panel  Q  is  one  of  the  panels  with  a  large  value  of 
ayjayp  and  hence  its  Mode  ffl  curve  is  higher.  Panels 
A2,  H,  and  Z  are  typical  of  most  of  the  panels,  with 
a  yield  stress  ratio  close  to  unity.  Panels  H  and  A6 
were  subjected  to  pure  in-plane  compression  and  they 
collapsed  due  to  plate  failure  (Mode  ID.  The  behavior 
of  stiffened  panels  that  are  subjected  to  a  purely  axial 
load  is  an  important  aspect  of  stiffened  panel  response 
and  merits  some  discussion.  Therefore  we  shall  use 
these  two  panels  as  examples  to  illustrate  and  explain 
this  behavior. 

When  the  loading  is  pure  in-plane  compression  the 
initial  eccentricity  has  a  strong  influence  on  whether 
coUapse  occurs  in  Mode  I  or  Mode  H.  Also,  as  dis¬ 
cussed  in  the  previous  section,  a  negative  initial  eccen¬ 
tricity  can  give  an  apparent  strengthening  against  a 
Mode  n  collapse;  This  phenomenon  occurred  m  panel 
H  Unfortunately,  neither  the  direction  nor  the  mag¬ 
nitude  of  the  initial  eccentricity  were  recorded  for 
panel  H.  However,  it  appears  clear  from  the  test  re¬ 
sults  that  the  initial  eccentricity  was  negative  because 
the  measured  values  of  lateral  deflection  are  negative 
at  all  stages  of  the  loading,  until  just  before  collaps^ 
At  an  applied  load  of  (T^  =  232  MPa  {cTaJ  cry  - 


TABLE  14.1  COMPARISON  OF  THEORY  AND  EXPERIMENT 


L 

LI 

Y 

Z 

N 

J 

H 

A6' 

Q 

A5 

A4 

A3 

A2 

A1 

P 


Panel  Characieristics 


Experiment 


Theon*" 


Discrepancy 


Panel 


Cf'ii 


315 

396 

363 

367 
332 
370 

378 

368 
318 
324 
372 
370 
383 
378 
324 


381 

389 


Uys 

(Typ 


378 
370  1.01 

396 


383 
397 
385 
361 
338  0.91 


371 

369 


(Ty 


270  0.86  305 

379  0.96  391 

1.04  367 


368 
1.19  346 

1.03  374 


1.03  381 


1.04  372 

1.25  335 


1.19  338 

0.97  369 


363 

354  0.92  378 

0.98  376 

1.14  334 


a/ A 


“3150 

1600 

1580 

6300 

930 

“1200*’ 

not 
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530 

890 

970 

1330 

1150 
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1400 
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Ai.ii.c 


0.79  (I) 
0.94  (1) 
0.94  a) 
0.93  (I) 
0.97  (I) 
0.95  (I) 


0 

0 

62 

123 

184 

235 


0.80  (H)  253 


i.2r(n) 
0.75  (H) 
0.75  01) 
0.76  01) 
0.82  (GH) 
0.84  (GH) 
0.85  (GH) 
0.88  (GH) 


153 
172 
170 

154 
139 
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93 

0 


W'ult 

(To.y 

a> 

CTom 

O-Y 
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9c 

“39.0 

0 

“0.89 

0 

“0.72 

0.17 

19.1 

“39.6 

0 

-0.71 

0 

“0.71 

0 

0 

“29.8 

0.17 

“0.65 

0.15 

“0.58 

0.08 

11.9 

“14.6 

0.33 

-0.31 

0.31 

“0.30 

0.02 

4.4 

“9.0 

0.53 

“0.19 

0.51 

-0.18 

0.02 

3.5 

0 

0.63 

0 

0.64 

0 

0.01 

1.6 

0 

0.66 

0 

0.59 

0 

0.07 

10.6 

0 

0.41 

0 

0.39 

0 

0.02 

4.9 

11.1 

0,51 

0.24 

0.51 

0.24 

0 

0 

22.8 

0.50 

0.49 

0.46 

0.45 

0.06 

8.6 

36.4 

0.42 

0.65 

0.38 

0.59 

0.07 

9.6 

37.0 

0.38 

0.69 

0.33 

0.59 

0.12 

15.2 

38.4 

0.37 

0.71 

0.31 

0.60 

0.12 

15.1 

42.0 

0.25 

0.79 

0.24 

0.74 

0.05 

6. 1 

46.8 

0 

1-02 

0 

1.00 

0 
0 
1  j 

2.0 

,  r  /A  -*-7^0  for  Mode  I*  for  Modes  U  and  m). 

^All  theoretical  results  are  for  a/A  “  -750  (  ror  mooc  a  =  -/7/12OO 

Tn  this  panel  the  eccentricit>  was  in  the  axial  load,  and  was  equivalent  to  A  / 

‘Xxjng  panel  (fl  =  4.93  m):  therefore  A  is  large. 
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p 


panels. 


232  /  381  =  0.61)  the  plating  began  to  fail  and  the 
panel  inunediately  deflected  in  the  positive  direction 
(i.e.,  toward  the  stiffeners).  On  further  increase  of 
load,  the  positive  deflection  grew  rapidly  and  a  Mode 
H  collapse  occurred  soon  after,  at  =  0.66. 

With  a  positive  eccentricity  (e.g.,  the  standard  design 
value  of  +a/750)  the  beam  column  theory  for  Mode 
n,  as  outlined  previously,  gives  a  collapse  load  ratio  of 
0.59.  Therefore  it  appears  that  the  panel  had  a  negative 
initial  eccentricity  which  artificially  strengthened  it. 
Obviously,  one  must  ignore  such  strengthening  when 
designing,  but  in  comparing  theory  and  experiment  it 
is  important  to  model  the  actual  situation  as  closely  as 
possible.  Therefore  a  small  negative  initial  eccen¬ 
tricity  (—a/ 2000)  was  assumed  and  the  theoretical 
collapse  loads  were  calculated  for  Mode  I  and  Mode  U 


(with  a  small  negative  value  of  A,  both  equations  have 
solutions).  The  resulting  values  of  cTa.ul'^Y  were  0.77 
and  0.67  for  Modes  I  and  II,  as  shown  in  Fig.  14. 10. 
Thus  when  the  negative  initial  eccentricity  which  ap¬ 
pears  to  have  been  present  is  represented  by  a  typical 
value,  the  theory  correctly  predicts  that  Mode  II  is  the 
collapse  mode,  and  the  predicted  value  agrees  closely 
with  the  experimental  v^ue  of  0.66.  For  comparison. 
Fig.  14. 10  also  shows  the  collapse  curves  correspond¬ 
ing  to  the  standard  design  values  of  initial  eccentricity 
(-a/750  for  Mode  I  and  H-a/750  for  Mode  H).  These 
arc  the  curves  which  would  in  practice  be  taken  as  the 
predicted  capability  of  the  panel.  Since  they  are  design 
curves  they  must  always  give  a  predicted  collapse  load 
which  is  less  than  the  actual  coUapse  load  for  any  panel 
whose  initial  eccentricity  is  within  the  permissible 
range:  -a/750  <  A  <  a/750.  At  the  same  time  they 
should  not  be  excessively  conservative.  As  shown  in 
the  figure,  the  curves  for  panel  H  fulfill  both  of  these 
requirements. 

Panel  A6  is  the  long  panel  and  consequently  it  has 
a  much  larger  slenderness  parameter:  Aj  =  1.46  and 
An  =  1.21,  compared  to  typical  values  of  0.95  and 
0.75  for  the  other  panels.  This  causes  a  substantial 
decrease  in  ultimate  strength  in  both  modes.  In  fact, 
the  stiffener  tensile  yield  curve  for  panel  A6  (the  curve 
GH  in  Fig.  14. 10)  is  so  close  to  the  origin  that  it  does 
not  intersect  the  Mode  II  curve,  and  the  Mode  m 
collapse  curve  is  a  straight  line  from  the  point 
Mo/Mp  =  I  to  the  point  where  the  GH  curve  crosses 
the  horizontal  axis.  As  noted  earlier,  this  is  conser¬ 
vative  because  in  such  cases  the  collapse  load  combi¬ 
nations  would  foUow  the  Mode  II  collapse  curve  up¬ 
ward  for  some  distance  (i.e. ,  up  to  moderate  values  of 
Mo),  But  the  conservativism  is  not  excessive  and  it 
emphasizes  the  fact  that  such  a  panel  does  have  gener¬ 
ally  low  values  of  nondimensional  collapse  loads. 

Panel  A6  also  provides  another  Illustration  of  the 
crucial  influence  which  the  initial  eccentricity  has 
when  there  is  no  lateral  load.  For  this  panel,  when  the 
standard  design  value  of  —a/750  is  assumed  for  Mode 
I  and  'ra/750  for  Mode  H,  the  two  values  of  cTa.u/ 
are  0.376  and  0.405  respectively,  as  shown  in  the 
figure.  Thus  if  we  were  considering  this  panel  in  an 
objective  design  context,  in  which  both  positive  and 
negative  eccentricity  must  be  allowed  for,  the  panel 
could  be  said  to  be  weaker  in  Mode  I  than  in  Mode  II. 
But  in  this  specific  case  the  panel  had  a  relatively  large 
positive  initial  eccentricity  ('ra/530)  and  therefore 
Mode  I  collapse  could  not  and  did  not  occur.  For  a 
Mode  I  collapse  to  occur,  A  would  have  to  be  negative 
(and  sufficiendy  negative  to  overcome  the  induced 
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eccentricity  effect).  When  the  measured  value  of 
-rfl/530  is  used  for  the  Mode  11  calculation,  the  pre¬ 
dicted  collapse  strength  is  0.390,  which  agrees  well 
with  the  experimental  result:  the  panel  in  fact  under¬ 
went  a  Mode  n  collapse  at  an  applied  load  ratio  of 
0.41. 

Since  the  tests  covered  a  full  range  of  load  combina¬ 
tions,  including  purely  lateral  loads  and  purely  in¬ 
plane  loads,  the  iscrepancy  between  theory  and  ex¬ 
periment  is  best  measured  along  a  radial  line  from  the 
origin  as  shown  for  panels  A2,  A5,  N,  and  Z  in  Fig. 
14.10.  The  relative  error  is  then  obtained  by  dividing 
this  radial  discrepancy  by  the  length  along  the  radial 
line  from  the  origin  to  the  experimental  collapse  point. 
This  has  been  done  for  all  fifteen  cases  and  the  results 
are  given  in  Table  14.1.  Also,  to  give  a  clearer  overall 
comparison.  Fig.  14.11  contains  all  of  the  radial  dis¬ 
crepancies  plotted  in  their  correct  positions  relative  to 
a  single  typical  collapse  curve.  The  calculations  were 
performed  using  the  standard  design  values  of  initial 
eccentricity,  ±a/750.  The  results  indicate  that  the 
beam  column  equations  presented  in  the  previous  sec¬ 
tions  are  of  satisfactory  accuracy,  and  that  the  value  of 
±a/750  for  the  initial  eccentricity  is  suitable  for  de¬ 
sign.  In  all  panels  except  A6  (the  long  panel)  the 
design  value  of  eccentricity  exceeded  the  actual  value 
(see  Table  14.1).  The  design  value  of  ±a/750  is  a 
maximum  permissible  value;  panels  with  larger  eccen¬ 
tricity  would  either  be  forbidden  by  the  fabrication 
specifications  or  would  be  designed  as  curved  panels. 
For  all  panels  in  which  the  initial  eccentricity'  was 
within  these  permissible  limits  the  actual  value  of  ulti¬ 
mate  strength  was  never  less  than  the  design  value,  and 
the  design  value  was  never  excessively  conservative. 
Moreover,  even  in  panel  A6,  in  which  the  initial  ec¬ 
centricity'  was  slightly  over  the  permissible  value,  the 
collapse  load  was  still  slightly  above  the  design  value. 
These  results  demonstrate  that  the  beam  column  ap¬ 
proach  is  well  suited  for  the  design  of  stiffened  panels, 
since  it  gives  an  accurate  estimate  of  strength  for  all 
load  combinations  and  for  any  typical  amount  of  initial 
eccentricity. 

Multispan  Panels  Without  Lateral  Load 

To  demonstrate  the  validity'  of  the  method  for  multi¬ 
span  panels  we  next  compare  its  predicted  values  of 
(Ta,u  and  £a,u  with  the  measured  values  obtained  by 
Dowling  et  al.  [5]  in  a  test  of  a  3.9-  m  long  box  girder 
subjected  to  a  sagging  bending  moment  (model  no.  2 
in  the  terminology  of  Ref.  5).  The  cross  section  is 
shown  in  Fig.  14.12.  The  upper  flange  (the  “strength 
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deck”)  consisted  of  five  in-line  panels  787  mm  in 
length.  The  panels  were  typical  in  ail  respects  to  ship 
panels;  details  are  given  in  the  figure  and  in  Ta¬ 
ble  14.2.  From  Ref.  5,  in  the  bay  where  the  collapse 
occurred  the  conditions  before  loading  were: 

maximum  initial  deflection  of  stiffeners: 
upward  (towards  plating) 

A  =  —0.010  in.  =  —0.25  mm  (a/3100) 
downward  (towards  stiffener) 

A  =  0.055  in.  =  1.40  mm  (a/ 563) 
average  longimdinal  residual  stress: 

CTr  =  O.li  cry 

Since  there  is  no  lateral  load,  both  Mode  I  and 
Mode  n  collapse  must  be  investigated.  In  order  to 
provide  a  worked  example  of  the  method  some  details 
of  the  calculation  will  be  given. 

Mode  I  collapse 

Beam-column  parameters: 


CT/r  =  — - —  (Typ  =  262  N/mm- 

ffom  (8.3.6):  Ci  =  .096,  /„  =  352000  mm\ 
p  =  15.73  mm 
C:  =  0.854,  yp,r  =  7.41  mm 

from  (14.2.18):  Ap  = 

=  (29. 5)297 

=  0.958  mm 

I/A  p-  (15.73)- 
from  (14.2.21):  A.  =  ^ 

=  0.55  mm 

from  (14.2.16):  A  =  0.63,  -q  =  0.060 
from  (14.2.19):  Pp  =  .0287 
from  (14.2.20):  Rn  =  0.899 
from  (14.2.22):  (Ta,u  =  202.3  N/mm- 
andfrom  (14.3.3):  =  1.128  X  10"’ 


from  (8.3.6):  C,  =  0.083,  I  =  326700  mm\ 
p  =  14.65  mm 

Cz  =  0.8775,  Vf  =  -44.6  mm 
from  (14.2.3):  A  =  0.623,  p  =  0.052 
from  (14.2.5):  R  =  0.925 

and  hence  the  value  of  cra.„  for  a  Mode  I  collapse  is 
0.925cr|.,  =  255.8  N/mm-.  From  (14.3.1)  the  value 
oiSa^  is  1.226  x  l0-\ 

Mode  11  collapse 

Plate  slenderness:  /3  =  1.89 
from  (14.2.8):  T  =  0.825, 


TABLE  14.2  PROPERTIES  OF  BOX  GIRDER  MODEL 


Component 

Dimensions 

(mm) 

cry 

(N/mra*) 

E 

(N/mm*) 

Compression 

4.88 

298.0 

208,500 

flange 

Tension 

4.88 

298.0 

208,500 

flange 

Web 

3:38 

211.6 

216,200 

Stiffeners 

50.8  X  15.9  X  4.8L 

276.5 

191,500 

Since  the  Mode  1  value  of  (Xa.u  was  larger  than  this 
Mode  II  value  the  latter  is  selected  as  the  predicted 
value  of  ultimate  strength. 

As  shown  in  Fig.  14.13  this  predicted  result  is  in 
quite  good  agreement  with  the  experimental  results: 


(N/mm^) 


Figure  14.13  Collapse  of  upper  tiange  of  box  girder  model. 
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Mode  n  collapse  a:  an  average  applied  stress  of  205.3 
N/mror  and  a  ‘•train  of  1.070  x  10*■^ 

It  should  be  noted  that  the  selection  of  the  predicted 
mode  of  collapse  must  be  based  on  cra,u  and  not  on  £a,u^ 
A  smaller  value  of  ea,u  for  Mode  I  does  not  mean  that 
Mode  I  collapse  would  occur  fint,  because  the  Mode 
n  collapse  actually  begins  much  earlier,  when  the 
plate  begins  to  fail,  and  from  this  point  onward  the 
beam-column  undergoes  ever-increasing  positive 
bending,  opposite  to  that  of  Mode  1.  Thus  a  Mode  I 
collapse  can  only  occur  if  the  value  of  {(Ta,u)\  is  smaller 
than  the  load  which  initiates  plate  failure.  Hence  we 
are  being  conserv^ative  in  saying  that  a  Mode  I  collapse 
only  requires  that  (cr^. be  less  than  and  then 

taking  this  smaller  value  to  be  the  ultimate  strength  of 
the  panel.  But  as  mentioned  earlier,  whenever  the  lat¬ 
eral  load  is  small  or  is  absent  the  occurrence  of  Mode 
I  or  Mode  n  depends  not  only  on  the  properties  and 
proportions  of  the  panel  but  also  on  the  pardcular 
circumstances  and  details  such  as  initial  deflections, 
load  sequence,  and  residual  stress,  which  cannot  be 
known  or  predicted  by  the  designer.  Hence  the  de¬ 
signer  must  always  be  conservative  and  assume  the 
worst,  which  means  choosing  whichever  mode  gives 
the  lower  value  of  strength. 


Multispan  Panels  With  Lateral  Load 

We  next  compare  the  method's  predictions  with  the 
results  of  tests  by  Smith  [6]  on  multispan  panels  sub¬ 
jected  to  a  uniform  lateral  pressure  in  addition  to  axial 
compression.  In  comparing  theory  and  experiment  it  is 
necessary^  to  have  full  information  about  the  initial 
conditions  (residual  stress,  preload  deflections,  and 
the  values  of  5o  and  Mq  due  to  the  lateral  load)  for  the 
particular  stiffener-plate  combination  which  triggered 
the  collapse.  Of  the  11  panels  tested  in  Ref.  6  panels 
3a  and  3b  are  chose  for  the  comparison  because  for 
these  panels  aU  of  the  required  information  was 
recorded  and  presented.  Also  these  two  panels  provide 
valuable  information  about  the  effect  of  lateral  pres¬ 
sure  because  they  were  a  matching  pair  and  one  of 
them  (3a)  had  a  lateral  pressure  of  3  psi  whereas  the 
other  had  no  lateral  load. 

Panel  3a 

This  panel  was  subjected  to  a  constant  lateral  pressure 
of  3  psi.  The  panel  collapsed  when  the  applied  axial 
stress  reached  0.69  ay,  where  ay  is  the  mean  value  of 


yield  stress  of  the  stiffeners  and  plating,  averaged  over 
the  panel  cross  section. 

This  case  is  worth  examining  in  some  detail,  es¬ 
pecially  in  regard  to  the  causal  sequence  of  failures, 
because  it  provides  an  example  of  how  a  laterally 
loaded  multispan  panel,  in  spite  of  the  rotational  re¬ 
straint  at  the  ends  of  each  span,  can  undergo  the  same 
type  of  collapse  as  a  laterally  loaded  pinned  panel. 

As  shown  in  Fig.  14.14  the  collapse  involved  up¬ 
ward  flexural  buckling  in  one  of  the  central  bays  and 
downward  bucking  of  an  adjacent  bay,  together  with  a 
stiffener  flange  mechanism  of  the  type  shown  in 
Fig.  14.4,  close  to  the  central  transverse  frame.  All 
three  failures  occurred  more  or  less  simultaneously, 
and  to  determine  the  precise  interactions  and  the  causal 
relationships  we  must  examine  the  state  of  loading  just 
prior  to  failure.  Figure  14.15  illustrates  the  bending 
moment  distribution  in  the  stiffeners  due  to  the  lateral 
pressure.  This  was  not  given  in  Ref.  6  but  was  calcu¬ 
lated  by  the  author  by  a  nonlinear  finite  element  anal¬ 
ysis  of  the  gross  panel.  The  total  compressive  stress  in 
the  stiffener  flange  is  a  maximum  at  the  midlength  of 
the  gross  panel,  because  of  the  peak  value  of  bending 
moment  there.  The  finite  element  analysis  showed  that 
when  the  applied  axial  stress  reached  the  experimental 
collapse  value  the  total  stress  in  the  stiffener  flange  at 
the  midlength  had  just  slightly  exceeded  the  yield 
stress  of  the  stiffener.  At  this  time  yield  had  not  oc¬ 
curred  anywhere  else  and  all  of  the  panels  were  still 
stable.  This  shows  clearly  that  it  was  the  stiffener 


Figure  14.14  Panel  3a  after  collapse.  (Photo  counesy  of  Admi¬ 
ralty  Marine  Technology  Establishment.  Dunfermline.  Scotland). 
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14.4  COMPARISON  OF  THEORY  .AND  EXPERLMENT 


Figure  14.15  Bending  moment  in  panel  3a. 


flange  mechanism  which  triggered  the  collapse,  and 
also  provides  the  explanation  as  to  why  the  panels 
buckled  at  the  same  time  as  the  formation  of  the  stiff¬ 
ener  flange  mechanism.  Such  a  mechanism  is  virtually 
the  same  as  a  hinge  (an  ordinary  hinge,  not  a  plastic 
hinge)  because  the  flange  has  lost  nearly  all  of  its  axial 
stiffness.  Therefore  the  stiffeners  rotated  about  this 
point  and  transferred  the  bending  moment  into  the 
spans.  The  new  bending  moment  in  each  span  is  ap¬ 
proximately  that  for  a  beam  which  is  pinned  at  one  end 
and  clamped  at  the  other,  with  a  midspan  value  of 
pba}l\l.  The  increased  bending  stress  in  the  plating 
was  sufficient  to  cause  a  Mode  II  collapse  of  span  AB, 
and  so  this  span  collapsed  immediately  after  the  flange 
niechanism  had  occurred,  with  no  further  increase  in 
axial  load.  Since  the  transverse  frame  provides  very 
Uttle  rotational  restraint  and  also  permits  at  least  a 
small  lateral  deflection,  it  is  the  hinge  rather  than  the 
frame  which  constitutes  the  endpoint  of  the  two  beam- 
columns.  and  since  AB  is  longer  than  BC  the  Mode  II 
collapse  occurred  in  the  former. 

This  case  demonstrates  that  in  spite  of  the  rotational 
restraint  which  exists  at  the  ends  of  each  span  in  a 
laterally  loaded  multispan  panel,  collapse  can  still  oc¬ 
cur  in  a  manner  very  similar  to  that  of  a  panel  with 
pinned  ends.  Therefore  the  pinned  beam-column  ap¬ 
proach,  with  allowance  for  the  values  of  and  Sj 
which  exist  within  the  span,  is  directly  relevant  to  such 
panels. 

In  the  rationally-based  design  process  of  Fig.  2. 19, 
the  hull  module  analysis,  member  limit  analysis  and 
constraint  evaluation  of  the  panels  (steps  3,  4,  and  5) 
will  indicate  if  a  stiffener  flange  mechanism  can  occur, 
and  if  so  the  use  of  any  reasonable  estimates  of  the 
post-mechanism  values  of  Mo  and  5o  will  give  a  good 
estimate  of  cra  „.  For  example,  using  the  approximate 
post-mechanism  values  of  pba~j\2  for  Mo  and 
T  (5pb<2'‘)/384£'/  (half  of  the  simply  supported  value) 
for  5o,  together  with  a  typical  design  value  of  a/750 


for  A,  gives  a  value  of  0.658a-  for  which  differs 
from  the  acmal  value  by  only  5%,  and  on  the  conser¬ 
vative  side. 

The  acmal  measured  value  of  A  for  the  span  where 
the  collapse  occurred  was  A  =  0.03  in.  (a/2000)  and 
when  this  value  is  used  the  predicted  value  of  cra  „  is 
0.69^y,  the  same  as  the  experimental  value. 

Panel  3b 

A  second  panel  which  was  virtually  identical  to  panel 
3a  collapsed  at  an  applied  stress  of  =  0.6  IcT).. 
Fig.  14.16  shows  the  panel  after  collapse.  Since  there 
is  no  lateral  load,  a  Mode  II  collapse  would  also  have 
been  possible  and  one  of  the  requirements  of  any 
method  for  estimating  ultimate  strength  is  that  it 
should  correctly  identify  the  critical  mode,  that  is,  the 
mode  with  the  lower  value  of  cr^„.  For  this  panel  the 
initial  deflection  of  the  collapsed  span  was  approxi¬ 
mately  0.15  in.  toward  the  plating  (A  =  -a/400). 
With  the  method  presented  herein  the  predicted  values 
of  for  Modes  I  and  II  are  Q.SlWy  and  0.69^ 
Thus  the  method  correctly  identifies  the  coOapse  mode 
and  predicts  a  value  for  which  difiers  trom  the 
experimental  value  by  only  7%  on  the  conservative 
side.  This  1%  margin  comes  mainly  from  the  con¬ 
servative  assumption  that  a  stiffener  flange  mechanism 
forms  as  soon  as  the  stress  in  the  midthickness  of  the 
stiffener  flange  reaches  yield,  whereas  in  reality  it  forms 
only  after  the  flange  is  fully  yielded. 


Figure  14.16  Panel  3b  after  collapse.  (Photo  courtesy  of  Admi¬ 
ralty  Marine  Technology  Establishment,  Dunfermline,  Scotland). 


•The  latter  is  obtained  using  the  largest  positive  value  of  A.  which 
was  0.22  in.  ((i/270). 
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ULTIMATE  STRENGTH  OF  STIFFENED  PANELS 


In  the  previous  multispan  examples  the  collapse  has 
been  Mo^  H,  and  hence  this  example  demonstrates 
the  accuracy  of  the  method  for  a  Mode  I  collapse  of  a 
multispan  panel. 

Conclusions 

The  analysis  of  panel  3b  shows  that  when  a  multispan 
panel  has  no  lateral  load  the  ultimate  strength  is  basic- 
aUy  that  of  a  simply  supported  subpanel.  For  reference 
purposes  we  shall  denote  this  value  as  (cra,«)o.  It  is 
always  the  lesser  of  two  values:  either  the  Mode  I 
value  from  (14.2.6)  or  the  Mode  n  value  from 
(14.2.22).  If  the  lateral  load  is  above  a  certain  mag¬ 
nitude,  each  subpanel  is  effectively  clamped  and  if  this 
remains  true  then  the  ultimate  strength  is  larger  than 
for  a  simply  supported  panel.  For  panels  3a  and  3b  the 
increase  was  from  0.61^^  to  0.69^^',  or  13%.  But  it  is 
extremely  difficult  to  predict  the  magnitude  of  lateral 
load  that  is  required  to  constrain  all  subpanels  to  be¬ 
have  as  clamped,  first,  because  the  behavior  of  each 
subpanel  is  strongly  influenced  by  local  factors  (eccen¬ 
tricity,  residual  stress,  nonuniformity  of  lateral  load, 
etc.)  and,  second,  because  each  subpanel  interacts 
with  its  neighbors.  Moreover,  although  the  lateral  load 
raav  be  sufficient  to  keep  each  subpanel  rotationally 
restrained  for  a  time,  such  that  the  value  of  cTc  can 
exceed  icra,u)o.  there  is  always  the  possibility  of  a  stiff¬ 
ener  flange  mechanism,  whereupon  two  or  more  of  the 
subpanels  suddenly  change  to  the  alternating  wave 
shape  and  immediately  undergo  collapse,  because  (Ja 
already  exceeds  (cra.u)o-  We  thus  have  four  factors  to 
consider: 

1.  The  possibility  that  the  lateral  load  may  become 
small. 

2.  The  difficulty  of  calculating  the  magnitude  of  the 
lateral  load  which  would  guarantee  end  clamping. 

3.  The  possibility  of  a  stiffener  flange  mechanism, 
which  becomes  more  likely  with  increasing  lateral 
load. 

4.  The  undesirability  of  sudden  collapse. 

For  all  of  these  reasons,  any  value  of  ultimate  strength 
larger  than  (cTo.Jo  should  not  be  used  for  design  pur¬ 
poses,  even  though  larger  values  may  in  fact  be 
achieved  when  there  is  a  lateral  load.  That  is,  (orfl.u)o  is 
an  upper  limit  or  “cutoff’  value.  In  a  rationally-based 
design  process  it  is  necessar}'  to  calculate  cr^.«  for  each 
load  case.  For  the  various  load  cases  involving  a  lateral 
load  it  would  be  pKDSsible  to  calculate  accurate  values 
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of  (7a,u  by  using  the  information  from  steps  3,  4,  and 
5  of  the  overall  design  procedure,  as  discussed  pre¬ 
viously.  But  for  design  purposes  this  is  not  necessary 
because  these  values  will  nearly  always  exceed  the 
cutoff  value,  (cr^.u)o.  Hence  for  design  purposes  it  is 
sufficient  to  use  the  simple  and  conservative  pinned 
beam-column  approach,  with  Mo  and  6o  equal  to  their 
simply  supported  values,  for  estimating  the  values  of 
for  the  various  load  cases  involving  a  lateral  load. 
Although  the  true  ultimate  strength  would  be  larger 
than  this  value  it  would  usually  be  ruled  out  because  it 
exceeded  (cr^.Jo-  Therefore  the  value  obtained  by  the 
pinned  beam-column  approach  will  not  be  unduly  con¬ 
servative.  For  example,  for  the  load  case  involving 
lateral  pressure  (panel  3a)  the  use  of  the  simply  sup- 
pxmed  values  of  Mo  and  6o  and  a  typical  value  of 
A  (fl/750)  gives  a  value  of  0.56o^y.  For  this  load  case 
the  true  ultimate  strength  is  0.69cry,  but  this  cannot  be 
used  for  design  purposes  because  it  exceeds  the  value 
for  zero  lateral  load,  (cTa.Jo  ~  0.61cry.  Hence,  taking 
O.SSay  as  the  ultimate  strength  for  the  laterally  loaded 
case  is  not  overly  conservative  since  it  is  only  8% 
below  the  experimentally  obtained  value  of  (cr^.Jo- 
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